(b)

> CONC := S(t) / VOL;

> MODEL := diff( S(t), t ) = 1 * 3 - 1 * CONC;

>

> IC := S(0) = 0.4 * 15;