double precision function myfunc(x, y, func)
        double precision x, y, z, v
        double precision func
        v = 1.5
        z = func(v)
        myfunc = x * y * z
        end