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.math.mathematica

Topic: System of differential equations
Replies: 3   Last Post: Dec 7, 2012 1:37 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
hazem.abdelhafiz@gmail.com

Posts: 2
Registered: 12/6/12
Re: System of differential equations
Posted: Dec 7, 2012 1:37 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


Thanks Bob,
Thanks for the fast comment
It now solves


On Fri, Dec 7, 2012 at 1:34 AM, Bob Hanlon <hanlonr357@gmail.com> wrote:

> The Sin functions must be entered as Sin vice sin. The range of t is
> too large to be handled.
>
> s = NDSolve[{
> ph1''[t] + ph1'[t] + Sin[ph1[t]] - mu[t] - mu'[t] == 0.005,
> ph2''[t] + ph2'[t] + Sin[ph2[t]] + mu[t] + mu'[t] == 0.005,
> 0.2 mu'[t] + mu[t] == ph1'[t] - ph2'[t], ph2[0] == 0, ph1[0] == 0,
> ph2'[0] == 0, ph1'[0] == 0, mu[0] == 0},
> {ph1[t], ph2[t], mu[t]},
> {t, 0, 240}];
>
> Plot[Evaluate[Tooltip[#, ToString[#]] & /@
> {ph1[t], ph2[t], mu[t]} /. s],
> {t, 0, 10},
> Frame -> True, Axes -> False]
>
> However, the range on t affects the results
>
> s = NDSolve[{
> ph1''[t] + ph1'[t] + Sin[ph1[t]] - mu[t] - mu'[t] == 0.005,
> ph2''[t] + ph2'[t] + Sin[ph2[t]] + mu[t] + mu'[t] == 0.005,
> 0.2 mu'[t] + mu[t] == ph1'[t] - ph2'[t], ph2[0] == 0, ph1[0] == 0,
> ph2'[0] == 0, ph1'[0] == 0, mu[0] == 0},
> {ph1[t], ph2[t], mu[t]},
> {t, 0, 100}];
>
> Plot[Evaluate[Tooltip[#, ToString[#]] & /@
> {ph1[t], ph2[t], mu[t]} /. s],
> {t, 0, 10},
> Frame -> True, Axes -> False]
>
>
> Bob Hanlon
>
>
> On Thu, Dec 6, 2012 at 4:58 AM, <hazem.abdelhafiz@gmail.com> wrote:

> > Dear all,
> > I want to solve this system of equations but I always get an error
> >
> > NDSolve::ndnum:Encountered non-numerical value for a derivative at t==0
> >
> > can any one help me on this
> >
> > s=NDSolve[{
> > ph1''[t]+ph1'[t]+sin[ph1[t]]-mu[t]-mu'[t]==0.005,
> > ph2''[t]+ph2'[t]+sin[ph2[t]]+mu[t]+mu'[t]==0.005,
> > 0.2mu'[t]+mu[t]==ph1'[t]-ph2'[t],
> > ph2[0]==0,
> > ph1[0]==0,
> > ph2'[0]==0,
> > ph1'[0]==0,
> > mu[0]==0
> > },
> > {ph1[t],ph2[t],mu[t]},
> > {t,0,5000}
> > ]
> >
> > thanks in advance,
> >

>





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.