import nsrunit; unit conversion on; unit s=1 second^1; unit nM=1E-6 meter^(-3)*mole^1; math main { //Warning: the following variables were set 'extern' or given // an initial value of '0' because the model would otherwise be // underdetermined: t, Ca extern real t s; extern real Ca nM; real M nM; M=6000; real Ntot nM; Ntot=1000; real K_mN nM; K_mN=535; // Var below replaced by constant in model eqns to satisfy unit correction // real n dimensionless; // n=2.92; real K_dN nM; K_dN=1760; real act_N dimensionless; real act_Napp dimensionless; real tau_actN s; tau_actN=0.070059; real tNinhib s; tNinhib=3600; // act_N=(Ca^2.92/(Ca^2.92+K_mN^2.92*(1+K_dN/M))); act_Napp=(if (t=tNinhib) act_N*2.718281828459045^((-1)*(t-tNinhib)/tau_actN) else 0); }