pure function minFunD_1(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = -funD_1(b) end function minFunD_1