Skip to content

Commit e52a177

Browse files
Update dae_premade_problems.jl
1 parent b33d1d8 commit e52a177

File tree

1 file changed

+17
-1
lines changed

1 file changed

+17
-1
lines changed

src/dae_premade_problems.jl

Lines changed: 17 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -13,5 +13,21 @@ f = function (r, yp, y, p, tres)
1313
end
1414
u0 = [1.0, 0, 0]
1515
du0 = [-0.04, 0.04, 0.0]
16-
"DAE residual form for the Robertson model"
16+
17+
@doc doc"""
18+
The Robertson biochemical reactions in DAE form
19+
20+
```math
21+
\\frac{dy₁}{dt} = -k₁y₁+k₃y₂y₃
22+
```
23+
```math
24+
\\frac{dy₂}{dt} = k₁y₁-k₂y₂^2-k₃y₂y₃
25+
```
26+
```math
27+
1 = y₁ + y₂ + y₃
28+
```
29+
where ``k₁=0.04``, ``k₂=3\\times10^7``, ``k₃=10^4``. For details, see:
30+
Hairer Norsett Wanner Solving Ordinary Differential Equations I - Nonstiff Problems Page 129
31+
Usually solved on ``[0,1e11]``
32+
"""
1733
prob_dae_resrob = DAEProblem(f,du0,u0,(0.0,100000.0))

0 commit comments

Comments
 (0)