1
1
# # Lotka-Volterra
2
2
3
- lotka = @ode_def_nohes LotkaVolterra begin
3
+ lotka = @ode_def_all LotkaVolterra begin
4
4
dx = a* x - b* x* y
5
5
dy = - c* y + d* x* y
6
6
end a b c d
@@ -19,7 +19,7 @@ prob_ode_lotkavoltera = ODEProblem(lotka,[1.0,1.0],(0.0,1.0),[1.5,1.0,3.0,1.0])
19
19
20
20
# # Fitzhugh-Nagumo
21
21
22
- fitz = @ode_def_nohes FitzhughNagumo begin
22
+ fitz = @ode_def_all FitzhughNagumo begin
23
23
dv = v - v^ 3 / 3 - w + l
24
24
dw = τinv* (v + a - b* w)
25
25
end a b τinv l
@@ -36,7 +36,7 @@ with initial condition ``v=w=1``
36
36
prob_ode_fitzhughnagumo = ODEProblem (fitz,[1.0 ;1.0 ],(0.0 ,1.0 ),(0.7 ,0.8 ,1 / 12.5 ,0.5 ))
37
37
38
38
# Van der Pol Equations
39
- van = @ode_def_noinvhes VanDerPol begin
39
+ van = @ode_def_all VanDerPol begin
40
40
dy = μ* ((1 - x^ 2 )* y - x)
41
41
dx = 1 * y
42
42
end μ
@@ -74,7 +74,7 @@ prob_ode_vanstiff = ODEProblem(van,[0;sqrt(3)],(0.0,1.0),1e6)
74
74
75
75
# ROBER
76
76
77
- rober = @ode_def_noinvjac Rober begin
77
+ rober = @ode_def_all Rober begin
78
78
dy₁ = - k₁* y₁+ k₃* y₂* y₃
79
79
dy₂ = k₁* y₁- k₂* y₂^ 2 - k₃* y₂* y₃
80
80
dy₃ = k₂* y₂^ 2
@@ -137,7 +137,7 @@ prob_ode_threebody = ODEProblem(threebody,[0.994, 0.0, 0.0, big(-2.0015851063790
137
137
138
138
# Rigid Body Equations
139
139
140
- rigid = @ode_def_noinvjac RigidBody begin
140
+ rigid = @ode_def_all RigidBody begin
141
141
dy₁ = I₁* y₂* y₃
142
142
dy₂ = I₂* y₁* y₃
143
143
dy₃ = I₃* y₁* y₂
@@ -276,7 +276,7 @@ f(y) = \\begin{pmatrix}
276
276
277
277
http://www.radford.edu/~thompson/vodef90web/problems/demosnodislin/Demos_Pitagora/DemoHires/demohires.pdf
278
278
"""
279
- hires = @ode_def_noinvjac Hires begin
279
+ hires = @ode_def_all Hires begin
280
280
dy1 = - p1* y1 + p2* y2 + p3* y3 + p4
281
281
dy2 = p1* y1 - p5* y2
282
282
dy3 = - p6* y3 + p2* y4 + p7* y5
@@ -316,7 +316,7 @@ where ``s=77.27``, ``w=0.161`` and ``q=8.375⋅10^{-6}``.
316
316
317
317
http://www.radford.edu/~thompson/vodef90web/problems/demosnodislin/Demos_Pitagora/DemoOrego/demoorego.pdf
318
318
"""
319
- orego = @ode_def Orego begin
319
+ orego = @ode_def_all Orego begin
320
320
dy1 = p1* (y2+ y1* (1 - p2* y1- y2))
321
321
dy2 = (y3- (1 + y1)* y2)/ p1
322
322
dy3 = p3* (y1- y3)
0 commit comments