Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

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

Search the Archive

Re: How to rectify the error for NDSolve ?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg127158] Re: How to rectify the error for NDSolve ?
  • From: Rahul Chakraborty <rahul.6sept at gmail.com>
  • Date: Mon, 2 Jul 2012 22:19:27 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com

Dear all,

 Thanks to Murray for your reply. I have made changes as per advice but still no result. The differential eqn is also changed to make it simpler. 

 Here is my code:

 Clear[x,\[Tau],Q,\[CapitalGamma],\[Phi]];
Q=100;
\[CapitalGamma]=50;
Subscript[\[Omega], 0]=10;
Subscript[V, p]=0.5;
Subscript[C, 0]=4.06 * 10^-4;
\[Phi]=90;
Subscript[V, 1](t)=Sqrt[2] Subscript[V, p]Sin(2Subscript[\[Omega], 0]t);
eqn= x''[t]+Subscript[\[Omega], 0]/Q x'[t]+Subscript[(\[Omega]^2), 0][1+Subscript[c, 0] Subscript[V, 1](t)]x[t]==\[CapitalGamma]Sin[Subscript[\[Omega], 0]t+\[Phi]];
sol=NDSolve[{eqn,x[0]==0.5,x'[0]==0},x[t],{t,0,250},MaxSteps-> 500000,WorkingPrecision-> 25][[1]]
ParametricPlot[Evaluate[{x[t],x'[t]}/.sol],{t,0,100},Frame->True,PlotRange-> All,AxesLabel-> {"t","Overscript[\[Theta], .]"}]

Regards,
rc



  • Prev by Date: Re: Rearranging terms - user-defined
  • Next by Date: Re: Numerical integration over an arbitrary 2D domain
  • Previous by thread: Re: How to rectify the error for NDSolve ?
  • Next by thread: Re: How to rectify the error for NDSolve ?