I'm solving a system of odes present in function1 for different set of initial conditions.
for i = 1:10
x0 = rand(10,1)
[t x] = ode15s(@(t,s) function1(t,s), tspan , x0);
field = "x" + num2str(i);
soln.(field) = x
I'd like to run these simulations parallely, i.e. for all i's at once.
Could someone suggest how this can be done?