function dF1dx(x) result(dFdx_val) real(rk), intent(in) :: x real(rk) :: dFdx_val dFdx_val = 15.0_rk * x**2 + 8.0_rk end function dF1dx