import nsrunit;
unit conversion on;
unit s=1 second^1;
unit per_s=1 second^(-1);

math main {
	realDomain time s;
	time.min=0;
	extern time.max;
	extern time.delta;
	real P(time) dimensionless;
	when(time=time.min) P=0;
	real Q(time) dimensionless;
	when(time=time.min) Q=0;
	real a per_s;
	a=0.1;
	real b per_s;
	b=0.1;
	real K dimensionless;
	K=0.2;
	real P_star dimensionless;
	real V_3 per_s;
	V_3=6;
	real V_4 per_s;
	V_4=2.5;
	real K_3 dimensionless;
	K_3=0.01;
	real K_4 dimensionless;
	K_4=0.01;
	real R(time) dimensionless;
	when(time=time.min) R=0;
	real V_1 per_s;
	V_1=1;
	real V_2 per_s;
	V_2=1.5;
	real K_1 dimensionless;
	K_1=0.01;
	real K_2 dimensionless;
	K_2=0.01;
	real R_star dimensionless;

	// <component name="environment">

	// <component name="P">
	P:time=(a*Q-b*P/(K+P));

	// <component name="P_star">

	// <component name="Q">
	Q:time=(V_1*(1-Q)/(K_1+(1-Q))-V_2*R*Q/(K_2+Q));

	// <component name="R">
	R:time=(P*V_3*(1-R)/(K_3+(1-R))-V_4*R/(K_4+R));

	// <component name="R_star">

	// <component name="law1_Pstar">
	P_star=(V_4/V_3*(1+2*K_3)/(1+2*K_4));

	// <component name="law1_Rstar">
	R_star=(V_1/V_2*(1+2*K_2)/(1+2*K_1));
}