MathGroup Archive 1999

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

Search the Archive

Re: How to NDSolve the differential equation

  • To: mathgroup at smc.vnet.net
  • Subject: [mg19640] Re: How to NDSolve the differential equation
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Tue, 7 Sep 1999 00:28:36 -0400
  • Organization: Universitaet Leipzig
  • References: <7qulo3$2r5@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Hi Chee,

your equation is singular at the origin. How ever a convergent power
series 
exist and a numerical solution is possible.

The problem ov NDSolve[] ist that it divedes the equation by x and get
Sin[x]/x but it does not calculate a Limit[] at x->0 it just find a
(not existing) singulaity. You cando the following:
Define the Sin[x]/x function


deqn = x y''[x] + Sin[x] y'[x] + 2 x y[x] == 0

sinc[0.] := 1
sinc[0] := 1
sinc[x_?NumericQ] := Sin[x]/x

deqn2 = (Simplify[#/x & /@ deqn]) /. Sin[x]/x -> sinc[x]

and 

NDSolve[{deqn2, y[0] == 0, y'[0] == 1}, y[x], {x, 0, 2}]

works as expected.

Hope that helps
  Jens


Chee Lim Cheung wrote:
> 
> Dear Mathematica gurus & users,
> 
> I encountered error messages with the tag Power::infy when I tried to do
> the following:
> 
> NDSolve[{x y''[x] + Sin[x] y'[x] + 2 x y[x] == 0,y[0]==0, y'[0]==1}, y,
> {x,0,2}]
> 
> Can anyone help me in getting Mathematica to produce an answer? I am using
> Mathematica 4.0 & Mathematica 3.0.
> 
> Thanks
> Chee


  • Prev by Date: Re: Nonlinear Fit with complex model
  • Next by Date: Re: Re: Double evaluation in recursive function
  • Previous by thread: Re: How to NDSolve the differential equation
  • Next by thread: crystal structure?