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

I have made the changes as suggested by Bob but still some errors are there:

 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;
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], 0]^2)[1+Subscript[c, 0]Subscript[V, 1][t]]x[t]==\[CapitalGamma] Sin[Subscript[\[Omega], 0]t+\[Phi]]//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]

Errors:

 Set::write: Tag Times in (0.707107 Sin[20 t])[t] is Protected. >>

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



Plz suggest what shall i do to get the solution.


Regards,
rc



  • Prev by Date: Re: how to export to Excel in Mathematica from a Parametricplot command
  • Next by Date: Re: how to export to Excel in Mathematica from a Parametricplot command
  • Previous by thread: Re: Re: How to rectify the error for NDSolve ?
  • Next by thread: Re: Default path for opening files