reactionm.h

用Monte Carlo方法

//单步时间
#define DT 1.0e-2	//min
#define DMIN 1.0

//反应速率常数
#define Ktr1 400
#define Ktr2 100
#define Ktr3 13
#define Ktr4 30
#define Ktr5 25
//每分钟

#define KTL 20
#define amRNA 0.5
#define aprot 0.03
#define ap(T) ((T>40)?1:0.5)
#define aFtsH 3.0
#define aHslVU(T) ((T>40)?2:0.7)
#define Kfold 15e3
#define K(T) ((T>40)?150:75)
#define enta(T) ((T>40)?1.75:0.35)
#define K1 1.0e9
#define K2 1.0e9
#define K3 1.0e6
#define K4 1.0e6
#define K5 1.0e6
#define K6 1.0e5
#define K7 1.0e5
#define K8 1.0e6
#define K9 1.0e9
#define K10 1.0e9
#define K11 1.0e5
#define K12 1.0e5