11# # Lotka-Volterra
22
3- lotka = @ode_def_nohes LotkaVolterra begin
3+ lotka = @ode_def_all LotkaVolterra begin
44 dx = a* x - b* x* y
55 dy = - c* y + d* x* y
66end 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])
1919
2020# # Fitzhugh-Nagumo
2121
22- fitz = @ode_def_nohes FitzhughNagumo begin
22+ fitz = @ode_def_all FitzhughNagumo begin
2323 dv = v - v^ 3 / 3 - w + l
2424 dw = τinv* (v + a - b* w)
2525end a b τinv l
@@ -36,7 +36,7 @@ with initial condition ``v=w=1``
3636prob_ode_fitzhughnagumo = ODEProblem (fitz,[1.0 ;1.0 ],(0.0 ,1.0 ),(0.7 ,0.8 ,1 / 12.5 ,0.5 ))
3737
3838# Van der Pol Equations
39- van = @ode_def_noinvhes VanDerPol begin
39+ van = @ode_def_all VanDerPol begin
4040 dy = μ* ((1 - x^ 2 )* y - x)
4141 dx = 1 * y
4242end μ
@@ -74,7 +74,7 @@ prob_ode_vanstiff = ODEProblem(van,[0;sqrt(3)],(0.0,1.0),1e6)
7474
7575# ROBER
7676
77- rober = @ode_def_noinvjac Rober begin
77+ rober = @ode_def_all Rober begin
7878 dy₁ = - k₁* y₁+ k₃* y₂* y₃
7979 dy₂ = k₁* y₁- k₂* y₂^ 2 - k₃* y₂* y₃
8080 dy₃ = k₂* y₂^ 2
@@ -137,7 +137,7 @@ prob_ode_threebody = ODEProblem(threebody,[0.994, 0.0, 0.0, big(-2.0015851063790
137137
138138# Rigid Body Equations
139139
140- rigid = @ode_def_noinvjac RigidBody begin
140+ rigid = @ode_def_all RigidBody begin
141141 dy₁ = I₁* y₂* y₃
142142 dy₂ = I₂* y₁* y₃
143143 dy₃ = I₃* y₁* y₂
@@ -276,7 +276,7 @@ f(y) = \\begin{pmatrix}
276276
277277http://www.radford.edu/~thompson/vodef90web/problems/demosnodislin/Demos_Pitagora/DemoHires/demohires.pdf
278278"""
279- hires = @ode_def_noinvjac Hires begin
279+ hires = @ode_def_all Hires begin
280280 dy1 = - p1* y1 + p2* y2 + p3* y3 + p4
281281 dy2 = p1* y1 - p5* y2
282282 dy3 = - p6* y3 + p2* y4 + p7* y5
@@ -316,7 +316,7 @@ where ``s=77.27``, ``w=0.161`` and ``q=8.375⋅10^{-6}``.
316316
317317http://www.radford.edu/~thompson/vodef90web/problems/demosnodislin/Demos_Pitagora/DemoOrego/demoorego.pdf
318318"""
319- orego = @ode_def Orego begin
319+ orego = @ode_def_all Orego begin
320320 dy1 = p1* (y2+ y1* (1 - p2* y1- y2))
321321 dy2 = (y3- (1 + y1)* y2)/ p1
322322 dy3 = p3* (y1- y3)
0 commit comments