double precision function f(x) implicit double precision (a-h, o-z) f = dsqrt(dabs(x-2.345d-01)) return end