When I do:
REAL FUNCTION f(x) REAL, INTENT(IN) :: x f = exp(-x) END FUNCTION f
it works but if I do:
Use x ** 2 to exponentiate. Or x*x in this case.
x ** 2
x*x