import nsrunit;
unit conversion on;
unit minute=60 second^1;
unit per_minute=.01666667 second^(-1);
unit per_litre=1E3 meter^(-3);
// unit micromolar predefined
unit micromolar_per_minute=1.6666667E-5 meter^(-3)*second^(-1)*mole^1;
unit minute=60 second^1;
unit per_minute=.01666667 second^(-1);
unit per_litre=1E3 meter^(-3);
// unit micromolar predefined
unit micromolar_per_minute=1.6666667E-5 meter^(-3)*second^(-1)*mole^1;
unit minute=60 second^1;
unit per_minute=.01666667 second^(-1);
unit per_litre=1E3 meter^(-3);
// unit micromolar predefined
unit micromolar_per_minute=1.6666667E-5 meter^(-3)*second^(-1)*mole^1;

math main {
	realDomain time minute;
	time.min=0;
	extern time.max;
	extern time.delta;
	real V_0 micromolar_per_minute;
	V_0=2;
	real V_1 micromolar_per_minute;
	V_1=2;
	real beta dimensionless;
	beta=0.6;
	real V_in micromolar_per_minute;
	real V_M2 micromolar_per_minute;
	V_M2=6;
	real Z(time) micromolar;
	when(time=time.min) Z=0.15;
	real K_2 micromolar;
	K_2=0.1;
	real V_2(time) micromolar_per_minute;
	real V2.zero micromolar;
	V2.zero=0;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real V2.one_dm dimensionless;
//	V2.one_dm=1;
	real V2.one_mm micromolar;
	V2.one_mm=1;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real V2.two dimensionless;
//	V2.two=2;
	real V_M3 micromolar_per_minute;
	V_M3=20;
	real K_Z micromolar;
	K_Z=0.5;
	real K_A micromolar;
	K_A=0.2;
	real K_Y micromolar;
	K_Y=0.2;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real m dimensionless;
//	m=2;
	real Y(time) micromolar;
	when(time=time.min) Y=1;
	real V3.A(time) micromolar;
	when(time=time.min) V3.A=0.42;
	real V_3(time) micromolar_per_minute;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real V3.two dimensionless;
//	V3.two=2;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real four dimensionless;
//	four=4;
	real V_M5 micromolar_per_minute;
	V_M5=5;
	real K_5 micromolar;
	K_5=1;
	real K_d micromolar;
	K_d=0.4;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real p dimensionless;
//	p=2;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real n dimensionless;
//	n=4;
	real V_5(time) micromolar_per_minute;
	real V5.zero micromolar;
	V5.zero=0;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real V5.one_dm dimensionless;
//	V5.one_dm=1;
	real V5.one_mm micromolar;
	V5.one_mm=1;
	real k per_minute;
	k=10;
	real k_f per_minute;
	k_f=1;
	real Y_init micromolar;
	Y_init=1;
	real one dimensionless;
	one=1;
	real epsilon per_minute;
	epsilon=0.1;
	real V_4 micromolar_per_minute;
	V_4=2;
	real InsP3_conc.A_init micromolar;
	InsP3_conc.A_init=0.42;

	// <component name="environment">

	// <component name="Vin">
	V_in=(V_0+V_1*beta);

	// <component name="V2">

	// <component name="V3">

	// <component name="V5">

	// <component name="cytosol">
	Z:time=(V_in-V_2+V_3+k_f*Y-k*Z);

	// <component name="internal_pool">

	// <component name="InsP3_conc">

	// <component name="V_2">
	V_2=(V_M2*Z^2/(K_2^2+Z^2)*V2.one_mm^1/(V2.zero^1+V2.one_mm^1)*V2.one_mm^1/(V2.zero^1+V2.one_mm^1));

	// <component name="V_3">
	V_3=(V_M3*Z^2/(K_Z^2+Z^2)*Y^2/(K_Y^2+Y^2)*V3.A^4/(K_A^4+V3.A^4));

	// <component name="V_5">
	V_5=(V_M5*V3.A^2/(K_5^2+V3.A^2)*Z^4/(K_d^4+Z^4)*V5.one_mm^1/(V5.zero^1+V5.one_mm^1));

	// <component name="internalpool">
	Y:time=(one*V_2-V_3-k_f*Y);

	// <component name="InsP3conc">
	V3.A:time=(beta*V_4-V_5-epsilon*V3.A);
}