function du = slow(t, u, V, mu)

du = mu*(0.2+V);