The Approximate Calculation of Definite Simple and Multiple Integrals

Similar documents
University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

The 7 th Balkan Conference on Operational Research BACOR 05 Constanta, May 2005, Romania

Departure Process from a M/M/m/ Queue

1 cos. where v v sin. Range Equations: for an object that lands at the same height at which it starts. v sin 2 i. t g. and. sin g

where v means the change in velocity, and t is the

Mathematical Economics MEMF e ME. Filomena Garcia. Topic 2 Calculus

NECESSARY AND SUFFICIENT CONDITIONS FOR ALMOST REGULARITY OF UNIFORM BIRKHOFF INTERPOLATION SCHEMES. by Nicolae Crainic

TEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by

MD. LUTFOR RAHMAN 1 AND KALIPADA SEN 2 Abstract

Projectile Motion. Parabolic Motion curved motion in the shape of a parabola. In the y direction, the equation of motion has a t 2.

Goal Programming Approach to Solve Multi- Objective Intuitionistic Fuzzy Non- Linear Programming Models

Confidence intervals for weighted polynomial calibrations

Department of Economics, Niigata Sangyo University, Niigata, Japan

Stat 543 Exam 2 Spring 2016

Stat 543 Exam 2 Spring 2016

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2013

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

ECE 2C, notes set 7: Basic Transistor Circuits; High-Frequency Response

Chapter 3 Differentiation and Integration

1. Statement of the problem

Spring Force and Power

A Satisfaction Degree of Optimal Value for. Grey Linear Programming*

Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS053) p.5185

PHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014

Spectral method for fractional quadratic Riccati differential equation

The gravitational field energy density for symmetrical and asymmetrical systems

Rectilinear motion. Lecture 2: Kinematics of Particles. External motion is known, find force. External forces are known, find motion

Chapter 5. Answers to Even Numbered Problems m kj. 6. (a) 900 J (b) (a) 31.9 J (b) 0 (c) 0 (d) 31.9 J. 10.

On the number of regions in an m-dimensional space cut by n hyperplanes

Chapter 10 Sinusoidal Steady-State Power Calculations

CONTINUOUS PARAMETER FREE FILLED FUNCTION METHOD

Complete Variance Decomposition Methods. Cédric J. Sallaberry

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

Conservation of Energy

Population element: 1 2 N. 1.1 Sampling with Replacement: Hansen-Hurwitz Estimator(HH)

A CLASS OF RECURSIVE SETS. Florentin Smarandache University of New Mexico 200 College Road Gallup, NM 87301, USA

Least Squares Fitting of Data

Lecture 22: Potential Energy

Lecture 09 Systems of Particles and Conservation of Linear Momentum

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

EE 330 Lecture 24. Small Signal Analysis Small Signal Analysis of BJT Amplifier

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES

8.6 The Complex Number System

TREATMENT OF THE TURNING POINT IN ADK-THEORY INCLUDING NON-ZERO INITIAL MOMENTA

Designing of Combined Continuous Lot By Lot Acceptance Sampling Plan

Algorithms for factoring

FUNDAMENTELE PROGRAMĂRII

Physics 2A Chapter 3 HW Solutions

PHYSICS 203-NYA-05 MECHANICS

Momentum. Momentum. Impulse. Momentum and Collisions

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Degrees of Freedom. Spherical (ball & socket) 3 (3 rotation) Two-Angle (universal) 2 (2 rotation)

STAT 511 FINAL EXAM NAME Spring 2001

Endogenous timing in a mixed oligopoly consisting of a single public firm and foreign competitors. Abstract

On the Calderón-Zygmund lemma for Sobolev functions

Numerical Transient Heat Conduction Experiment

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Regularity Properties of p-distance. Transformations in Image Analysis

The Distribution of Energy Within. Baryonic Sub-Atomic Particles.

Chapter 21 - The Kinetic Theory of Gases

A Radon-Nikodym Theorem for Completely Positive Maps

Physics 2A Chapters 6 - Work & Energy Fall 2017

General Results of Local Metric Dimensions of. Edge-Corona of Graphs

9. Complex Numbers. 1. Numbers revisited. 2. Imaginary number i: General form of complex numbers. 3. Manipulation of complex numbers

Using Genetic Algorithms in System Identification

The Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions

+ E 1,1.k + E 2,1.k Again, we need a constraint because our model is over-parameterized. We add the constraint that

Chapter 3 and Chapter 4

10/24/2013. PHY 113 C General Physics I 11 AM 12:15 PM TR Olin 101. Plan for Lecture 17: Review of Chapters 9-13, 15-16

COS 511: Theoretical Machine Learning

FIVE INTEGER NUMBER ALGORITHMS TO SOLVE LINEAR SYSTEMS

Osp(1 2M)-invariant higher-spin systems

Applied Stochastic Processes

STATISTICS QUESTIONS. Step by Step Solutions.

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

The path of ants Dragos Crisan, Andrei Petridean, 11 th grade. Colegiul National "Emil Racovita", Cluj-Napoca

{ } MATH section 7.2 Volumes (Washer Method) Page 1 of 8. = = 5 ; about x-axis. y x y x. r i. x 5= interval: 0. = x + 0 = x + = + = +

AN ALGORITHM FOR CALCULATING THE CYCLETIME AND GREENTIMES FOR A SIGNALIZED INTERSECTION

Gaussian Mixture Models

total If no external forces act, the total linear momentum of the system is conserved. This occurs in collisions and explosions.

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

The middle point of each range is used to calculated the sample mean and sample variance as follows:

Constraints on multiparticle production in scalar field theory from classical simulations

International Journal of Mathematical Archive-4(12), 2013, Available online through ISSN

a. (All your answers should be in the letter!

EE215 FUNDAMENTALS OF ELECTRICAL ENGINEERING

A NOTE ON THE DISCRETE FOURIER RESTRICTION PROBLEM

Physics 4C. Chapter 19: Conceptual Questions: 6, 8, 10 Problems: 3, 13, 24, 31, 35, 48, 53, 63, 65, 78, 87

(b) i(t) for t 0. (c) υ 1 (t) and υ 2 (t) for t 0. Solution: υ 2 (0 ) = I 0 R 1 = = 10 V. υ 1 (0 ) = 0. (Given).

Subspace Algorithms in Modal Parameter Estimation for Operational Modal Analysis: Perspectives and Practices

Computational and Statistical Learning theory Assignment 4

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Physical Chemistry I for Biochemists. Lecture 18 (2/23/11) Announcement

Practical Newton s Method

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Bayesian decision theory. Nuno Vasconcelos ECE Department, UCSD

Complex Numbers Alpha, Round 1 Test #123

Math1110 (Spring 2009) Prelim 3 - Solutions

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Transcription:

BULETNUL Unverstăţ Petrol Gae dn Ploeşt Vol. LV No. /6 7 - Sera Mateatcă - noratcă - Fcă The Aroate Calculaton o Dente Sle and Multle nterals Tănase Dnu Unverstatea Petrol-Gae dn Ploeşt Bd. Bucureşt 9 Ploeşt Catedra de Mateatcă e-al: tdnu@al.u-loest.ro Abstract Ths aer resents the local scheata and the roras n turbo-ascal or the aroate calculaton ethods o the sle dente nteral [] o the double dente nteral [] and o the trle dente nteral []. At the end there are analed the calculaton errors corresondn to these ethods that use the olnoals o nterolaton to ve nodes. Ke words: dvson nterolaton alorth ntroducton Paer [] resents ethod B or the aroate calculaton o sle dente nteral where :[ a b] R s a nterable uncton on [ ] t s consdered an equdstant dvson wth Δ : wth the dstance o the dvson densons: b ( ) d () a a b. n nodes. ( a < <... < < < < b)... n b a h then the aroate value o the nteral n n n h 7( ( a) ( b) ) ( ) ( ) 5 n n ( ) ( ) () or whch the calculaton error s: where: e M 5 ( b a) 6 τ h ()

8 Tănase Dnu su ( 5) ( ) 5 M cu C [ ab ]. [ a b] n aer [] the BB ethod descrbes the aroate calculaton o the double dente nteral: wth : D R nterable where a b. wth ( ) ( ) or an [ ] D ( ) dd () D R and t s bordered b both ( ) The aroate value o the nteral s ven b the orula: n n n h 7 n 5 where: ( ϕ ( ) ϕ( )) ϕ( ) ϕ( ) ϕ( ) wth ( ) n ( ) ϕ (5) h ϕ( ) ( ( ) ( )) ( ) ( ) 7 5 ( ) ( ) (6) b a where the nor o the nterval [a b] s: h and the nor o the nterval [ ( ) ( )] n ( ) ( ) s: h wth {... n}. Paer [] resents the ethod BBB or the aroate calculaton o the trle dente nterals. wth : V R nterable on V ( ) ddd (7) V R. The doan V s dened b suraces: ( ) and ( ) wth ( ) ( ) or ( )( ) D where D s the roecton o the doan V on the O lan and t s dened b the curves: ( ) and ( ) wth ( ) ( ) or ( ) [ a b] where [a b] s the roecton o D on the as O. The aroate value o the nteral s ven b the orula: n n ( ( ) ( )) ( ) ( ) ( ) n h 7 ϕ ϕ n ϕ ϕ ϕ 5 n ( ) ϕ (8) b a h ν Δ and n where the nterval [a b] has the nor ( )

The Aroate Calculaton o Dente Sle and Multle nterals 9 ( ) ( ) ( ) ( ) ( ) ( ) 7 5 h ϕ ( ) ( ) (9) where the nterval ( ) ( ) [ ] has the nor ( ) ( ) ( ) h Δ ν and ( ) ( ) ( ) ( ) ( ) 7 5 h ( ) ( ) ( ) () or { } { } n...... and the nterval ( ) ( ) [ ] has the nora o the equdstant dvson ( ) ( ) ( ) h Δ ν. Contents A. Proras rora sle_nteral; te vecarra[..5] o real; var a:vec; absssht:real; n:nteer; uncton su(var a:vec; :nteer; h:real):real; var :nteer; ssssss:real; ben s:; s:; s:; s:; or : to - do ben s:sa[*]; s:sa[*]; s:sa[*]; s:sa[*] s:s-a[]; ss:7*(a[] a[*n]); ss:ss*s*s*s*s; ss:(ss**h)/5; su:ss

Tănase Dnu uncton un(var :nteer; ah:real):real; var :real; ben :a*h; un:(*)/(**) ben {clrscr;} wrte('a'); readln(a); wrte('b'); readln(b); wrte('n'); readln(n); h:(b-a)/(*n); or : to *n do a[]:un(ah); s:su(anh); wrteln('nteral 's::6); {reeat untl eressed} end. rora double_nteral; te vecarra[..5] o real; var c:vec; absshabh:real; n:nteer; uncton su(var a:vec;:nteer;h:real):real; var :nteer; ssssss:real; ben s:; s:; s:; s:; or : to - do ben s:sa[*]; s:sa[*]; s:sa[*]; s:sa[*] s:s-a[]; ss:7*(a[] a[*n]); ss:ss*s*s*s*s; ss:(ss**h)/5; su:ss uncton un(:nteer; ahah:real):real; var :real; ben :a*h; :a*h; un:/sqrt(**)

The Aroate Calculaton o Dente Sle and Multle nterals uncton (var :nteer; ah:real):real; var t:real; ben :a*h; t:sqrt(abs(-*)); :t; uncton (var :nteer; ah:real):real; var t:real; ben :a*h; t:-sqrt(abs(-*)); :t ben {clrscr; } wrte('a'); readln(a); wrte('b'); readln(b); wrte('n'); readln(n); wrte(''); readln(); h:(b-a)/(*n); or : to *n do ben a:(ah); b:(ah); h:(b-a)/(*); or : to * do c[]:un(ahah); []:su(ch); ss:su(nh); wrteln('nteral ' ss::6) end. rora treble_nteral; {uses crt;} te vecarra[..5] o real; var c:vec; abhabhabhss:real; n:nteer; uncton su(var a:vec;:nteer;h:real):real; var :nteer; ssssss:real; ben s:; s:; s:; s:; or : to - do ben s:sa[*]; s:sa[*];

Tănase Dnu s:sa[*]; s:sa[*] s:s-a[]; ss:7*(a[] a[*n]); ss:ss*s*s*s*s; ss:(ss**h)/5; su:ss; uncton un(var :nteer; ahahah:real):real; var t:real; ben :a*h; :a*h; :a*h; t:sqrt(***); un:t uncton (var :nteer; ab:real):real; var :real; ben :a*b; :sqrt(abs(9-*)) uncton (var :nteer; ab:real):real; var :real; ben :a*b; :-sqrt(abs(9-*)) uncton (var :nteer; ahah:real):real; var t:real; ben :a*h; :a*h; t:; :t uncton (var :nteer; ahah:real):real; var :real; ben :a*h; :a*h; :(**)/9 ben {clrscr;} wrte('a'); readln(a); wrte('b'); readln(b); wrte('n'); readln(n);

The Aroate Calculaton o Dente Sle and Multle nterals wrte(''); readln(); wrte(''); readln(); h:(b-a)/(*n); or : to *n do ben a:(ah); b:(ah); h:(b-a)/(*); or : to * do ben a:(abab); b:(abab); h:(b-a)/(*); or : to * do c[]:un(ahahah); []:su(ch); []:su(h); ss:su(nh); wrteln('nteral ' ss::6) end. B. Alcatons Let us calculate the ollown dente nterals: D ( ) dd ; d where D ( ) { 9}; ddd where V ( ) 9 and [ ] V usn the anteror roras. Soluton Drectl calculated b eans o the analtcal ethods we obtan: ln 96; 58599 { } 988 Aroatel calculated b eans o B BB BBB ethods and usn the roraes resented above we obtan: 9 585985 988 or n a b ; or n a b ; or n a b.

Tănase Dnu Reerences. Dnu T. - Anală nuercă Edtura U.P.G. dn Ploeşt. D n u T. - Calculul aroatv al nteralelor sle Buletnul U.P.G. dn Ploeşt nr./. D n u T. - Calculul aroatv al nteralelor duble Buletnul U.P.G. dn Ploeşt nr./5. D n u T. - The Rouh Calculus o Treble Dened nteral Buletnul U.P.G. dn Ploeşt nr./6. -6 5. Grore Gh. - Lecţ de Anală nuercă Facultatea de Mateatcă Bucureşt (curs ltoraat) 6. Marnescu Gh. - Anală nuercă Edtura Acadee Roâne Bucureşt 97 Calculul aroatv al nteralelor dente sle ş ultle Reuat În această lucrare sunt reentate scheele loce ş roraele în Turbo-Pascal entru etodele de calcul aroatv ale nteralelor dente sle [] ale nteralelor dente duble [] ş ale nteralelor dente trle []. În nal sunt analate erorle de calcul coresunătoare acestor etode care olosesc olnoaele de nterolare cu cnc nodur