tinit = 7.d0 tstep = 0.2d0 nt = 21 rinit = 0.d0 rstep = 0.2 nr = 61 n = 2