Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

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


Math Forum » Discussions » Software » comp.soft-sys.matlab

Topic: how to Solve coupled ODEs
Replies: 5   Last Post: Feb 12, 2013 11:07 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Torsten

Posts: 1,450
Registered: 11/8/10
Re: how to Solve coupled ODEs
Posted: Feb 7, 2013 3:15 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"Henry" wrote in message <keujnd$h64$1@newscl01ah.mathworks.com>...
> Hi there,
>
> I'm having a hard time solving a set of 2 coupled ODEs:
>
> dx/ds = 1/x * [y* (d + y/s) - a*x*f]
> dy/ds = 1/x * [-y * (b + y) * f] - y/s - c
>
> where f = sqrt(x^2 + y^2*(e + y)), and all constants are vectors of (100 x 1). The independent variable is "s".
>
> %----------------------------------------------------------------------
> global a b c d e
> y = ode45(@odeeqns, s, [1 1])
>
> function dyds = odeeqns(s,y)
> global a b c d e
> dyds(:,1) = 1./y(1) .*(y(2) .* (d + s ./ y(2)) - a .* y(1) .* sqrt(y(1).^2 + (e + y(2)).^2));
> dyds(:,2) = - (y(2) .* (b + y(2)) .* sqrt(y(1).^2 + (e + y(2)).^2))./y(1) - y(2)./s - c;
> return
> %----------------------------------------------------------------------
>
> Any help please?
> Thanks!


Take a look at example 1 under
http://www.mathworks.de/de/help/matlab/ref/ode23.html
and you'll recognize the two (or three) errors you made in your small piece of code
from above.

Best wishes
Torsten.



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.