TITLE One Compartment Model, with Dosing

PARAMETER
{
        decay_time = 5.0         (hr)
        dose = 10
        dose_interval = 10       (hr)
        first_dose_time = 0      (hr)
}

INDEPENDENT
{
        time FROM 0 TO 50 WITH 200      (hr)
}

ASSIGNED
{
        next_dose_time
}

STATE
{
	y FROM 0 TO 12
}

DERIVATIVE rate
{                               
        y' = - y / decay_time
}

INITIAL
{
        next_dose_time = first_dose_time
}

BREAKPOINT
{
        if (time >= next_dose_time) {
                y = y + dose
                next_dose_time = next_dose_time + dose_interval
        }
	SOLVE rate 
}