function func5(x) result(f) real(rk), dimension(:), intent(in) :: x real(rk) :: f f = x(1)**2 + 0.5_rk*x(2)**2 end function func5