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: [mg127172] Re: How to rectify the error for NDSolve ?
• From: Rahul Chakraborty <rahul.6sept at gmail.com>
• Date: Thu, 5 Jul 2012 06:11:34 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com

Following is the code as modified:

 Clear[x,Q,\[CapitalGamma],\[Phi]];
Q=100;
\[CapitalGamma]=50;
Subscript[\[Omega], 0]=10;
Subscript[V, p]=1/2;
Subscript[C, 0]=406 * 10^-6;
\[Phi]=90;
\[Gamma]=1/Q;
Subscript[V, 1][t]=Sqrt[2] Subscript[V, p]Sin[2Subscript[\[Omega], 0]t];
eqn= x''[t]+Subscript[\[Gamma]\[Omega], 0]x'[t]+(Subscript[\[Omega], 0]^2)[1+Subscript[c, 0]Subscript[V, 1][t]]x[t]-\[CapitalGamma] Sin[Subscript[\[Omega], 0]t+\[Phi]]==0//Simplify;
sol=NDSolve[{eqn,x[0]==1/2,x'[0]==0},x[t],{t,0,25}][[1]]
ParametricPlot[Evaluate[{x[t]/.sol,D[x[t]/.sol,t]}],{t,0,25},Frame->True,AxesLabel->{"t",Overscript[\[Theta],"."]},AspectRatio->1]


 ERROR showing as follows:

 NDSolve::ndnum: Encountered non-numerical value for a derivative at t == 0.`. >>

 Plz advice me how to correct this error.
  
Regards,

 rc

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