0 delt nt phi0 delphi0 1.d-1 10000 1.d-1 1.d-5 0 eps h1 hmin 1.d-14 1.d-10 1.d-16