import nsrunit;
unit conversion on;
unit s=1 second^1;
unit uM=1E-3 meter^(-3)*mole^1;
unit um3_uM=1E-21 mole^1;
unit um2=1E-12 meter^2;
unit um3=1E-18 meter^3;
unit per_um=1E6 meter^(-1);
unit per_um2=1E12 meter^(-2);
unit per_um2_per_s=1E12 meter^(-2)*second^(-1);
unit per_uM_per_s=1E3 meter^3*second^(-1)*mole^(-1);
unit per_s=1 second^(-1);
unit uM_per_s=1E-3 meter^(-3)*second^(-1)*mole^1;
unit per_uM_per_um3=1E21 mole^(-1);
unit uM_um2=1E-15 meter^(-1)*mole^1;
unit um2_per_s=1E-12 meter^2*second^(-1);

math main {
	realDomain t s;
	t.min=0;
	extern t.max;
	extern t.delta;
	real Ca(t) uM;
	when(t=t.min) Ca=0.1;
	real L(t) uM;
	real Ls uM;
	Ls=0.1;
	real ts s;
	ts=10;
	real Gd(t) per_um2;
	when(t=t.min) Gd=10000;
	real Gt(t) per_um2;
	when(t=t.min) Gt=0;
	real R(t) per_um2;
	when(t=t.min) R=13.9;
	real Rl(t) per_um2;
	when(t=t.min) Rl=0;
	real Rg(t) per_um2;
	when(t=t.min) Rg=5.06;
	real Rlg(t) per_um2;
	when(t=t.min) Rlg=0;
	real Rlgp(t) per_um2;
	when(t=t.min) Rlgp=0;
	real IP3(t) uM;
	when(t=t.min) IP3=0.015;
	real PIP2 per_um2;
	PIP2=4000;
	real Pc(t) per_um2;
	when(t=t.min) Pc=9.09;
	real Pcg(t) per_um2;
	when(t=t.min) Pcg=0;
	real P(t) per_um2;
	when(t=t.min) P=90.9;
	real Pg(t) per_um2;
	when(t=t.min) Pg=0;
	real J1(t) per_um2_per_s;
	real kf1 per_uM_per_s;
	kf1=0.0003;
	real kr1 per_s;
	real Kd1 uM;
	Kd1=3e-5;
	real J2(t) per_um2_per_s;
	real kf2 um2_per_s;
	kf2=2.75e-4;
	real kr2 per_s;
	real Kd2 per_um2;
	Kd2=27500;
	real J3(t) per_um2_per_s;
	real kf3 um2_per_s;
	kf3=1;
	real kr3 per_s;
	kr3=0.001;
	real J4(t) per_um2_per_s;
	real kf4 per_uM_per_s;
	kf4=0.3;
	real kr4 per_s;
	real Kd4 uM;
	Kd4=3e-5;
	real J5(t) per_um2_per_s;
	real kf5 per_s;
	kf5=0.0004;
	real J6(t) per_um2_per_s;
	real kf6 per_s;
	kf6=1;
	real J7(t) per_um2_per_s;
	real kf7 per_s;
	kf7=0.15;
	real J8(t) per_um2_per_s;
	real kf8 per_uM_per_s;
	kf8=0.0167;
	real kr8 per_s;
	kr8=0.0167;
	real J9(t) per_um2_per_s;
	real kf9 um2_per_s;
	kf9=0.0042;
	real kr9 per_s;
	kr9=1;
	real J10(t) per_um2_per_s;
	real kf10 um2_per_s;
	kf10=0.042;
	real kr10 per_s;
	kr10=1;
	real J11(t) per_um2_per_s;
	real kf11 per_uM_per_s;
	kf11=0.0334;
	real kr11 per_s;
	real Kd11 uM;
	Kd11=0.1;
	real J12(t) per_um2_per_s;
	real kf12 per_s;
	kf12=6;
	real J13(t) per_um2_per_s;
	real kf13 per_s;
	kf13=6;
	real J14(t) per_um2_per_s;
	real kf14 per_s;
	kf14=0.444;
	real Km14 uM;
	Km14=19.8;
	real J15(t) per_um2_per_s;
	real kf15 per_s;
	kf15=3.8;
	real Km15 uM;
	Km15=5;
	real J16(t) uM_per_s;
	real kf16 per_s;
	kf16=1.25;
	real Cpc uM_um2;
	real Cc uM;
	real Cp per_um2;
	real Vc um3;
	Vc=2550;
	real Rpc per_um;
	Rpc=4.61;

	// <component name="all">
	Cc=(1/(Vc*(602.2 per_uM_per_um3)));
	Cp=(1/(Vc*Rpc));
	Cpc=(Cc/Cp);
	P:t=(J13-(J9+J8));
	Pg:t=(J9-(J11+J13));
	Pc:t=(J8+J12-J10);
	Pcg:t=(J10+J11-J12);
	J13=(kf13*Pg);
	J12=(kf12*Pcg);
	kr11=(kf11*Kd11);
	J11=(kf11*Pg*Ca-kr11*Pcg);
	J10=(kf10*Pc*Gt-kr10*Pcg);
	J8=(kf8*P*Ca-kr8*Pc);
	J9=(kf9*P*Gt-kr9*Pg);
	J16=(kf16*IP3);
	J14=(kf14*Pc*PIP2/(Km14/Cpc+PIP2));
	J15=(kf15*Pcg*PIP2/(Km15/Cpc+PIP2));
	IP3:t=(Cpc*(J14+J15)-J16);
	Gd:t=(J7+J13+J12-(J2+J3));
	Gt:t=(J6-(J7+J9+J10));
	J7=(kf7*Gt);
	Ca:t=(Cpc*((-1)*1)*(J8+J11));
	L=(if ((t<(ts+(.15 s))) and (t>=ts)) Ls/(1+exp((-1)*(80 per_s)*(t-ts-(.05 s)))) else if (t>=(ts+(.15 s))) Ls else (0 uM));
	kr1=(kf1*Kd1);
	J1=(kf1*R*L-kr1*Rl);
	kr2=(kf2*Kd2);
	J2=(kf2*R*Gd-kr2*Rg);
	R:t=((-1)*1*(J1+J2));
	J3=(kf3*Rl*Gd-kr3*Rlg);
	Rl:t=(J1+J6-J3);
	kr4=(kf4*Kd4);
	J4=(kf4*L*Rg-kr4*Rlg);
	Rg:t=(J2-J4);
	J5=(kf5*Rlg);
	Rlgp:t=J5;
	J6=(kf6*Rlg);
	Rlg:t=(J3-J5+J4-J6);
}