Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Topic: (yet unsolved) how to Solve coupled ODEs
Replies: 4   Last Post: Feb 13, 2013 4:11 AM

 Messages: [ Previous | Next ]
 Torsten Posts: 1,717 Registered: 11/8/10
Re: (yet unsolved) how to Solve coupled ODEs
Posted: Feb 13, 2013 4:11 AM

"Henry" wrote in message <kfeavo\$2f\$1@newscl01ah.mathworks.com>...
> I modified the code to:
>
> tspan = linspace(1, 2,100);
> y = 0.0001*ones(1,200);
> [T,sols] = ode45(@yprime, tspan,y);
> and
> dyds = [dyds1' ; dyds2'];
>
> and it now runs, thanks.
>
> But I have this question: My odes are time-independent, function of "s" which is a distance, but since the ode45 will not work without a tspan, I came up with an arbitrary tspan. It looks as if ode45 interprets that "s" is a time? I'm confused as to how interpret the result of the routine which has dimensions of 100x200. I was expecting 1 x 200, i.e. independent of time.
> thanks

By the way:
I suspect that you want to set different values of a,b,c and d for different values of s and solve for only _two_ solutions y(1) and y(2). That's not what you do in your code above.
You will have to use MATLAB's interp1 function to get values for a,b,c and d at the distance s where the solver needs the values of the derivatives dyds(1) and dyds(2).

Or do you want to make a parameter study, i.e. different solutions for different
combinations of the parameters a,b,c and d ? Then your code above
is also incorrect because you always refer to y(1) and y(2) in the evaluation of the derivatives, thus the solution for only the _first_ parameter set.