We simulated and analyzed Goos–Hänchen (GH) shifts of 633 nm polarized light through a two-prism structure, consisting of a right triangle prism and an isosceles triangle prism with Kretschmann–Raether configuration, by comparing the results from COMSOL Multiphysics (CM) simulation software with that of a stationary-phase analysis (SPA). For this two-prism structure, using a gold film that of thickness 45 nm, the maximum positive GH shift, obtained using SPA at the resonance angle of 44.1°, was 354 μm. Using CM at an incident angle of 43.8°, we found the maximum positive GH shift of 9.45 μm. The results obtained using CM are in agreement with those obtained by the SPA around the resonance angle, although the enhancement effect from CM is much less than that of SPA. This is because SPA depends on the differentiation of the phase shift with respect to the incident angle, while a drastic phase shift occurs at the resonance angle. These results are useful for designing high-sensitivity SPR sensors based on GH shift measurement and for application in waveguide-type SPR devices, with sizes in the order of micro millimeter.