MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Error with NDSolve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg50949] Re: Error with NDSolve
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Wed, 29 Sep 2004 03:15:15 -0400 (EDT)
  • Organization: Universitaet Leipzig
  • References: <cjaru4$nni$1@smc.vnet.net>
  • Reply-to: kuska at informatik.uni-leipzig.de
  • Sender: owner-wri-mathgroup at wolfram.com

Hi,

sol = NDSolve[{si2'[t] == 2, si2[0] == 0, th2'[t] == Sin[si2[t]]/r2[t], 
        th2[0] == 0, r2'[t] == Cos[si2[t]], r2[0] == 1}, {si2, r2, th2},
{t, 
        0, 25}];

r[t_] = r2[t] /. First[sol];
th[t_] = th2[t] /. First[sol];
si[t_] = si2[t] /. First[sol];

no error. Avoid the usage of % in the FrontEnd.

Regards
  Jens

"Narasimham G.L." wrote:
> 
> There is error in si[t] even after the interpolating function is
> obtained. Thanks for fix.
> 
> NDSolve[  {si2'[t] == 2, si2[0] == 0,
> th2'[t] == Sin[si2[t]]/r2[t], th2[0] == 0, r2'[t] == Cos[si2[t]],
> r2[0] == 1}, {si2, r2, th2}, {t, 0, 25}]; r [t_] = r2[t] /. First[%];
> th[t_] = th2[t] /. First[%];
> (* error -> *)  si[t_] = si2[t] /. First[%];


  • Prev by Date: Re: Hyperbolic function identity
  • Next by Date: Re: Error with NDSolve
  • Previous by thread: Re: Error with NDSolve
  • Next by thread: Re: Error with NDSolve