function F1(x) result(F_val) real(rk), intent(in) :: x real(rk) :: F_val F_val = 5.0_rk * x**3 + 8.0_rk * x - 5.0_rk end function F1