       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 := 1
sinc[x_?NumericQ] := Sin[x]/x

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

and

NDSolve[{deqn2, y == 0, y' == 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, y'==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?