Re: NDSolve::ndsz question
- To: mathgroup at smc.vnet.net
- Subject: [mg63719] Re: NDSolve::ndsz question
- From: "Flip" <fliptomato at gmail.com>
- Date: Wed, 11 Jan 2006 06:49:42 -0500 (EST)
- References: <dpnqoa$6rh$1@smc.vnet.net><dpvmjv$ltf$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Pratik, Jean, Brian, and Jens--thanks to all for your replies. They've
been very helpful. I was wondering if I could impose and ask a slightly
more general question:
My goal in doing this calculation was, given this ODE, to determine the
behavior of y[x] as one varies k.
Thus, given this ODE:
Derivative[1][y][x] == (-k)*(1/(Sqrt[x]*E^(epsR*x)))*
(y[x]^2 - (b*x^2*BesselK[2, x])^2)
b = 23/10000
epsR = 3/125
k = variable
With initial condition: (slightly different from previous post):
y[1/10] == b*(1/10)^2*BesselK[2, 1/10]
And a condition:
y[infinity] = some number (say, 10^-10)
How do I determine which k yields this value? Brian seems to have
solved the ODE analytically (which I'm curious about, though it is
somewhat outside the scope of this forum), is there a clever way to get
mathematica to do this computaiton?
Thanks for all of your help,
Flip