0 nk uninitial delta u number of time steps 1000 -15.d0 5.00d-2 600 0 gamma^2 = - nu^2 = m^2-1/4 8.75d0 0 estimated error initial step size minimum step size 1.d-11 1.d-10 1.d-15