TD1_lotka_volterra.out.ref 922 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
license GNU GPL 2, http://lifeware.inria.fr/biocham4/
k1*A*B for A+B=>2*B.
k2*A for A=>2*A.
k3*B for B=>_.
present(A,a).
present(B,b).
parameter(
  a = 1,
  b = 1,
  k1 = 2,
  k2 = 2,
  k3 = 1
).
<img src="graph0.png" alt="Simulation results">
<img src="graph1.png" alt="Simulation results">
<img src="ode.tex" alt="Simulation results">
No complex invariant found
<img src="plot-0.csv" alt="Simulation results">
<img src="plot-1.csv" alt="Simulation results">
<img src="plot-2.csv" alt="Simulation results">
<img src="plot-3.csv" alt="Simulation results">
22
23
present({}).
present({A}).
24
25
26
27
28
29
30
31
32
33
reachable(stable(A))
reachable(stable(not A))
reachable(stable(not B))
reachable(steady(B))
checkpoint2(B,not A)
checkpoint2(A,not B)
AG(f)
EG(f)
EF(EG(s))
not EU(not x,y)
SOLIMAN Sylvain's avatar
SOLIMAN Sylvain committed
34
not y/\EF(y)/\not EU(not x,y)
35
36
37
reachable(checkpoint2(not A,not B)) is true
reachable(checkpoint2(not B,not A)) is false
<img src="plot-4.csv" alt="Simulation results">