out1(n) = out1(n)/nd/dt*2.0d0