import nsrunit;
// Warning: unit conversion turned off due to unit errors in 1 equation(s)
unit conversion off;
unit minute=60 second^1;
unit first_order_rate_constant=.01666667 second^(-1);
unit flux=1.6666667E-8 second^(-1)*mole^1;
// unit micromolar predefined

math main {
	//Warning:  the following variables were set 'extern' or given
	//  an initial value of '0' because the model would otherwise be
	//  underdetermined:  X, M
	realDomain time minute;
	time.min=0;
	extern time.max;
	extern time.delta;
	real C.C(time) micromolar;
	when(time=time.min) C.C=0.01;
	real X(time) dimensionless;
	//Warning:  Assuming zero initial condition; nothing provided in original CellML model.
	when(time=time.min) X=0;
	real kd first_order_rate_constant;
	kd=0.01;
	real Kd micromolar;
	Kd=0.02;
	real vi flux;
	vi=0.025;
	real vd flux;
	vd=0.25;
	real M(time) dimensionless;
	//Warning:  Assuming zero initial condition; nothing provided in original CellML model.
	when(time=time.min) M=0;
	real M_star(time) dimensionless;
	real V1(time) first_order_rate_constant;
	real V2 first_order_rate_constant;
	V2=1.5;
	real K1 dimensionless;
	K1=0.005;
	real K2 dimensionless;
	K2=0.005;
	real M_init dimensionless;
	M_init=0.01;
	real X_init dimensionless;
	X_init=0.01;
	real X_star(time) dimensionless;
	real V3(time) first_order_rate_constant;
	real V4 first_order_rate_constant;
	V4=0.5;
	real K3 dimensionless;
	K3=0.005;
	real K4 dimensionless;
	K4=0.005;
	real Kc micromolar;
	Kc=0.5;
	real VM1 first_order_rate_constant;
	VM1=3;
	real VM3 first_order_rate_constant;
	VM3=1;

	// <component name="environment">

	// <component name="C">
	C.C:time=(vi-vd*X*C.C/(Kd+C.C)-kd*C.C);

	// <component name="M">

	// <component name="M_star">

	// <component name="X">

	// <component name="X_star">

	// <component name="parameters">
	V1=(VM1*C.C/(Kc+C.C));
	V3=(M*VM3);

	// <component name="law1_M">
	M:time=(V1*M_star/(K1+M_star)-V2*M/(K2+M));

	// <component name="law2_M">
	M_star=(1-M);

	// <component name="law1_X">
	X:time=(V3*X_star/(K3+X_star)-V4*X/(K4+X));

	// <component name="law2_Xstar">
	X_star=(1-X);
}