NDSolve::ndsz question
- To: mathgroup at smc.vnet.net
- Subject: [mg63599] NDSolve::ndsz question
- From: "Flip" <fliptomato at gmail.com>
- Date: Sat, 7 Jan 2006 02:29:23 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
Greetings, I'm trying to solve a differential equation but I'm getting
the following error:
NDSolve::ndsz: At x == 18.84025621252942`, step size is effectively
zero; \
singularity or stiff system suspected
I'm not sure how to deal with this--is there a way to allow smaller
step sizes? Here is my complete code:
k := 400000000;
b := 0.0023;
epsR := 0.024;
solution := NDSolve[
{y'[x] == -k x^(-0.5) \[ExponentialE]^(-epsR x) (y[x]^2 - (b x^2
BesselK[2,x])^2),
y[0.1] == b*BesselK[2,0.1]/(0.1)}, y, {x,0.1,5000}]
Evaluate[y[5000] /. solution]
Any thoughts would be appreciated! Thanks,
Flip
- Follow-Ups:
- Re: NDSolve::ndsz question
- From: Pratik Desai <pdesai1@umbc.edu>
- Re: NDSolve::ndsz question