Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

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?


  • Prev by Date: Apply function to parts of a list
  • Next by Date: Rotation of 3D objects
  • Previous by thread: Re: How do I reformulate my NDSolve program
  • Next by thread: OptimizeExpression