CC BY
0EURASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES
Innovative Academy Research Support Center UIF = 8.3 | SJIF = 7.906 www.in-academy.uz
APPROXIMATE CALCULATION OF THE EXACT INTEGRAL
IRODAKHON RAJABOVA
Nukus State Pedagogical Institute named after Ajiniyaz, 3rd year
student SUVHANOV AG'AJON Nukus State Pedagogical Institute named after Ajiniyaz, 3rd year
student
https://doi.org/10.5281/zenodo.10976795
ABSTRACT
ARTICLE INFO
Received: 08th April 2024 Accepted: 15th April 2024 Online: 16th April 2024
KEYWORDS Infinitesimal, area, volume of a solid, integral, indefinite integral, surface, surface integral, sign of indefinite integral.
The article provides detailed information about the concept of integral and the history of its origin. In addition to the Riemann integral, which is covered in the school program, a number of other integrals (Riemann-Stiltes, Daniell, Alfred Haar, Henstock-Kurzweil, Ito and Stratonovich, Chokvet) are described. The difference between integrals is explained with the help of concrete examples. The concept of definite integral and its application in practical problems are widely covered, counterexamples of some theorems are presented.
XU
/(x) dx in the above-mentioned methods of calculating the definite integral consisted in finding an initial function F(x) of the function f(x) and calculating its value. But for some definite integrals, we may encounter the following problems when applying these methods:
1) Finding the initial function F(x) is complicated;
2) the initial function F(x) is complex, and its F(a) and F(b) values are difficult to calculate;
3) the function F(x) is not expressed in elementary functions;
4) The function f(x) under the integral is given in table form.
In such cases, it is necessary to calculate the exact integral approximately. Various formulas have been found to solve this problem, and they are generally called quadrature formulas. Below are some of these formulas.
1. Formula of a rectangle. To derive this formula, we first divide the section [a,b] into n equal parts with points a = x0 < xx < x2 <... < xn = b In this case, the length of each piece is
equal to Ax = (Figure 1)
We denote the values of the function f(x) under the integral at the points x0, xx, x2,... xn-1, xn by y0, ylt y2,... yn-i, yn and the following we make sums:
yo • Ax + yx • Ax + y2 • Ax + - + yn-iAx; yoAx + yiAx + y2Ax + + ■■■ + ynAx.yo •
Ax + y1 • Ax + y2 • Ax + —+ yn-1Ax; y0Ax + y1Ax + y2Ax + +-----+ ynAx. y0 • Ax + y1 •
Ax + y2 • Ax + —+ yn-1Ax; y0Ax + y1Ax + y2Ax + +-----+ ynAx.
EURASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES
Innovative Academy Research Support Center UIF = 8.3 | SJIF = 7.906 www.in-academy.uz
Figure 1
Each of these sums is an integral sum for the function f(x) in the section [a,b], and the following approximate formulas can be written for them:
I
f(x)dx
b
a
u
I
f(x)dx
n
b — a
( Ï0 +yi +72 + ••• + yn-i)
n
■( yi+y2+y3 + - + yn)
These formulas are called rectangular formulas.
II. Formula of trapezoids. In the formula for rectangles considered above, we replaced the curve y=f(x) with stair lines. If we replace y=f(x) with dashed lines inside, we get a more accurate value of the definite integral. In this case, the curved trapezoid aABb consists of a set of trapezoids bounded from above by AA1,A1A2, ...,An-1B vectors (Figure 2). The first in this
Figure 2
<17 J.17 IT I ^ 7
the face of a trapezoid ——1 • Ax, the other's face ——2
To' 1
•Ax and so on j^f(x)dx
(y°+y±^x + yi+yi
\ -) o
Ax + • b
+ yn-1+yn • ax) or
^ b —a fyo+yn \ J f(x)dx « - (—--+ yi+y2 + - + yn-i)
n \ 2
will be. This formula is called the trapezoidal formula for the approximation of the definite integral. Here the number n is chosen arbitrarily. The larger the number n, the more accurate the value of the integral.
III. Parabola formula (Simpson's formula). We divide the section [a,b] into n=2m equal parts. [x0,XjJ and [xx,x2] are the faces of the curved trapezoids bounded by the curve y=f(x)
a
a
EURASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES
Innovative Academy Research Support Center UIF = 8.3 | SJIF = 7.906 www.in-academy.uz
from the points M0(x0,y0), M1(x1,y1), M2(x2,y2) we replace it with a curved trapezoid bounded by a passing parabola. Such a curved trapezoid is called a parabolic trapezoid (Figure 3).
The equation of a parabola whose axis is parallel to the 0y axis It consists of y = Ax2 + 5x + C.
- 1
y \
,v/
Mi 3
0 = = a = b
Figure 3
Coefficients A, B, C are found from the condition that the parabola passes through the given three points. We make parabolas as above for the rest of the sections. The sum of the faces of the formed parabolic trapezoids gives an approximate value of the integral. It consists of the following formula:
Jai,/(x)dx « ^ [yo + y2rn + 2(y2 + y4 + - + y2m-2) + 4(^1 + 73 + + ■■■ + y2m-l)L This formula is called Simpson's formula.
If the limit on the right side of the above equation exists and is finite, then the eigenintegral is called convergent, otherwise, it is called divergent.
In many cases, it is not necessary to know the exact value of the eigenintegral, it is enough to estimate whether it is converging or receding, and if it is converging, its value is sufficient.
I. If a < x < ro in the infinite semiinterval 0 < /(x) < $(x) and Ja+ro^(x) dx is an eigenintegral approximant, then Ja+ro/(x)dx The characteristic integral is also converging and the following inequality holds:
+œ
+œ
J/WdxsJffWdx
If the function f(x) is defined in the interval its characteristic integral of type I
over this interval is expressed as follows by the characteristic integrals introduced above.
+œ
+œ
J /(x) dx = J /(x) dx + J /(x) dx = lim J /(x) dx +
-œ — œ c a
b
//m* .
+ lim
Ô^+œ
a
a
c
c
c
EURASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES
Innovative Academy Research Support Center UIF = 8.3 | SJIF = 7.906 www.in-academy.uz
Now we generalize the concept of definite integral for unbounded functions. The given function y=f(x) is not bounded on the half-interval (a;b), but for arbitrary se(0,b-a], let this function be bounded and integrable on the section [a+8,b]. In this case, f(ex')f(e) —
fa+J(x) dx, se(0, b — a] function can be viewed.
If the function y=f(x) is not bounded at the point x=a, then [a;+ro) or (-^;a] is determined by the following mixed-type proper integrals over infinite semiintervals:
+œ
+œ
j f(x) dx = j f(x) dx + j f(x) dx (a < b < +œ)
a a b
a c a
j f(x) dx = j f(x) dx + j f(x) dx (— ™ < c < a)
Examples:
-1 dx
1. J0 YT^i' calculate the integral using the rectangular formula.
Solving. For this, we divide the integration section [0;1] into n=10 parts and present the results of calculations in the following table:
b
— œ
— œ
C
i xi = 0.1i 1 + x¡ 1 f(xô = 1 + x? ^f(xt)
1 0.1 1.01 0.9901 0.9901
2 0.2 1.04 0.9615 1.9516
3 0.3 1.09 0.9174 2.8690
4 0.4 1.16 0.8621 3.7311
5 0.5 1.25 0.8000 4.5311
6 0.6 1.36 0.7353 5.2664
7 0.7 1.49 0.6711 5.9375
8 0.8 1.64 0.6098 6.5473
9 0.9 1.81 0.5525 7.0998
10 1.0 2 0.5000 7.5998
l — 0
Since Ax = = 0,1 in our example, we get the following result based on the formula of a right rectangle.
EURASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES
Innovative Academy Research Support Center UIF = 8.3 | SJIF = 7.906 www.in-academy.uz
1 dx
0,1 • 7,5998 = 0,75998.
f
Jo i+x2
2. J^tt^t approximate the integral using the trapezoidal formula. Also take n=10.
1+x2
c i ■ ö—a 1—0 1 n t
Solving.-= — = — = 0,1;
ö n 10 10
yp+yn 2
1+0.5 2
= 0,75 va y1 + y2 + +y3 + - + y9 = 7,0998
y1 + y2 + +y3 + —+ y9 = 7,0998 because
-1 dx
i
0 1+x2 9
0,1 • [0.75 + 7,0998] = 0,78498.
3. J V6x- 5 dx calculate the integral using Simpson's formula. In this case, take n=8.
a 9—1
Solving. h =-= — = 1. calculate the integral using Simpson's formula. In this case,
take n=8.
x0 = 1,y0 = V6 - 1- 5 = VÎ = 1,0000, X1 = 2,y1 = V6- 2 - 5 = V7 « 2.6458, x2 = 3,y2 = V6- 3 - 5 = VÏ3 « 3.6056, x3 = 4,y3 = V6 • 4 — 5 = VÎ9 « 4.3589, x4 = 5,y4 = V6^ 5 - 5 = V25 = 5.0000,
5.5678, « 6.0828,
x7 = 8,y7 = V6 • 8 - 5 = V43 « 6.5574, x8 = 9,y8 = V6^ 9 - 5 = V49 = 7.0000.
x5 = 6, y5 = V6 • 6 - 5 = V3T x6 = 7,y6 = V6 • 7 - 5 = V37
These values found
¿Vto dx « ^ [y0 + y2rn + 2(y2 + y4 + - + y2m—2) + 4(y1 + y3 + +- + y2m—1)]
fc—a
These values found
We have 2m = 8, ¿ — a = 9 — 1 = 8 because ^ = — = 1 So, J9V6x — 5 dx
6m 24 3 J1
1
i[1 + 7 + 2(3.6056 + 5.0000 + 6.0828) + +4(2.6458 + 4.3589 + 5.5678 + 6.5574)] =
1
1(8 + 2 • 14.6884 + 4 • • 19.1299) = 37.9655. integral JQ+ro e—x dx .
4.Calculate
the
characteristic
Solving.
+œ
J
e x dx = lim
fc^+œ
J
dx = - lim (-e x) | = = - lim (e x) |
fc^+œ '0 fc^+œ '0
0
1
J+œe Xdx = lim ifc e Xdx = - lim (-e x) = = - lim (e x) =
= - lim e—+ lim e0 = - lim — + 1 = 1.
fc^+œ fc^+œ fc^+œ
rb , ... __ ifc
lim e + lim e
fc^+œ 0 0
+œ
0
lim -1 + 1 = 1. i+œ e x dx = lim if e x dx =
fc^+œ fc^+œ fc^+œ eD fc^+œ 0
lim (-e—*) = = - lim (e—*) |fc = - lim e—+ lim e0 = - lim 1 + 1 = 1.
fc^+œ e
So,
the given property is an integral approximant.
_ + QQ (¿^
5. Calculate the characteristic integral .
e
0
EURASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES
Innovative Academy Research Support Center UIF = 8.3 | SJIF = 7.906 www.in-academy.uz
_ . . r+œ dx r0 dx r+œ dx r0 dx
Solving. I —- = I —- + I -- = lim I -- +
О j-œ 1+x2 J-œ 1+x2 J0 1+x2 a^-œ Ja 1+x2
-œ 1+x2
'0 1+x2
+ lim ff-^ = lim arctqx + lim arctqx I =
b^+œ 0 1+x2 a^-œ ]a b^+œ '0
= lim(arctq0 — arctqa) + lim (arctqb — arctqO) =
a^-œ b^+œ
lim arctqa + lim arctqb = —arctq(—œ) + arctq(+œ) =
г^—œ b^+œ
= —(—П)+П = П
V 2j 2 2
+
n
- = П.
2
Hence, the given eigenvalue is an integral approximant. 6. Calculate the characteristic integral f^ ^
Solving. At x = 0 the function under the integral has an infinite discontinuity. So, according to the definition
fidx = Um f1dx = Um inx I1 + lim (Inl — Ina) =
J0 x a^+0 a x a^+0 'a a^+0
= — lim Ina = —(—m) = m. Hence, the given integral is regressive.
2 d x
7. Calculate the characteristic integral f , -3,--
5 J-1 Mo—)2
Solving. Here the function under the integral [-1;2] has an infinite discontinuity at the internal point x=1 of the integration section. Therefore, by definition,
2 12 a
dx f dx f dx f dx
: + I - = lim I - +
J M(x — i)2 a^1 J '
Г dx Г
J 3J(x — l)2 = J
[ 3(x — l)2 _{ 3(x — l)2 J 3Kx-1)
d x
+lim I , — = lim I (
a^J 33(x — l)2 a^1 J
lim \(x — l) 3d(x — l) +
3(x — l)2 2
f —— 31 I a 3 I 12
+lim \(x — l) 3d(x — l) = 3 lim Vx — l I + 3lim Vx — l I =
a^1 J a^1 '-1 a^1 'a
a
= 3lim(Va—l + V2) + 3lim(VT — Va—l) = 3^V2 + 3 = a a =3(V2 + l).
a^1
References:
1. Соминский Сборник задач по вышей алебре, Москва, «Наука» 1972 Д.К. Фаддеев , И.С.
2. Abdalimov B. Oliy matematika. -Toshkent: O'qituychi,1994.
3. Abdalimov B. va boshqalar. Oliy matematikadan masalalar yechish bo'yicha qo'llanma. -Toshkent: O'qituvchi,1985.
4. Abdalimov B.,Salixov Sh. Oliy matematika qisqa kursi. -Toshkent: O'qituychi,1983.
5. Zaytsev I. Visshaya matematika. -Moskva: Visshaya shkola, 1991
6. Slobodskaya V.A. Kratkiy kurs visshey matematiki. -Moskva:Visshaya shkola,1969.
2
a
