@@ -70,7 +70,8 @@ with ``μ=1.0`` and ``u_0=[x => \\sqrt{3}, y => 0]``
70
70
71
71
Non-stiff parameters.
72
72
"""
73
- prob_ode_vanderpol = ODEProblem (van, [y => 0 , x => sqrt (3 ), μ => 1.0 ], (0.0 , 1.0 ), jac= true , eval_module = @__MODULE__ )
73
+ prob_ode_vanderpol = ODEProblem (van, [y => 0 , x => sqrt (3 ), μ => 1.0 ], (0.0 , 1.0 ),
74
+ jac = true , eval_module = @__MODULE__ )
74
75
75
76
"""
76
77
Van der Pol Equations
@@ -86,7 +87,8 @@ with ``μ=10^6`` and ``u_0=[x => \\sqrt{3}, y => 0]``
86
87
87
88
Stiff parameters.
88
89
"""
89
- prob_ode_vanderpol_stiff = ODEProblem (van, [y => 0 , x => sqrt (3 ), μ => 1e6 ], (0.0 , 1.0 ), jac= true , eval_module = @__MODULE__ )
90
+ prob_ode_vanderpol_stiff = ODEProblem (van, [y => 0 , x => sqrt (3 ), μ => 1e6 ], (0.0 , 1.0 ),
91
+ jac = true , eval_module = @__MODULE__ )
90
92
91
93
# ROBER
92
94
@parameters k₁ k₂ k₃
@@ -116,7 +118,9 @@ Hairer Norsett Wanner Solving Ordinary Differential Equations I - Nonstiff Probl
116
118
117
119
Usually solved on ``[0,1e11]``
118
120
"""
119
- prob_ode_rober = ODEProblem (rober, [[y₁, y₂, y₃] .=> [1.0 ; 0.0 ; 0.0 ]; [k₁, k₂, k₃] .=> (0.04 , 3e7 , 1e4 )], (0.0 , 1e11 ), jac = true , eval_module = @__MODULE__ )
121
+ prob_ode_rober = ODEProblem (
122
+ rober, [[y₁, y₂, y₃] .=> [1.0 ; 0.0 ; 0.0 ]; [k₁, k₂, k₃] .=> (0.04 , 3e7 , 1e4 )],
123
+ (0.0 , 1e11 ), jac = true , eval_module = @__MODULE__ )
120
124
121
125
# Three Body
122
126
const threebody_μ = big (0.012277471 );
@@ -200,8 +204,10 @@ or Hairer Norsett Wanner Solving Ordinary Differential Equations I - Nonstiff Pr
200
204
201
205
Usually solved from 0 to 20.
202
206
"""
203
- prob_ode_rigidbody = ODEProblem (rigid, [[y₁, y₂, y₃] .=> [1.0 , 0.0 , 0.9 ]; [I₁, I₂, I₃] .=> (- 2.0 , 1.25 , - 0.5 )], (0.0 , 20.0 ),
204
- jac = true , eval_module = @__MODULE__ )
207
+ prob_ode_rigidbody = ODEProblem (
208
+ rigid, [[y₁, y₂, y₃] .=> [1.0 , 0.0 , 0.9 ]; [I₁, I₂, I₃] .=> (- 2.0 , 1.25 , - 0.5 )],
209
+ (0.0 , 20.0 ),
210
+ jac = true , eval_module = @__MODULE__ )
205
211
206
212
# Pleiades Problem
207
213
@@ -365,8 +371,9 @@ u0[1] = 1
365
371
u0[8 ] = 0.0057
366
372
367
373
u0 = [y1, y2, y3, y4, y5, y6, y7, y8] .=> u0
368
- p = [p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12] .=> (1.71 , 0.43 , 8.32 , 0.0007 , 8.75 ,
369
- 10.03 , 0.035 , 1.12 , 1.745 , 280.0 , 0.69 , 1.81 )
374
+ p = [p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12] .=>
375
+ (1.71 , 0.43 , 8.32 , 0.0007 , 8.75 ,
376
+ 10.03 , 0.035 , 1.12 , 1.745 , 280.0 , 0.69 , 1.81 )
370
377
371
378
"""
372
379
Hires Problem (Stiff)
@@ -390,7 +397,8 @@ where ``f`` is defined by
390
397
Reference: [demohires.pdf](http://www.radford.edu/~thompson/vodef90web/problems/demosnodislin/Demos_Pitagora/DemoHires/demohires.pdf)
391
398
Notebook: [Hires.ipynb](http://nbviewer.jupyter.org/github/JuliaDiffEq/DiffEqBenchmarks.jl/blob/master/StiffODE/Hires.ipynb)
392
399
"""
393
- prob_ode_hires = ODEProblem (hires, [u0; p], (0.0 , 321.8122 ), jac = true , eval_module = @__MODULE__ )
400
+ prob_ode_hires = ODEProblem (
401
+ hires, [u0; p], (0.0 , 321.8122 ), jac = true , eval_module = @__MODULE__ )
394
402
395
403
@parameters p1 p2 p3
396
404
@variables y1 (t) y2 (t) y3 (t)
@@ -418,4 +426,6 @@ where ``s=77.27``, ``w=0.161`` and ``q=8.375⋅10^{-6}``.
418
426
Reference: [demoorego.pdf](http://www.radford.edu/~thompson/vodef90web/problems/demosnodislin/Demos_Pitagora/DemoOrego/demoorego.pdf)
419
427
Notebook: [Orego.ipynb](http://nbviewer.jupyter.org/github/JuliaDiffEq/DiffEqBenchmarks.jl/blob/master/StiffODE/Orego.ipynb)
420
428
"""
421
- prob_ode_orego = ODEProblem (orego, [[y1, y2, y3] .=> [1.0 , 2.0 , 3.0 ]; [p1, p2, p3] .=> [77.27 , 8.375e-6 , 0.161 ]], (0.0 , 30.0 ), jac = true , eval_module = @__MODULE__ )
429
+ prob_ode_orego = ODEProblem (
430
+ orego, [[y1, y2, y3] .=> [1.0 , 2.0 , 3.0 ]; [p1, p2, p3] .=> [77.27 , 8.375e-6 , 0.161 ]],
431
+ (0.0 , 30.0 ), jac = true , eval_module = @__MODULE__ )
0 commit comments