This is my current homework problem. I have completed parts a and b, but am now stuck on part c. I'm not sure what it means by plot the trajectories with different initial conditions, nor am I 100% sure how to do that. I know because of my positive eigenvalues that the solution will be unstable but after that I'm lost. Any advice would be really great.
$\begingroup$
$\endgroup$
2
-
2$\begingroup$ Welcome to MSE. Take a look at this. $\endgroup$– Rohit NamjoshiFeb 11, 2020 at 3:32
-
$\begingroup$ This problem is quite similar to Lorenz system. Please take a look at the "Mathematica simulation" part. $\endgroup$– Αλέξανδρος ΖεγγFeb 11, 2020 at 10:27
Add a comment
|
1 Answer
$\begingroup$
$\endgroup$
You could try NDSolve
then ParametricPlot3D
Manipulate[
Module[{a = 0.25, b = 0.5, ode1, ode2, ode3, x, y, z, t},
ode1 = x'[t] == -y[t] - z[t];
ode2 = y'[t] == x[t] + a y[t];
ode3 = z'[t] == b + z[t] (x[t] - c);
sol = NDSolve[{ode1, ode2, ode3, x[0] == x0, y[0] == y0,
z[0] == z0}, {x[t], y[t], z[t]}, {t, 0, maxTime}, MaxSteps -> Infinity];
ParametricPlot3D[{x[t], y[t], z[t]} /. sol, {t, 0, maxTime},
PlotRange -> {{-15, 15}, {-15, 15}, {-10, 30}},
AxesLabel -> {"x(t)", "y(t)", "z(t)"}, BaseStyle -> 12,
ImageSize -> 400, PerformanceGoal -> "Quality", Mesh -> None]
],
{{x0, .1, "x(0)"}, 0, 2, .1, Appearance -> "Labeled", ImageSize -> Tiny},
{{y0, .1, "y(0)"}, 0, 1, .1, Appearance -> "Labeled", ImageSize -> Tiny},
{{z0, .5, "z(0)"}, 0, 2, .1, Appearance -> "Labeled", ImageSize -> Tiny},
{{c, 6.5, "c"}, 0, 7, .1, Appearance -> "Labeled", ImageSize -> Tiny},
{{maxTime, 200, "time"}, 0.1, 250, .1, Appearance -> "Labeled", ImageSize -> Tiny},
TrackedSymbols :> {x0, y0, z0, c, maxTime},
ContinuousAction -> True,
Alignment -> Center,
SynchronousUpdating -> True,
SynchronousInitialization -> True,
FrameMargins -> 1, ImageMargins -> 1,
ControlPlacement -> Left
]