Re: How do I reformulate my NDSolve program
- To: mathgroup at smc.vnet.net
- Subject: [mg86217] Re: How do I reformulate my NDSolve program
- From: Alex Cloninger <acloninger at wustl.edu>
- Date: Thu, 6 Mar 2008 02:58:31 -0500 (EST)
Sorry. It was supposed to print epsilon, my general parameter, but I guess it didn't carry over. You could just say instead
solution = NDSolve[{x'[t] == 2p[t], x[0] == 0,
p'[t] == I*(2 + 1/2)(I*x[t])^(1 + 1/2), p[0] == 1}, {x, p}, {t,0,10}, WorkingPrecision -> 40, MaxSteps ->Infinity][[1]];
And:
ParametricPlot[{Re[p[t]] /. solution, Im[p[t]] /. solution}, {t,0,10},PlotRange -> {{-2, 2}, {-2, 2}}]
Sorry about that. Could you help me out now?