Integration : Numerical Method


Numerical Integration (Trapezium Rule)



Curve တစ်ခု $y=f(x)$ ထောα€€်α€›ှိ $a$ α€”ှင့် $b$ α€€ြားပိုင်းထတွင်း area α€€ို $\displaystyle\int_{a}^{b} f(x) d x$ ဆိုတဲ့ integration α€”α€Š်းα€œα€™်းα€–ြင့် α€›ှာα€–ွေα€”ိုင်α€€ြောင်း α€€ို α€’ီ post α€™ှာ α€–ော်ပြခဲ့ပြီး α€–ြα€…်ပါα€α€š်။


α€žို့α€žော် integration method α€€ို α€žုံး၍ area α€›ှာα€–ွေα€›α€”် α€™α€œွα€š်α€€ူα€žော function α€™ျားα€…ွာα€›ှိပါα€α€š်။ α€’ီα€œိုထခြေထနေα€™ှာတော့ exact area α€€ို α€›ှာα€–ွေα€›α€”် α€™α€œွα€š်α€€ူတော့ပဲ Numerical method α€€ို α€žုံးပြီး estimated area α€€ိုα€žာ α€›ှာα€–ွေα€”ိုင်α€™ှာ α€–ြα€…်ပါα€α€š်။ ထဆိုပါ Estimated area α€€ိုα€›ှာα€–ွေα€žα€Š့်α€”α€Š်းα€€ို Trapezium Rule α€œို့ ခေါ်ပါα€α€š်။


Trapezium Rule α€€ို α€žုံးα€›α€”် ထောα€€်ပါ ထချα€€်α€™ျားα€€ို α€œုပ်ဆောင်α€›α€™ှာ α€–ြα€…်ပါα€α€š်။


  • ပေးထားα€žော curve ထောα€€်α€›ှိ area α€›ှာα€œိုα€žော α€€ြားပိုင်းထတွင်း ထကျα€š်တူα€Šီα€žော ထောင်α€œိုα€€် ထစိပ်ထပိုင်းα€™ျား ပိုင်းα€–ြတ်ရပါα€™α€š်။
  • Trapezium တစ်ခုခြင်းα€…ီ၏ area α€™ျားထားα€œုံးပေါင်းခြင်း α€–ြင့် α€œိုထပ်α€žော α€§α€›ိα€šာ၏ ခန့်α€™ှα€”်းတန်α€–ိုးα€€ို α€›α€›ှိα€™ှာ α€–ြα€…်ပါα€α€š်။
  • Trapezium ၏ area α€€ို ( $\dfrac{1}{2}$ × α€•ြိုင်α€žောထနားα€”ှα€…်ခုပေါင်းα€œα€’် × α€‘α€€ျα€š်) α€–ြင့်α€›ှာα€”ိုင်ပါα€α€š်။

ထောα€€်ပါပုံα€€ို ဆက်α€œα€€်α€œေ့α€œာα€€ြα€Š့်ပါ။ Curve $y=f(x)$ ထောα€€်α€›ှိ $a$ α€”ှင့် $b$ α€€ြားပိုင်းထတွင်း area α€€ို α€›ှာα€›α€”်ထတွα€€် ထကျα€š်တူα€Šီα€žော ထောင်α€œိုα€€် ထပိုင်း $6$ ပိုင်းပိုင်းထားပါα€α€š်။ ခန့်α€™ှα€”်း α€§α€›ိα€šာ = ထောင်α€œိုα€€် ထပိုင်း $6$ ပိုင်း α€§α€›ိα€šာα€™ျား ပေါင်းα€œα€’်α€Ÿု α€‘α€œွα€š်တကူ α€žိα€”ိုင်ပါα€α€š်။



α€€ြားပိုင်း $a$ α€”ှင့် $b$ α€€ြား ထပိုင်း $6$ ပိုင်း, ပိုင်းထားα€›α€™ှာ တစ်ပိုင်းα€…ီ၏ ထကျα€š်α€€ို $h$ α€Ÿုα€žα€်α€™ှတ်ပါα€™α€š်။ α€’ါ့α€€ြောင့် $h$ α€›ဲ့ α€‘α€œျားα€€ $a$ α€”ှင့် $b$ α€€ြား ထကွာထဝေးα€€ို $6$ α€”ှင့်α€…ားခြင်းα€–ြင့် α€›α€›ှိα€”ိုင်α€™ှာ α€–ြα€…်ပါα€α€š်။ $$ h = \dfrac{b-a}{6}$$ Trapezium တစ်ခုချင်းα€…ီ၏ α€§α€›ိα€šာα€€ို ထောα€€်ပါထတိုင်း α€›ှာα€šူα€”ိုင်ပါα€α€š်။


ပထမ Trapezium ၏ α€§α€›ိα€šာ $\dfrac{h}{2}\left(y_0+y_1\right)$
α€’ုတိα€š Trapezium ၏ α€§α€›ိα€šာ $\dfrac{h}{2}\left(y_1+y_2\right)$
တတိα€š Trapezium ၏ α€§α€›ိα€šာ $\dfrac{h}{2}\left(y_2+y_3\right)$
စတုတ္ထ Trapezium ၏ α€§α€›ိα€šာ $\dfrac{h}{2}\left(y_3+y_4\right)$
ပဉ္α€…α€™ Trapezium ၏ α€§α€›ိα€šာ $\dfrac{h}{2}\left(y_4+y_5\right)$
ဆဋ္α€Œα€™ Trapezium ၏ α€§α€›ိα€šာ $\dfrac{h}{2}\left(y_5+y_6\right)$

ထို့α€€ြောင့် α€–ော်ပြပါ ပုံ၏ α€žα€်α€™ှတ်ထားα€žော α€€ြားပိုင်းထတွင်းα€›ှိ α€§α€›ိα€šာα€€ို ထောα€€်ပါထတိုင်း α€›ှာα€šူα€”ိုင်ပါα€α€š်။


$\begin{aligned} \text{Estimated Area } &= \dfrac{h}{2}\left(y_0+y_1+y_1+y_2+y_2+y_3+y_3+y_4+y_4+y_5+y_5+y_6\right)\\\\ &= \dfrac{h}{2}\left[\left(y_0+y_6+ 2(y_1+y_2+y_3+y_4+y_5\right)\right]\\\\ \text{where } h = \dfrac{b-a}{6}& \end{aligned}$


α€‘α€€α€š်၍ α€žα€်α€™ှတ်α€€ြားပိုင်း $a$ α€”ှင့် $b$ ထတွင်း ထောင်α€œိုα€€် ထစိပ်ထပိုင်း $n$ ခု ပိုင်းα€–ြတ်ထားα€žα€Š် ဆိုပါα€…ို့။ ထိုထခါ ထပိုင်း တစ်ခု၏ ထကျα€š်α€™ှာ $$h = \dfrac{b-a}{n}$$ α€–ြα€…်α€™α€Š်။ ထိုထခါ ခန့်α€™ှα€”်း α€§α€›ိα€šာ၏ တန်းα€–ိုးα€™ှာ


$$A \approx \dfrac{h}{2}\left[\left(y_0+y_n+ 2(y_1+y_2+y_3+\ldots+y_{n-1}+y_n\right)\right]\\\\$$

ဆိုα€žော ပုံα€žေα€”α€Š်း α€–ြင့်α€›ှာα€šူα€”ိုင်α€™α€Š်။


Numerical Integration (Trapezium Rule)
$\displaystyle\int_{a}^{b} f(x) dx \approx \dfrac{h}{2}\left[y_{0}+y_{n}+2\left(y_{1}+y_{2}+y_{3}+\ldots+y_{n-1}\right)\right]$ where $h=\dfrac{b-a}{n}. $

α€‘α€œွα€š်α€™ှတ်α€›α€”်

$\displaystyle\int_{a}^{b} f(x)dx \approx$ half width of strip $\times$ (first $+$ last $+$ twice the sum of all the others).

α€‘α€€α€š်၍ Trapezium α€™ျားα€žα€Š် Curve ၏ ထပေါ်α€žို့ α€€ျော်α€”ေပါα€€ estimated area α€žα€Š် ထမှα€”်α€›ှိα€žော တန်α€–ိုးထက် α€€ြီးα€”ေပြီး ၎င်းထခြေထနေα€€ို over estimate α€Ÿုခေါ်α€žα€Š်။ Fig (1) α€€ို α€€ြα€Š့်ပါ။

Fig (1)

α€‘α€€α€š်၍ Trapezium α€™ျားα€žα€Š် Curve ၏ ထောα€€်တွင်α€žာ α€›ှိပါα€€ estimated area α€žα€Š် ထမှα€”်α€›ှိα€žော တန်α€–ိုးထောα€€် α€„α€š်α€”ေပြီး ၎င်းထခြေထနေα€€ို under estimate α€Ÿုခေါ်α€žα€Š်။ Fig (2) α€€ို α€€ြα€Š့်ပါ။

Fig (2)

ဆက်α€œα€€်၍ trapezium rule α€€ို α€‘α€žုံးပြုα€žα€Š့် numerical integration ဆိုင်α€›ာ ပုα€…္ဆာα€™ျားα€€ို တင်ပြပါα€™α€Š်။


Example (1)


The diagram shows the part of the curve $y=\dfrac{\ln x}{x}$ for $0 <x \leq 4$.
The curve cuts the $x$-axis at $A$ and its maximum point is $M$.
(i) Write down the co-ordinates of $A$.
(ii) Show that the $x$ co-ordinate of $M$ is e, and write down the $y$ co-ordinate of $M$ in terms of $e$.
(iii) Use the trapezium rule with three intervals to estimate the value of $$ \int_{1}^{4} \dfrac{\ln x}{x} dx $$ correct to $2$ decimal places.
(iv) State, with a reason, whether the trapezium rule gives an underestimate or an overestimate of the true value of the integral in part (iii).

Solution

$\begin{aligned} y=\dfrac{\ln x}{x}& \\\\ \text { When } y &=0, \\\\ \dfrac{\ln x}{x} &=0 \\\\ \text { Since } x & \neq 0, \\\\ \ln x &=0 \\\\ \therefore x &=1\\\\ \end{aligned}$
$\therefore$ Whe point $A$ is $(1,0)$.
$\begin{aligned} &\\ \dfrac{d y}{d x} &=\dfrac{x\left(\dfrac{1}{x}\right)-\ln x}{x^{2}} \\\\ &=\dfrac{1-\ln x}{x^{2}} \\\\ \dfrac{d y}{d x} &=0 \text { when } \\\\ \dfrac{1-\ln x}{x^{2}} &=0 \\\\ \ln x &=1 \\\\ x &=e \\\\ \therefore y &=\dfrac{\ln e}{e}=\dfrac{1}{e}\\\\ \end{aligned}$
$\therefore$ The point $M$ is $\left(e, \dfrac{1}{e}\right)\\\\ $.
To estimate $\displaystyle\int_{1}^{4} \dfrac{\ln x}{x} d x$ using trapezinm rule, with three intervals,
$\begin{aligned} &\\ a &=1, b=4, n=3 \\\\ \therefore h &=\dfrac{b-a}{n} \\\\ &=\dfrac{4-1}{3} \\\\ &=1 \end{aligned}$

$\begin{array}{|c||c|c|c|c|} \hline x & 1 & 2 & 3 & 4 \\ \hline y & 0 & 0.35 & 0.37 & 0.35 \\ \hline \end{array}$


$\begin{aligned} \therefore \displaystyle\int_{1}^{4} \dfrac{\ln x}{x} &\approx\dfrac{1}{2}[(0+0.35)+2(0.35+0.37)] \\\\ &=0.90\\\\ \end{aligned}$
It can be seen from the diagram that this is an under-estimate since the top edges of the strips all lie below the curve.


Example (2)

Use the trapezium rule with three intervals to find an approximation to $\displaystyle\int_{1}^{4} \frac{6}{1+\sqrt{x}}dx$. Give your answer correct to $5$ significant figures.

Solution

Let $y=\dfrac{6}{1+\sqrt{x}}$
For $\displaystyle\int_{1}^{4} y d x$,
$\begin{aligned} &\\ a &=1, b=4, n=3 \text { (given) } \\\\ \therefore \quad h &=\frac{b-a}{n} \\\\ &=\frac{4-1}{3} \\\\ &=1\\\\ \end{aligned}$

$\begin{array}{|c||c|c|c|c|} \hline x & 1 & 2 & 3 & 4 \\ \hline y & 3 & 2.4853 & 2.1962 & 2 \\ \hline \end{array}$

$\begin{aligned} & \displaystyle\int_{1}^{4} \frac{6}{1+\sqrt{x}} d x \\\\ \approx&\ \frac{1}{2}[(3+2)+2(2.4853+2.1962)] \\\\ =&\ 7.1815 \end{aligned}$

Example (3)

The diagram shows part of the curve $y=\dfrac{\mathrm{e}^{x}}{2 x}$. Use the trapezium rule with 4 intervals to estimate the area of the shaded region, giving your answer correct to 2 decimal places. State, with a reason, whether the trapezium rule gives an under-estimate or an over-estimate of the true value of the shaded area.

Solution

$\begin{aligned} y &=\frac{e^{x}}{2 x} \\\\ a &=1, b=3, n=4 \\\\ \therefore h &=\frac{b-a}{n} \\\\ &=\frac{3-1}{4} \\\\ &=\frac{1}{2}\\\\ \end{aligned}$

$\begin{array}{|c||c|c|c|c|c|} \hline x & 7 & 1.5 & 2 & 2.5 & 3 \\ \hline y & 1.359 & 1.494 & 1.847 & 1.437 & 3.348 \\ \hline \end{array}$


$\begin{aligned} & \displaystyle\int_{1}^{3} \frac{e^{x}}{2 x} d x \\\\ \approx&\ \frac{1}{4}[1.359+3.348+2(1.494+1.847+1.437)] \\\\ =&\ 3.566 \\\\ =&\ 3.57 \quad \text { (two decimal places) }\\\\ \end{aligned}$
It can be seen from the diagram that this is an under-estimate since the top edges of the strips all lie below the curve.

α€…ာဖတ်α€žူ၏ ထမြင်α€€ို α€œေးα€…ားα€…ွာα€…ောင့်α€™ျှော်α€œျα€€်!

Post a Comment

To be published, comments must be reviewed by the administrator *

Previous Post Next Post
πŸ’¬ 1
TM
Target Mathematics
Usually replies instantly
TM
Target Mathematics α€™ှ α€€ူα€Šီα€›α€”် α€‘α€žα€„့်α€›ှိပါα€α€š်။ α€˜ာα€™ျား α€žိα€›ှိချင်ပါα€žα€œဲ။ Target Mathematics Facebook Page α€™ှာα€œဲ တိုα€€်α€›ိုα€€် α€™ေးα€™ြα€”်းα€”ိုင်ပါα€α€š်