A METHOD FOR THE RAPID NUMERICAL CALCULATION OF PARTIAL SUMS OF GENERALIZED HARMONICAL SERIES WITH PRESCRIBED ACCURACY

Similar documents
In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is

Area and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ]

CURVE FITTING LEAST SQUARES METHOD

COMPUTERISED ALGEBRA USED TO CALCULATE X n COST AND SOME COSTS FROM CONVERSIONS OF P-BASE SYSTEM WITH REFERENCES OF P-ADIC NUMBERS FROM

MTH 146 Class 7 Notes

Review of Linear Algebra

Chapter 7. Bounds for weighted sums of Random Variables

DATA FITTING. Intensive Computation 2013/2014. Annalisa Massini

APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:

Sequences and summations

Patterns of Continued Fractions with a Positive Integer as a Gap

Department of Statistics, Dibrugarh University, Dibrugarh, Assam, India. Department of Statistics, G. C. College, Silchar, Assam, India.

Introduction to mathematical Statistics

Fibonacci and Lucas Numbers as Tridiagonal Matrix Determinants

ICS141: Discrete Mathematics for Computer Science I

, we would have a series, designated as + j 1

6.6 Moments and Centers of Mass

Solutions Manual for Polymer Science and Technology Third Edition

19 22 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y Write the given sum or difference as a single integral in

20 23 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y 2 2

Numerical differentiation (Finite difference)

Sebastián Martín Ruiz. Applications of Smarandache Function, and Prime and Coprime Functions

Chapter Linear Regression

SUMMARY OF THE ZETA REGULARIZATION METHOD APPLIED TO THE CALCULATION OF DIVERGENT SERIES

1. (25 points) Use the limit definition of the definite integral and the sum formulas to compute. [1 x + x2

Analytical Approach for the Solution of Thermodynamic Identities with Relativistic General Equation of State in a Mixture of Gases

On Solution of Min-Max Composition Fuzzy Relational Equation

MATRIX AND VECTOR NORMS

Available online through

Chapter Simpson s 1/3 Rule of Integration. ( x)

MULTIOBJECTIVE NONLINEAR FRACTIONAL PROGRAMMING PROBLEMS INVOLVING GENERALIZED d - TYPE-I n -SET FUNCTIONS

European Journal of Mathematics and Computer Science Vol. 3 No. 1, 2016 ISSN ISSN

MATLAB PROGRAM FOR THE NUMERICAL SOLUTION OF DUHAMEL CONVOLUTION INTEGRAL

An Alternative Method to Find the Solution of Zero One Integer Linear Fractional Programming Problem with the Help of -Matrix

On Several Inequalities Deduced Using a Power Series Approach

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

TiCC TR November, Gauss Sums, Partitions and Constant-Value Codes. A.J. van Zanten. TiCC, Tilburg University Tilburg, The Netherlands

An Extended Mixture Inverse Gaussian Distribution

10.2 Series. , we get. which is called an infinite series ( or just a series) and is denoted, for short, by the symbol. i i n

Probabilistic approach to the distribution of primes and to the proof of Legendre and Elliott-Halberstam conjectures VICTOR VOLFSON

this is the indefinite integral Since integration is the reverse of differentiation we can check the previous by [ ]

PATTERNS IN CONTINUED FRACTION EXPANSIONS

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

Regression. By Jugal Kalita Based on Chapter 17 of Chapra and Canale, Numerical Methods for Engineers

14.2 Line Integrals. determines a partition P of the curve by points Pi ( xi, y

CS321. Introduction to Numerical Methods

On the Pell p-circulant sequences

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

PubH 7405: REGRESSION ANALYSIS REGRESSION IN MATRIX TERMS

The definite Riemann integral

Some Unbiased Classes of Estimators of Finite Population Mean

The linear system. The problem: solve

More Regression Lecture Notes CE 311K - McKinney Introduction to Computer Methods Department of Civil Engineering The University of Texas at Austin

COMPLEX NUMBERS AND DE MOIVRE S THEOREM

The z-transform. LTI System description. Prof. Siripong Potisuk

Z b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...

Chapter 2 Intro to Math Techniques for Quantum Mechanics

Preliminary Examinations: Upper V Mathematics Paper 1

Debabrata Dey and Atanu Lahiri

APPROXIMATE INTEGRATION

3. REVIEW OF PROPERTIES OF EIGENVALUES AND EIGENVECTORS

Chapter 12-b Integral Calculus - Extra

Mu Sequences/Series Solutions National Convention 2014

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

Math 2414 Activity 16 (Due by end of class August 13) 1. Let f be a positive, continuous, decreasing function for x 1, and suppose that

THE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK SUMMER EXAMINATION 2005 FIRST ENGINEERING

Chapter Gauss-Seidel Method

Mathematically, integration is just finding the area under a curve from one point to another. It is b

Keywords: Heptic non-homogeneous equation, Pyramidal numbers, Pronic numbers, fourth dimensional figurate numbers.

In an algebraic expression of the form (1), like terms are terms with the same power of the variables (in this case

The Regulated and Riemann Integrals

Abstract inner product spaces

24 Concept of wave function. x 2. Ae is finite everywhere in space.

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Stats & Summary

International Journal of Scientific and Research Publications, Volume 3, Issue 5, May ISSN

POWER SERIES R. E. SHOWALTER

lower lower median upper upper proportions. 1 cm= 10 mm extreme quartile quartile extreme 28mm =?cm I I I

INTRODUCTION TO INTEGRATION

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

ELEG 3143 Probability & Stochastic Process Ch. 5 Elements of Statistics

EECE 301 Signals & Systems

The Z-Transform in DSP Lecture Andreas Spanias

Math 1B, lecture 4: Error bounds for numerical methods

Multiple Choice Test. Chapter Adequacy of Models for Regression

ELEMENTS OF NUMBER THEORY. In the following we will use mainly integers and positive integers. - the set of integers - the set of positive integers

CS321. Numerical Analysis

Integration by Parts for D K

Math 153: Lecture Notes For Chapter 1

Numerical Differentiation and Integration

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

FORMULAE FOR FINITE DIFFERENCE APPROXIMATIONS, QUADRATURE AND LINEAR MULTISTEP METHODS

INTRODUCTION ( ) 1. Errors

On a class of analytic functions defined by Ruscheweyh derivative

Complex Variables. Chapter 19 Series and Residues. March 26, 2013 Lecturer: Shih-Yuan Chen

Bond Additive Modeling 5. Mathematical Properties of the Variable Sum Exdeg Index

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

D KL (P Q) := p i ln p i q i

Systems of second order ordinary differential equations

Transcription:

UPB c Bull, eres D, Vol 8, No, 00 A METHOD FOR THE RAPD NUMERAL ALULATON OF PARTAL UM OF GENERALZED HARMONAL ERE WTH PRERBED AURAY BERBENTE e roue o etodă ouă etru clculul rd l suelor rţle le serlor roce geerlzte, cu o recze restbltă Folosd tegrăr e tervle lese vtjos se obţ exres sle de clcul cre reduc sure cu ulte orde de ăre De seee, etod erte găsre uor exres sle etru rg le erorlor sufcet de strâse Reducere tulu de clcul fţă de rutele extse e clcultorele ctule este sefctvă Oe rooses ew ethod for the rd clculto of rtl sus of geerlzed hrocl seres wth rescrbed ccurcy By usg tegrto o dvtgeously selected tervls oe obts sle exressos of clculus reducg the suto to servl orders of gtude As well, the gve ethod llows fdg sle exressos for error bouds restrcted tervls The reducto of coutto te s cored to the exstg rutes s sgfct Proble forto Let be the su: ( ), R () For fro () oe obtes the hrocl geerlzed seres tht for > s coverget We cll ( ) N rtl su The roble s to ccurte d rd coute ths su for very lrge, rtculr for coverget seres To ths, oe cosders the fucto: Assst, Det of stee eroutce, Uversty POLTEHNA of Buchrest

Berbete f x ( x), x > 0, R, () whch s fte dfferetble The we wrte Tylor forul of the for: f ( x) ( x ) ( x ) ( ) * ( x ) R ( x,, ), N, () where the rest R hs the exresso: ( x ) R ( x,, ) ; ( x) ( x) ( θ) θ x, θ ( 0,); () Further, oe tegrtes the fucto f (x) wth the lts [ -, ], (0;) coveetly chose By deotg ths tegrl, d tkg () to ccout, oe yelds: 5 dx, (, ) x 5 (5) For the rest R, we hve led the e vlue theore O the other hd, by drect tegrto of f(x), oe obts: l, - for ; [( ) ( ) ], for () We wt sle exresso for the su:, The, ccordg to (), oe c get reducto of ters, f there s uber k N, such s: k ; k, (0;) (7)

A ethod for the rd uercl clculto of rtl sus 5 Fro (7), oe gets two ossble vlues: / ; (8) Now we wrte the relto (5) for d : ;, ; 5 (), ; 80 (0) By eltg the ter /, oe gets: 05 0 8 () Now we tke the su of () fro, to,, to obt:, () were stds for the error ter: 0 ; 05 5 () By deotg the su:, () oe gets the forul for clculto of the rtl su, :, (5)

Berbete By usg (), oe yelds: ` ( ) / (( ) ), ( ) /, () the lt, 0, oe gets: () ( ) ( )( ) l l ( ) (-) For 0, oe recover the trvl cse ( 0 ) ( 0, d ) 0 For,, fro () oe obts: ( ) ( ), < 0 ; (-b) (-c) The error evluto ( ) order to fd boud for the error (see ), we cosder, fro () d (0), the ost desdvtgeous cobto of d to get xu ostve vlue: ( ) 5 x ( ) 05, () e the xu dfferece of ters rthess

A ethod for the rd uercl clculto of rtl sus 7 Further, we observe tht for y twce dfferetble fucto g(x), oe obts, by usg Tylor forul d tegto the tervl [, ], the exreso: g ( x) dx g g g" ( ), (, ), e () g e beg the e vlue of g(x) the tervl (, ) The, for g (x) > 0, oe obt the egulty: () g e < g ; g" > 0, tht s, for g > 0, the e vlue s sller th the fucto vlue the ddle of the tervl Now, cosder the fucto: 05 x ( x), x, g > x () whch s ostve together wth ts dervtve g (x) for x > / Tkg the tervl [, ], > /, we fd the ddle the vlue rethess fro (): g ( ) 05 (5) fct, the dex tkes lrger vlues ( 7 ) Thus, we get the error evluto: ( ) g ( x) dx g ( x) x 5 ( )( ) ( G ( ) G ( ) ) 0 5 dx ()

8 Berbete where the fucto G(x) coes fro tegrto hvg the exresso: G 05 x ( x) x (7) By deotg x the lt for : ( ) ( )( ) 0 05 x ( ), (8) oe gets the bsolute error evluto: Results orsos ( ) ( ) ( ) x (8), x Tbles,, re gve the sus for 05, (hrocl seres) d (coverget seres) desred ccurcy of t lest 5; d 7 exct decl fgures resectvely Tble 05 ( 05) ( 05) Tble exct,x 7 0788 05 E - 5 58 085 E - 807 007 E - 7 Exele lculte su 0,5 00000 wth t lest seve exct decl fgures, by usg Aswer Oe clcultes ( 05) by usg the relto (): 00 000 ( 05) 00 000 7008;

A ethod for the rd uercl clculto of rtl sus the: ( 05) ( 05) ( 05) 00000 000000 0758 exct, x 8 77857 08 E - 5 007 0807 E - 0 5777 08 E - 7 Tble Exle By usg Tble, clculte the su of ters of the hrocl seres foud betwee oe llo d bllo, wth t lest seve exct decl fgures Aswer Oe hs to clculte the dfferece (see ): ( ) ( ) 0 0 0 l 0 0 ( 0,exct ( 0 ) ( 0 ) 0,exct 077758 ) 0 l 0 0 Rerk ths cse, we hve dfferece error: d oe c exect eve sller errors th tht,0,0 0 reseted Tble By usg Petu couter t 500 MHz, the coutto te s bout 8 utes Tble exct, x 7 577 05 E - 5 57 08 E - 00 0 E - 7 Exle Becuse for, oe gets coverget seres we c clculte the su (see b; c): 00 007 0

0 Berbete Tble, the xu error vrto for, ( ), x s reseted: fter cresg o the tervl [0; 075], the xu error s decresg costtly Error vrto wth Tble 0 0 05 050 075 0 5, x 0 0 07 0580 085 080 0807 05 087 R E F E R E N E Ncolescu M, Dculescu N, Mrcus, Alză tetcă vol, EDP, Bucureşt, tăăşlă O, Alză tetcă, EDP, Bucureşt, 8 Berbete, Mtr, Zcu, Metode uerce de clcul Ed Tehcă, Bucureşt, 77 Berbete, Berbete, Metode ctttve, vol, Ed Bre, Bucureşt, 00