Services & Resources / Wolfram Forums / MathGroup Archive

MathGroup Archive 2013

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

Search the Archive

Problem in solving Differential Equation

  • To: mathgroup at smc.vnet.net
  • Subject: [mg120032] Problem in solving Differential Equation
  • From: Rahul Chakraborty <rahul.6sept at gmail.com>
  • Date: Tue, 5 Mar 2013 22:14:41 -0500 (EST)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-newout@smc.vnet.net
  • Delivered-to: mathgroup-newsend@smc.vnet.net
Dear all,

 Following differential equation seems to have some error in coding by me. kindly let me know where i have gone wrong.

Clear[x];
 \[Omega]:=-2;
eqn=x''[t]+ Sin  x[t]-\[Omega]^2  Sin  x[t]Cos  x[t]==0//Simplify;
sol=NDSolve[{eqn,x[0]==1/2,x'[0]==0},x[t],{t,0,1000}][[1]]
ParametricPlot[Evaluate[{x[t]/.sol,D[x[t]/.sol,t]}],{t,0,25},Frame->True,AxesLabel->{"x",Overscript[x,"."]},AspectRatio->1,PlotStyle->{{Red,AbsoluteThickness[2]}},TextStyle->{FontFamily->"Courier",FontWeight->"Bold",FontSize->16}]
ParametricPlot[Evaluate[{t,x[t]/.sol}],{t,0,10},Frame->True,AxesLabel->{"t","x"},AspectRatio->.5,PlotStyle->{{Green,AbsoluteThickness[3]}},TextStyle->{FontFamily->"Courier",FontWeight->"Bold",FontSize->16}]

Regards,

 rahul
  • Prev by Date: Re: Mathematica and Lisp
  • Next by Date: Re: Import of URL no longer works in V9
  • Previous by thread: Re: Alternative to Table
  • Next by thread: Re: Problem in solving Differential Equation