df(t)/dt = p * (f(t) + A)/(2(t+N)) + (1 - p) f(t)/(2t)