import nsrunit;
unit conversion on;
unit s=1 second^1;
unit nM=1E-6 meter^(-3)*mole^1;
unit per_s=1 second^(-1);
unit nM_per_s=1E-6 meter^(-3)*second^(-1)*mole^1;

math main {
	//Warning:  the following variables were set 'extern' or given
	//  an initial value of '0' because the model would otherwise be
	//  underdetermined:  t
	extern real t s;
	real stimEnd s;
	stimEnd=1800;
	real Ca nM;

	// <component name="calcium">
	Ca=(if ((t>=(0 s)) and (t<=(180 s))) ((459 nM)-(576 nM))*t/(180 s)+(576 nM) else if ((t>(180 s)) and (t<=(360 s))) ((406 nM)-(459 nM))*(t-(180 s))/(180 s)+(459 nM) else if ((t>(360 s)) and (t<=(540 s))) ((377 nM)-(406 nM))*(t-(360 s))/(180 s)+(406 nM) else if ((t>(540 s)) and (t<=(720 s))) ((348 nM)-(377 nM))*(t-(540 s))/(180 s)+(377 nM) else if ((t>(720 s)) and (t<=(900 s))) ((336 nM)-(348 nM))*(t-(720 s))/(180 s)+(348 nM) else if ((t>(900 s)) and (t<=(1080 s))) ((334 nM)-(336 nM))*(t-(900 s))/(180 s)+(336 nM) else if ((t>(1080 s)) and (t<=stimEnd)) (334 nM) else (60 nM));
}