Skip to content

Commit c77a9ee

Browse files
add some examples
1 parent aa793d8 commit c77a9ee

File tree

1 file changed

+18
-5
lines changed

1 file changed

+18
-5
lines changed

README.md

Lines changed: 18 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -20,14 +20,27 @@ the documentation, which contains the unreleased features.
2020
The following are valid ODE definitions.
2121

2222
```julia
23-
f = @ode_def begin
23+
using DifferentialEquations
24+
25+
# Non-Stiff ODE
26+
27+
lotka_volterra = @ode_def begin
28+
d🐁 = α*🐁 - β*🐁*🐈
29+
d🐈 = -γ*🐈 + δ*🐁*🐈
30+
end α β γ δ
31+
32+
p = [1.5,1.0,3.0,1.0]; u0 = [1.0;1.0]
33+
prob = ODEProblem(lotka_volterra,u0,(0.0,10.0),p)
34+
sol = solve(prob,Tsit5(),reltol=1e-6,abstol=1e-6)
35+
36+
# Stiff ODE
37+
38+
rober = @ode_def begin
2439
dy₁ = -k₁*y₁+k₃*y₂*y₃
2540
dy₂ = k₁*y₁-k₂*y₂^2-k₃*y₂*y₃
2641
dy₃ = k₂*y₂^2
2742
end k₁ k₂ k₃
2843

29-
f = @ode_def begin
30-
d🐁 = α*🐁 - β*🐁*🐈
31-
d🐈 = -γ*🐈 + δ*🐁*🐈
32-
end α β γ δ
44+
prob = ODEProblem(rober,[1.0,0.0,0.0],(0.0,1e5),[0.04,3e7,1e4])
45+
sol = solve(prob)
3346
```

0 commit comments

Comments
 (0)