$fast=10; $slow=1; $CycD_del=0; $Cdc20_del=0; $Rb_del=0; // Initial conditions // A.istate = 1; // if nothing, initial condition is randomly chosen CycD.istate=1; CycE.istate=0; CycA.istate=0; CycB.istate = 0; Rb.istate=1; E2F.istate=0; p27.istate=1; Cdc20.istate=0; UbcH10.istate=0; cdh1.istate=1; Rb.is_internal=TRUE; E2F.is_internal = TRUE; p27.is_internal=TRUE; Cdc20.is_internal = TRUE; UbcH10.is_internal = TRUE; cdh1.is_internal=TRUE; // number of trajectories sample_count = 50000; // maximum time of a trajectory - // after max_time is reached, the computation stops max_time = 1000; // time of an interval on which probabilities are computed and averaged time_tick = .2; // if discrete time is 1, the computation is fully discrete // in this case, the time interval is equal to time tick // if discrete time is 0, the computation is continuous // in this case, the time interval is chosen randomly discrete_time = 0; use_physrandgen = TRUE; seed_pseudorandom = 100; statdist_traj_count = 1000; statdist_cluster_threshold = 0.98; thread_count=4;