pure function minFunOverB_2(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = -funOverB(b, 2.0_dp) end function minFunOverB_2