Services & ResourcesWolfram ForumsMathGroup Archive

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

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

• Follow-Ups:
  • Re: How to rectify the error for NDSolve ?
    • From: Bob Hanlon <hanlonr357@gmail.com>