function func1(x) result(f) complex(rk), intent(in) :: x complex(rk) :: f f = x**2 + 2.0_rk*x end function func1