2 ndim 1.0e0 gravitational constant 0.e0 0.e0 0.5d0 x0,y0,r0 2.e0 1.e0 hin, hout