import nsrunit;
unit conversion on;
unit month=2592000 second^1;
unit first_order_rate_constant=3.8580247E-7 second^(-1);
unit month=2592000 second^1;
unit first_order_rate_constant=3.8580247E-7 second^(-1);
unit fraction=1  dimensionless;
unit month=2592000 second^1;
unit flux=3.8580247E-7 second^(-1);
unit first_order_rate_constant=3.8580247E-7 second^(-1);
unit flux=3.8580247E-7 second^(-1);
unit flux=3.8580247E-7 second^(-1);

math main {
	realDomain time month;
	time.min=0;
	extern time.max;
	extern time.delta;
	real L(time) fraction;
	when(time=time.min) L=0.33;
	real L_init fraction;
	L_init=0.33;
	real mu_T month;
	mu_T=2.4;
	real mu_L month;
	mu_L=2;
	real mu_A month;
	mu_A=3.6;
	real L.epsilon first_order_rate_constant;
	L.epsilon=0;
	real L.M(time) fraction;
	when(time=time.min) L.M=0;
	real T(time) fraction;
	when(time=time.min) T=0.33;
	real L_0 dimensionless;
	real A(time) fraction;
	when(time=time.min) A=0.34;
	real A_init fraction;
	A_init=0.34;
	real A_0 dimensionless;
	real A.eps_zero first_order_rate_constant;
	A.eps_zero=0;
	real A.M_zero fraction;
	A.M_zero=0;
	real T_init fraction;
	T_init=0.33;
	real T_0 dimensionless;
	real T.eps_zero first_order_rate_constant;
	T.eps_zero=0;
	real T.M_zero fraction;
	T.M_zero=0;

	// <component name="environment">

	// <component name="L">

	// <component name="A">

	// <component name="T">

	// <component name="M">
	L.M:time=(L.epsilon*T);

	// <component name="reaction_constants">

	// <component name="Ai">
	A:time=(1/mu_L*L-1/mu_A*A+A.eps_zero*A.M_zero);

	// <component name="Ti">
	T:time=(1/mu_A*A-1/mu_T*T+T.eps_zero*T.M_zero);

	// <component name="Li">
	L:time=(1/mu_T*T-1/mu_L*L+L.epsilon*L.M);

	// <component name="A_0i">
	A_0=(mu_L/(mu_A+mu_T+mu_L));

	// <component name="T_0i">
	T_0=(mu_L/(mu_A+mu_T+mu_L));

	// <component name="L_0i">
	L_0=(mu_L/(mu_A+mu_T+mu_L));
}