MathGroup Archive 2007

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

Search the Archive

Re: ndsolve ndsv problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg79853] Re: ndsolve ndsv problem
  • From: chuck009 <dmilioto at comcast.com>
  • Date: Tue, 7 Aug 2007 01:25:40 -0400 (EDT)

I don't see where you're supplying the end points for the numeric calculations.  The code I wrote below just sets the boundary point at infinity at 1000 and does generate a solution which matches the two boundary points you supplied somewhat ok:

sol = NDSolve[0.0005635632304299039*psi[R] == 
     (2*Derivative[1][psi][R])/R + Derivative[2][psi][R] && 
    Derivative[1][psi][1] == 7.017487714753473 && 
    psi[1000] == 0, psi, {R, 1, 1000}, WorkingPrecision -> 30]
f[x_] := Evaluate[psi[x] /. Flatten[sol]]; 
Plot[f[x] /. Evaluate[Flatten[sol]], {x, 1, 100}]

In[146]:=
f[1000]
fd[x_] = D[f[x], x]
fd[1]

Out[146]=
\!\(0``40.78712463722495\)

Out[147]=
InterpolatingFunction[{{1.00000000000000000000000000000,1000.\
00000000000000000000000000}},<>][x]

Out[148]=
7.0324851897180024334773666



> Steps to reproduce:
> 1. Issue this command (given to you in InputFrom due
> to a copy/paste problem)
> 
> NDSolve[0.0005635632304299039*psi[R] ==
>    (2*Derivative[1][psi][R])/R +
>     Derivative[2][psi][R] &&
>   Derivative[1][psi][1] == 7.017487714753473 &&
>   psi[Infinity] == 0, psi, R]
>


  • Prev by Date: Projectile motion
  • Next by Date: Re: Solving differential equations with parameters
  • Previous by thread: ndsolve ndsv problem
  • Next by thread: Using ShowLegend