viewbarcode.com

Series and Their Convergence in Software Maker PDF-417 2d barcode in Software Series and Their Convergence




How to generate, print barcode using .NET, Java sdk library control with example project source code free download:
5.3 Series and Their Convergence using none tomake none in asp.net web,windows application Java programming language can sometimes displ none for none ay convergence that is only apparent, not real. The remedy is to try two different transformations as a check. Since convergence acceleration is so much more dif cult for a series of positive terms than for an alternating series, occasionally it is useful to convert a series of positive terms into an alternating series.

Van Wijngaarden has given a transformation for accomplishing this [6]:. 1 X rD1 vr D 1 X rD1 . 1/r (5.3.20).

where wr vr C 2v2 none none r C 4v4r C 8v8r C (5.3.21) Equations (5.

3.20) and (5.3.

21) replace a simple sum by a two-dimensional sum, each term in (5.3.20) being itself an in nite sum (5.

3.21). This may seem a strange way to save on work! Since, however, the indices in (5.

3.21) increase tremendously rapidly, as powers of 2, it often requires only a few terms to converge (5.3.

21) to extraordinary accuracy. You do, however, need to be able to compute the vr s ef ciently for random values r. The standard updating tricks for sequential r s, mentioned above following equation (5.

3.1), can t be used. Once you ve generated the alternating series by Van Wijngaarden s transformation, the Levin d transformation is particularly effective at summing the series [8].

This strategy is most useful for linearly convergent series with close to 1. For logarithmically convergent series, even the transformed series (5.3.

21) is often too slowly convergent to be useful numerically. As an example of how to call the routines Epsalg or Levin, consider the problem of evaluating the integral Z 1 x J0 .x/ dx D K0 .

1/ D 0:4210244382 : : : (5.3.22) I D 1 C x2 0 Standard quadrature methods such as qromo fail because the integrand has a long oscillatory tail, giving alternating positive and negative contributions that tend to cancel.

A good way of evaluating such an integral is to split it into a sum of integrals between successive zeros of J0 .x/: Z I D. f .x/ dx D 1 X j D0 (5.3.23).

where Z Ij D f .x/ dx;. f .xj / D 0;. j D 0; 1; : : :. (5.3.24).

We take x 1 equal t o the lower limit of the integral, zero in this example. The idea is to evaluate the relatively simple integrals Ij by qromb or Gaussian quadrature, and then accelerate the convergence of the series (5.3.

23), since we expect the contributions to alternate in sign. For the example (5.3.

22), we don t even need accurate values of the zeros of J0 .x/. It is good enough to take xj D .

j C 1/ , which is asymptotically correct. Here is the code:. levex.h 5. Evaluation of Functions Doub func(const Dou none none b x) Integrand for (5.3.22).

{ if (x == 0.0) return 0.0; else { Bessel bess; return x*bess.

jnu(0.0,x)/(1.0+x*x); } } Int main_levex(void) This sample program shows how to use the Levin u transformation to evaluate an oscillatory integral, equation (5.

3.22). { const Doub PI=3.

141592653589793; Int nterm=12; Doub beta=1.0,a=0.0,b=0.

0,sum=0.0; Levin series(100,0.0); cout << setw(5) << "N" << setw(19) << "Sum (direct)" << setw(21) << "Sum (Levin)" << endl; for (Int n=0; n<=nterm; n++) { b+=PI; Doub s=qromb(func,a,b,1.

e-8); a=b; sum+=s; Doub omega=(beta+n)*s; Use u transformation. Doub ans=series.next(sum,omega,beta); cout << setw(5) << n << fixed << setprecision(14) << setw(21) << sum << setw(21) << ans << endl; } return 0; }.

Setting eps to 1 10 8 in qromb, we get 9 signi cant digits with about 200 function evaluations by n D 8. Replacing qromb with a Gaussian quadrature routine cuts the number of function evaluations in half. Note that n D 8 corresponds to an upper limit in the integral of 9 , where the amplitude of the integrand is still of order 10 2 .

This shows the remarkable power of convergence acceleration. (For more on oscillatory integrals, see 13.9.

). CITED REFERENCES AN none for none D FURTHER READING: Weniger, E.J. 1989, Nonlinear Sequence Transformations for the Acceleration of Convergence and the Summation of Divergent Series, Computer Physics Reports, vol.

10, pp. 189 371.[1] Abramowitz, M.

, and Stegun, I.A. 1964, Handbook of Mathematical Functions (Washington: National Bureau of Standards); reprinted 1968 (New York: Dover); online at http://www.

nr. com/aands, 3.6.

[2] Mathews, J., and Walker, R.L.

1970, Mathematical Methods of Physics, 2nd ed. (Reading, MA: W.A.

Benjamin/Addison-Wesley), 2.3.[3] Numerical Recipes Software 2007, Derivation of the Levin Transformation, Numerical Recipes Webnote No.

6, at http://www.nr.com/webnotes 6 [4] Smith, D.

A., and Ford, W.F.

1982, Numerical Comparisons of Nonlinear Convergence Accelerators, Mathematics of Computation, vol. 38, pp. 481 499.

[5] Goodwin, E.T. (ed.

) 1961, Modern Computing Methods, 2nd ed. (New York: Philosophical Library), 13 [van Wijngaarden s transformations].[6].

Copyright © viewbarcode.com . All rights reserved.