1
1
Getting Started
2
2
================================================================================
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).
4
6
5
7
Mathematical formulation
6
- -----------------------------
8
+ ------------------------
7
9
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
10
12
11
13
.. math ::
12
14
@@ -95,9 +97,10 @@ diagram by visualizing all spatial-components of the state vector.
95
97
problem.plotPhaseSpace(sol);
96
98
97
99
.. 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
101
104
102
105
``OTP `` also supports animations for the computed solution:
103
106
@@ -106,9 +109,10 @@ diagram by visualizing all spatial-components of the state vector.
106
109
% Create a movie of the solution
107
110
problem.movie(sol);
108
111
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
112
116
113
117
Changing the parameters
114
118
------------------------
@@ -125,9 +129,12 @@ For example, changing the parameter :math:`\rho` in the Lorenz system leads to a
125
129
sol = problem.solve('MaxStep' , 1e-4);
126
130
problem.movie(sol);
127
131
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
+
131
138
Changing the solver
132
139
-------------------
133
140
@@ -141,10 +148,11 @@ following code:
141
148
142
149
sol = ode23s(problem.RHS.F, problem.TimeSpan, problem.Y0, ...
143
150
odeset('Jacobian', problem.RHS.Jacobian));
151
+
144
152
This is particularly useful when you want to compare the performance of different solvers on the same problem.
145
153
146
154
Next steps
147
155
------------------------
148
156
Explore different problems available in OTP by browsing the Problems Gallery in the sidebar. You can define your custom
149
157
problems 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