import nsrunit;
unit conversion on;
unit per_second=1 second^(-1);
unit mM=1 meter^(-3)*mole^1;
unit per_mM_per_second=1 meter^3*second^(-1)*mole^(-1);
unit mM_per_second=1 meter^(-3)*second^(-1)*mole^1;

math main {
	realDomain time second;
	time.min=0;
	extern time.max;
	extern time.delta;
	real k_on per_mM_per_second;
	k_on=100000;
	real k_1 per_second;
	k_1=4506534.5;
	real k_2 per_second;
	k_2=16975046.7;
	real k_3o per_second;
	k_3o=100000;
	real k_3i per_second;
	k_3i=531067.4;
	real k_4 per_second;
	k_4=174094707.5;
	real k_fe per_second;
	k_fe=100000;
	real k_be per_second;
	real k_ff per_second;
	k_ff=8747;
	real k_bf per_second;
	k_bf=67210;
	real E_T dimensionless;
	E_T=1;
	real Ce_Na(time) mM;
	when(time=time.min) Ce_Na=0;
	real Ce_K mM;
	Ce_K=10;
	real Ce_Cl mM;
	Ce_Cl=110;
	real Ci_Na mM;
	Ci_Na=7;
	real Ci_K mM;
	Ci_K=124;
	real Ci_Cl mM;
	Ci_Cl=9;
	real E_0(time) dimensionless;
	when(time=time.min) E_0=0.0495937;
	real E_1(time) dimensionless;
	when(time=time.min) E_1=0.10666;
	real E_2(time) dimensionless;
	when(time=time.min) E_2=0.0691;
	real E_3(time) dimensionless;
	when(time=time.min) E_3=0.6873;
	real E_4(time) dimensionless;
	when(time=time.min) E_4=0.04342;
	real E_5(time) dimensionless;
	when(time=time.min) E_5=0.0001297;
	real E_6(time) dimensionless;
	when(time=time.min) E_6=0.0003047;
	real E_7(time) dimensionless;
	when(time=time.min) E_7=0.005334;
	real E_8(time) dimensionless;
	when(time=time.min) E_8=0.0001985;
	real E_9_temp(time) dimensionless;
	when(time=time.min) E_9_temp=0.038;
	real E_9(time) dimensionless;
	real J_23(time) per_second;
	real J_32(time) per_second;
	real J_34(time) per_second;
	real J_43(time) per_second;
	real J_45(time) per_second;
	real J_54(time) per_second;
	real J_56(time) per_second;
	real J_65(time) per_second;
	real J_67(time) per_second;
	real J_76(time) per_second;
	real J_78(time) per_second;
	real J_24(time) per_second;
	real J_42(time) per_second;
	real J_25(time) per_second;
	real J_52(time) per_second;
	real J_26(time) per_second;
	real J_62(time) per_second;
	real J_27(time) per_second;
	real J_72(time) per_second;
	real J_28(time) per_second;
	real J_Rb_norm(time) dimensionless;

	// <component name="environment">

	// <component name="parameters">
	k_be=(k_ff*k_3i*k_fe/(k_bf*k_3o));

	// <component name="ion_concentrations">
	Ce_Na:time=((97 mM_per_second)-(0 mM_per_second));

	// <component name="kinetic_equations">
	E_0:time=((-1)*k_on*Ce_Na*E_0+k_1*E_1-k_be*E_0+k_fe*E_9);
	E_1:time=((-1)*k_on*Ce_Cl*E_1+k_2*E_2-k_1*E_1+k_on*Ce_Na*E_0);
	E_2:time=((-1)*k_on*Ce_K*E_2+k_3o*E_3-k_2*E_2+k_on*Ce_Cl*E_1);
	E_3:time=((-1)*k_on*Ce_Cl*E_3+k_4*E_4-k_3o*E_3+k_on*Ce_K*E_2);
	E_4:time=(k_on*Ce_Cl*E_3-k_4*E_4+k_bf*E_5-k_ff*E_4);
	E_5:time=(k_on*Ci_Na*E_6-k_1*E_5-k_bf*E_5+k_ff*E_4);
	E_6:time=((-1)*k_on*Ci_Na*E_6-k_2*E_6+k_1*E_5+k_on*Ci_Cl*E_7);
	E_7:time=((-1)*k_on*Ci_Cl*E_7+k_2*E_6-k_3i*E_7+k_on*Ci_K*E_8);
	E_8:time=((-1)*k_on*Ci_K*E_8+k_3i*E_7-k_4*E_8+k_on*Ci_Cl*E_9);
	E_9_temp:time=((-1)*k_on*Ci_Cl*E_9+k_4*E_8-k_fe*E_9+k_be*E_0);
	E_9=(E_T-(E_0+E_1+E_2+E_3+E_4+E_5+E_6+E_7+E_8));

	// <component name="flux_equations">
	J_23=(k_on*Ce_K*E_2);
	J_32=(k_3o*E_3);
	J_34=(k_on*Ce_Cl*E_3);
	J_43=(k_4*E_4);
	J_45=(k_ff*E_4);
	J_54=(k_bf*E_5);
	J_56=(k_1*E_5);
	J_65=(k_on*Ci_Na*E_6);
	J_67=(k_2*E_6);
	J_76=(k_on*Ci_Cl*E_7);
	J_78=(k_3i*E_7);
	J_24=(J_23*J_34/(J_34+J_32));
	J_42=(J_32*J_43/(J_34+J_32));
	J_25=(J_24*J_45/(J_45+J_42));
	J_52=(J_42*J_54/(J_45+J_42));
	J_26=(J_25*J_56/(J_56+J_52));
	J_62=(J_52*J_65/(J_56+J_52));
	J_27=(J_26*J_67/(J_67+J_62));
	J_72=(J_62*J_76/(J_67+J_62));
	J_28=(J_27*J_78/(J_78+J_72));
	J_Rb_norm=(J_28/(500 per_second));
}