The phase is nonnegative on the interval and vanishes only at its two endpoints. An interval away from both endpoints contributes exponentially little, so both endpoint contributions must be expanded before their orders are ranked.
At the left endpoint,
The Laplace's method scale is . Hence its leading contribution is
The substitution evaluates the integral by the Gamma function. Since there is no quartic term in the phase, the next left-end correction is relatively , and thus absolutely .
At the right endpoint put . There
so gives
The leading right-end integral is a Gaussian integral. Its first correction is also absolutely . The mixed cubic and quadratic Laplace endpoints therefore produce the first two global terms
Keeping only the dominant cubic endpoint would miss the second term, which is larger than that endpoint's own first correction.