TITLE Diffusion model, lumped representation
:  National Biomedical Simulation Resource, 3-20-91
:  Ten-segment lumped model of diffusion in a leaky cell.
:  Default output is a series of rising, saturating
:    curves, with maximum amplitude decreasing as the
:    distance from the input end increases.

PARAMETER
{
	kc = 3.		(/sec)
	kl = .3		(/sec)
	delay = 0	(sec)
	height = 1
}

INDEPENDENT
{
	time FROM 0 TO 10 WITH 100	(sec)
}

PLOT C1,C3,C5,C7,CA VS time

STATE
{
	C1 C2 C3 C4 C5 C6 C7 C8 C9 CA
}

KINETIC diffuse
{
~	C1 << (step(time, delay, height))
~	C1 <-> C2	(kc,kc)
~	C1 ->		(kl)
~	C2 <-> C3	(kc,kc)
~	C2 ->		(kl)
~	C3 <-> C4	(kc,kc)
~	C3 ->		(kl)
~	C4 <-> C5	(kc,kc)
~	C4 ->		(kl)
~	C5 <-> C6	(kc,kc)
~	C5 ->		(kl)
~	C6 <-> C7	(kc,kc)
~	C6 ->		(kl)
~	C7 <-> C8	(kc,kc)
~	C7 ->		(kl)
~	C8 <-> C9	(kc,kc)
~	C8 ->		(kl)
~	C9 <-> CA	(kc,kc)
~	C9 ->		(kl)
~	CA ->		(kl)
}

BREAKPOINT
{
	SOLVE diffuse
}