pure function minFunOverB_0_R2(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = -Root2 * funOverB_0(b) end function minFunOverB_0_R2