|
| 1 | +ctmc |
| 2 | + |
| 3 | +const int var_bound = 20; |
| 4 | + |
| 5 | +const double k1 = 0.1; |
| 6 | +const double k2 = 0.0002; |
| 7 | +const double k3 = 1; |
| 8 | +const double k4 = 0.0002; |
| 9 | +const double k5 = 1; |
| 10 | +const double k6 = 0.0002; |
| 11 | +const double k7 = 0.01; |
| 12 | +const double k8 = 0.1; |
| 13 | +const double k9 = 0.02; |
| 14 | +const double k10 = 0.1; |
| 15 | +const double k11 = 0.01; |
| 16 | +const double k12 = 0.1; |
| 17 | + |
| 18 | +module motility |
| 19 | + |
| 20 | +codY : [0..var_bound] init 1; |
| 21 | +flache : [0..var_bound] init 1; |
| 22 | +SigD_hag : [0..var_bound] init 1; |
| 23 | +hag : [0..var_bound] init 1; |
| 24 | +CodY_flache : [0..var_bound] init 1; |
| 25 | +CodY_hag : [0..var_bound] init 1; |
| 26 | +CodY : [0..var_bound] init 10; |
| 27 | +SigD : [0..var_bound] init 10; |
| 28 | +Hag : [0..var_bound] init 10; |
| 29 | + |
| 30 | +[R1] codY>0 & CodY<var_bound -> (k1*codY) : (CodY'=CodY+1); |
| 31 | +[R2] CodY>0 -> (k2*CodY) : (CodY'=CodY-1); |
| 32 | +[R3] flache>0 & SigD<var_bound -> (k3*flache) : (SigD'=SigD+1); |
| 33 | +[R4] SigD>0 -> (k4*SigD) : (SigD'=SigD-1); |
| 34 | +[R5] SigD_hag>0 & SigD<var_bound & hag<var_bound & Hag<var_bound-> (k5*SigD_hag) : (SigD_hag'=SigD_hag-1)&(SigD'=SigD+1)&(hag'=hag+1)&(Hag'=Hag+1); |
| 35 | +[R6] Hag>0 -> (k6*Hag) : (Hag'=Hag-1); |
| 36 | +[R7] SigD>0 & hag>0 & SigD_hag<var_bound-> (k7*SigD*hag) : (SigD'=SigD-1)&(hag'=hag-1)&(SigD_hag'=SigD_hag+1); |
| 37 | +[R8] SigD_hag>0 & SigD<var_bound & hag<var_bound -> (k8*SigD_hag) : (SigD_hag'=SigD_hag-1)&(SigD'=SigD+1)&(hag'=hag+1); |
| 38 | +[R9] CodY>0 & flache>0 & CodY_flache<var_bound -> (k9*CodY*flache) : (CodY'=CodY-1)&(flache'=flache-1)&(CodY_flache'=CodY_flache+1); |
| 39 | +[R10] CodY_flache>0 & CodY<var_bound & flache<var_bound -> (k10*CodY_flache) : (CodY_flache'=CodY_flache-1)& (CodY'=CodY+1) & (flache'=flache+1); |
| 40 | +[R11] CodY>0 & hag>0 & CodY_hag<var_bound -> (k11*CodY*hag) : (CodY'=CodY-1)&(hag'=hag-1)&(CodY_hag'=CodY_hag+1); |
| 41 | +[R12] CodY_hag>0 & CodY<var_bound & hag<var_bound -> (k12*CodY_hag) : (CodY_hag'=CodY_hag-1)& (CodY'=CodY+1) & (hag'=hag+1); |
| 42 | + |
| 43 | + |
| 44 | + |
| 45 | +endmodule |
| 46 | + |
| 47 | + |
| 48 | + |
| 49 | + |
| 50 | + |
| 51 | + |
0 commit comments