Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'parameter scan for NDSolve' topicStudent Support Forum > General > Archives > "parameter scan for NDSolve"

Next Comment >Help | Reply To Topic
Author Comment/Response
Steve
08/06/12 05:14am

Hi,

I have a system of ordinary differential equations set up in a function so they can take their parameter values from an array e.g {k1,k2,k3}.

I then plot the output of NDSolve for the equations (a time course simulation).

I would like to be able to "scan" the parameter values and observe the steady state value of the ordinary differential equations.

At the moment, I have something like this:

parameters={k1,k2,k3}

sim[X_,parameters_,tmax_:72[:=
solution=NDSolve[equation1'[t]= ..... , equation2'[t]= ........, equation1[0],equation2[0]/.parameters,{equation1,equation2},{t,0,tmax}]

I would like to be able to input a range of values for the parameters k1,k2,k3 and plot the steady state of equation[1] and equation[2] in response to the changes.

Is the possible in Mathematica? I was wondering about creating some kind of heat map, but I am very new to Mathematica and struggling with the code.

Thanks for any assistance.

URL: ,

Subject (listing for 'parameter scan for NDSolve')
Author Date Posted
parameter scan for NDSolve Steve 08/06/12 05:14am
Re: parameter scan for NDSolve Bill Simpson 08/07/12 5:36pm
Next Comment >Help | Reply To Topic