MathGroup Archive 2010

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

Search the Archive

Re: sketch a phase portrait for the nonlinear system

Hi Hongtao,

this is more a math question than a Mathematica question ;)

The system has 3 stationary points: saddle at (0,0) and 2 spirals at
(-1,0) and (1,0). The initial conditions (+-0.1, +- 0.16) are chosen
such that they lie very close to the stable and unstable manifolds of
the saddle point at (0,0). That's why you see the nice cross at

I don't know how much you tried to change them, but as long as you
keep them inside the plot range -3<=x0<=3 and -3<=y0<=3, you should
see the trajectory in the plot.


On Jan 15, 12:17=C2 am, =E6=B6=9B =E6=B4=AA <hongtao... at> wro=
> Hi,everyone,
> I am new in Mathematica.When I try to sketch a phase portrait for the non=
linear system, I don't know how to plan the input values for the ODEs to ge=
t the right FIGURE.
> My example is from a book as follows:
> ode1[x0_,y0_]:=NDSolve[{x=E2=80=99[t]==y[t],y=E2=80=99[=
> sol[1]=ode1[0.1,0.16];sol[2]=ode1[0.1,-0.06];
> sol[3]=ode1[-0.1,-0.16];sol[4]=ode1[-0.1,0.06];
> sol[5]=ode1[1.3,0];sol[6]=ode1[-1.3,0];
> p=ParametricPlot[Evaluate[Table[{x[t],y[t]}/.sol[i],{i,6}]],{t,-10,10},=
> Show[p,PlotRange->{{-3,3},{-3,3}},AxesLabel->{"x","y"}]
> Would you like to tell me why we choose the input values for the ODEs suc=
h as ode1[0.1,0.16]; ode1[0.1,-0.06];...? You see, when I tried to change t=
he values ,I can't see the figure any more.
> Best Regards,
> hongtao

  • Prev by Date: Re: More /.{I->-1} craziness
  • Next by Date: Re: Re: restricting interpolating functions to be positive
  • Previous by thread: sketch a phase portrait for the nonlinear system
  • Next by thread: Converting date strings to DateList format - a need for speed