a = 12.d0 z = 6.d0 tinit = 7.d0 tstep = 1d0 nt = 6 rinit = 0.d0 rstep = 1.d0 nr = 13 n = 2