next up previous contents
Next: About this document Up: Notes on numerical integration Previous: Richardson Extrapolation

Change of variables

 

Change of variables is an extremely powerful method for performing integrals not only analytically but also numerically. A change of variables can considerably improve the accuracy of regular-interval techniques for functions with rapid variations in particular regions of the integration domain and can allow one to perform integrals which would otherwise be impossible, such as for functions with integrable singularities or integrals over infinite domains.

Change of variables begins with the familiar analytic formula,

  equation427

where r(x) is any function which maps the unit interval tex2html_wrap_inline1029 to the interval tex2html_wrap_inline1031 . The most trivial example of such a map is the linear map

eqnarray434

which transforms integration over any interval to integration over the unit interval. Figure 4(a) illustrates such a map along with the distribution of resulting sample points tex2html_wrap_inline1033 .

   figure439
Figure 4: Mapping of regular sampling of unit interval under change of variables: (a) linear map, (b) map to semi-infinite interval, (c) map to semi-infinite interval with concentration of points at origin.

Maps which transform integration over an infinite interval to integration over a finite interval represent a far more important example. The map

  eqnarray451

for instance, makes possible the evaluation of integrals over the interval tex2html_wrap_inline1035 with the techniques discussed above. As tex2html_wrap_inline1037 , tex2html_wrap_inline1039 , and this map spaces the sample points tex2html_wrap_inline1033 regularly near the origin. However, As tex2html_wrap_inline1043 , however, the sample points eventually tend toward infinity. (See Figure 4(b).) Under this map, the integral

displaymath1045

becomes

displaymath1047

which Figure 5 illustrates to be much more manageable. In practice, such changes of variables require care to ensure that the proper limiting values for f(r(x)) (dr/dx) are used at the endpoints where the change of variables factor tex2html_wrap_inline1051 becomes singular.

   figure476
Figure: Integrand of tex2html_wrap_inline1053 : (a) direct integration, (b) under change of variables in Eq. (21).

Finally, as an example of a problem with an integrable singularity consider the integral

displaymath1055

As tex2html_wrap_inline1057 , the integrand approaches tex2html_wrap_inline1059 and diverges toward infinity. (See Figure 6(a).) Heuristically, the solution is to concentrate sample points tex2html_wrap_inline1033 near the singularity at the origin. A map which both extends the range of integration to infinity and concentrates points near the origin is

  eqnarray493

(See Figure 4(c).) Figure 6(b) shows that the resulting integrand f(r(x)) (dr(x)/dx) is now well behaved and quite suited to regular interval techniques.

   figure506
Figure: Integrand of tex2html_wrap_inline1065 : (a) direct integration, (b) under change of variables in Eq. (22).


next up previous contents
Next: About this document Up: Notes on numerical integration Previous: Richardson Extrapolation

Tomas Arias
Sun Mar 5 13:48:24 EST 2000