NDSolve question
- To: mathgroup at smc.vnet.net
- Subject: [mg44408] NDSolve question
- From: "Emmanuel Dechenaux" <dechenau at purdue.edu>
- Date: Sat, 8 Nov 2003 04:50:51 -0500 (EST)
- Organization: Purdue University
- Sender: owner-wri-mathgroup at wolfram.com
Hi,
I have what probably is a very basic question. I'm trying to solve the
complicated first order differential equation below. There is a singularity
at zero (which is my initial condition). The equation is also highly
non-linear. I'm not quite sure what to do. I used different initial
conditions as well, but Mathematica stops very quickly at x=3...... even if
I increase MaxSteps. Again, I'm not quite sure what to do. If someone could
help me, that would be great. Thanks.
Emmanuel
(My apology if the mathematica input didn't paste cleanly).
r = 0.7
solution =
NDSolve[{ (y' )[x] == ( ((x - y[x]) )^ ((r) ) - (( (-y[x] )) )^ ((r) ) ) /
(x* ((r) ) ((x - y[x]) )^ ((r - 1) ) + ((1000 - x) ) (( (-y[x] )) )^ ((r -
1) ) ), y[0] == 0}, y, {x, 0, 1000}, MaxSteps -> 100000]