11Getting Started
22================================================================================
3- This guide shows you how to create, solve and visualize different problems in ODE Test Problems (``OTP ``).
3+
4+ This guide shows you how to create, solve and visualize different problems in
5+ ODE Test Problems (OTP).
46
57Mathematical formulation
6- -----------------------------
8+ ------------------------
79
8- All test problems in `` OTP `` are considered as a first-order
9- differential-algebraic equation of the form
10+ All test problems in OTP are considered as a first-order differential-algebraic
11+ equation of the form
1012
1113.. math ::
1214
@@ -95,9 +97,10 @@ diagram by visualizing all spatial-components of the state vector.
9597 problem.plotPhaseSpace(sol);
9698
9799 .. image :: ../images/Lorenz-solution-plot.png
98- :width: 200
99- .. image :: ../images/Lorenz-Phase-plot.png
100- :width: 200
100+ :width: 256
101+
102+ .. image :: ../images/Lorenz-phase-plot.png
103+ :width: 256
101104
102105``OTP `` also supports animations for the computed solution:
103106
@@ -106,9 +109,10 @@ diagram by visualizing all spatial-components of the state vector.
106109 % Create a movie of the solution
107110 problem.movie(sol);
108111
109- .. video :: ../_static/Lorenz-Original-Canonical.webm
110- :loop:
111- :width: 200
112+ .. video :: ../images/animations/Lorenz-original-canonical.webm
113+ :autoplay:
114+ :loop:
115+ :width: 256
112116
113117Changing the parameters
114118------------------------
@@ -125,9 +129,12 @@ For example, changing the parameter :math:`\rho` in the Lorenz system leads to a
125129 sol = problem.solve('MaxStep' , 1e-4);
126130 problem.movie(sol);
127131
128- .. video :: ../_static/Lorenz-Alternate-Canonical.webm
129- :loop:
130- :width: 200
132+ .. video :: ../images/animations/Lorenz-alternate-canonical.webm
133+ :autoplay:
134+ :loop:
135+ :width: 256
136+ :class: center
137+
131138Changing the solver
132139-------------------
133140
@@ -141,10 +148,11 @@ following code:
141148
142149 sol = ode23s(problem.RHS.F, problem.TimeSpan, problem.Y0, ...
143150 odeset('Jacobian', problem.RHS.Jacobian));
151+
144152 This is particularly useful when you want to compare the performance of different solvers on the same problem.
145153
146154Next steps
147155------------------------
148156Explore different problems available in OTP by browsing the Problems Gallery in the sidebar. You can define your custom
149157problems by creating a new class that inherits from the ``otp.Problem `` class.
150- See the `Contributing Guide <../contributing >`_ for more details.
158+ See the `Contributing Guide <../contributing >`_ for more details.
0 commit comments