0.1d0 10.d0 ! a, b 200 ! maxit 1.d-4 ! tol