@@ -135,14 +135,15 @@ threebody = (du,u,p,t) -> begin
135
135
du[3 ] = u[1 ] + 2 u[4 ] - threebody_μ′* (u[1 ]+ threebody_μ)/ D₁ - threebody_μ* (u[1 ]- threebody_μ′)/ D₂
136
136
du[4 ] = u[2 ] - 2 u[3 ] - threebody_μ′* u[2 ]/ D₁ - threebody_μ* u[2 ]/ D₂
137
137
end
138
- """
138
+
139
+ @doc doc"""
139
140
The ThreeBody problem as written by Hairer: (Non-stiff)
140
141
141
142
```math
142
- \f rac{dy₁}{dt} = y₁ + 2\f rac{dy₂}{dt} - \b ar{μ}\\ frac{y₁+μ}{D₁} - μ\ \ frac{y₁-\b ar{μ}}{D₂}
143
+ \f rac{dy₁}{dt} = y₁ + 2\f rac{dy₂}{dt} - \b ar{μ}\f rac{y₁+μ}{D₁} - μ\f rac{y₁-\b ar{μ}}{D₂}
143
144
```
144
145
```math
145
- \f rac{dy₂}{dt} = y₂ - 2\f rac{dy₁}{dt} - \b ar{μ}\\ frac{y₂}{D₁} - μ\ \ frac{y₂}{D₂}
146
+ \f rac{dy₂}{dt} = y₂ - 2\f rac{dy₁}{dt} - \b ar{μ}\f rac{y₂}{D₁} - μ\f rac{y₂}{D₂}
146
147
```
147
148
```math
148
149
D₁ = ((y₁+μ)^2 + y₂^2)^{3/2}
@@ -220,14 +221,15 @@ pleiades = (du,u,p,t) -> begin
220
221
end
221
222
end
222
223
end
223
- """
224
+
225
+ @doc doc"""
224
226
Pleiades Problem (Non-stiff)
225
227
226
228
```math
227
- \f rac{d^2xᵢ}{dt^2} = \\ sum_{j≠i} mⱼ(xⱼ-xᵢ)/rᵢⱼ
229
+ \f rac{d^2xᵢ}{dt^2} = \s um_{j≠i} mⱼ(xⱼ-xᵢ)/rᵢⱼ
228
230
```
229
231
```math
230
- \f rac{d^2yᵢ}{dt^2} = \\ sum_{j≠i} mⱼ(yⱼ-yᵢ)/rᵢⱼ
232
+ \f rac{d^2yᵢ}{dt^2} = \s um_{j≠i} mⱼ(yⱼ-yᵢ)/rᵢⱼ
231
233
```
232
234
233
235
where
@@ -281,19 +283,19 @@ y₆(0) = -4
281
283
y₇(0) = 4
282
284
```
283
285
284
- and with ``\f rac{dxᵢ(0)}{dt}=\f rac{dyᵢ(0)}{dt}=0`` except for
286
+ and with ``\\ frac{dxᵢ(0)}{dt}=\f rac{dyᵢ(0)}{dt}=0`` except for
285
287
286
288
```math
287
- \f rac{dx₆(0)}{dt} = 1.75
289
+ \\ frac{dx₆(0)}{dt} = 1.75
288
290
```
289
291
```math
290
- \f rac{dx₇(0)}{dt} = -1.5
292
+ \\ frac{dx₇(0)}{dt} = -1.5
291
293
```
292
294
```math
293
- \f rac{dy₄(0)}{dt} = -1.25
295
+ \\ frac{dy₄(0)}{dt} = -1.25
294
296
```
295
297
```math
296
- \f rac{dy₅(0)}{dt} = 1
298
+ \\ frac{dy₅(0)}{dt} = 1
297
299
```
298
300
299
301
From Hairer Norsett Wanner Solving Ordinary Differential Equations I - Nonstiff Problems Page 244
0 commit comments