We employ the method of comparison equations to study the propagation of a massless minimally coupled scalar field on the Schwarzschild background. In particular, we show that this method allows us to obtain explicit approximate expressions for the radial modes with energy below the peak of the effective potential which are fairly accurate over the whole region outside the horizon. This case can be of particular interest, for example, for the problem of black hole evaporation.

We apply the method of comparison equations to study cosmological perturbations during inflation, obtaining the full power spectra of scalar and tensor perturbations to first and to second order in the slow-roll parameters. We compare our results with those derived by means of other methods, in particular the Green's function method and the improved WKB approximation, and find agreement for the slow-roll structure. The method of comparison equations, just as the improved WKB approximation, can however be applied to more general situations where the slow-roll approximation fails.

We extend the WKB method for the computation of cosmological perturbations during inflation beyond leading order and provide the power spectra of scalar and tensor perturbations to second order in the slow-roll parameters. Our method does not require that the slow-roll parameters be constant. Although leading and next-to-leading results in the slow-roll parameters depend on the approximation technique used in the computation, we find that the inflationary theoretical predictions obtained may reach the accuracy required by planned observations. In two technical appendices, we compare our techniques and results with previous findings.

We study the WKB approximation beyond leading order for cosmological perturbations during inflation. To first order in the slow-roll parameters, we show that an improved WKB approximation leads to analytical results agreeing to within 0.1% with the standard slow-roll results. Moreover, the leading WKB approximation to second order in the slow-roll parameters leads to analytical predictions in qualitative agreement with those obtained by the Green's function method.

Improved Wentzel-Kramers-Brillouin (WKB)-type approximations are presented in order to study cosmological perturbations beyond the lowest order. Our methods are based on functions which approximate the true perturbation modes over the complete range of the independent (Langer) variable, from sub-horizon to super-horizon scales, and include the region near the turning point. We employ both a perturbative Green's function technique and an adiabatic (or ``semiclassical'') expansion (for a linear turning point) in order to compute higher order corrections. Improved general expressions for the WKB scalar and tensor power spectra are derived for both techniques. We test our methods on the benchmark of power-law inflation, which allows comparison with exact expressions for the perturbations, and find that the next-to-leading order adiabatic expansion yields the amplitude of the power spectra with excellent accuracy, whereas the next-to-leading order with the perturbative Green's function method does not improve the leading order result significantly. However, in more general cases, either or both methods may be useful.