double precision function f(x) implicit double precision (a-h, o-z) f = sin(x) return end