C(t+delt) = ((C(t) - (P(t)/L(t)))*exp(-L(t)*delt)) + (P(t)/L(t))

where we are trying to solve:

dC/dt = P - L*C