tinit = 7.d0 tfin = 1.1d1 tstep = 1.0d0 rinit = 0.d0 rfin = 1.4d1 rstep = 0.1d0 n = 2.d0