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]