An Analytic Algorithm for Generalized Abel Integral Equation
An Analytic Algorithm for Generalized Abel Integral Equation
An Analytic Algorithm for Generalized Abel Integral Equation
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>An</strong>alytic algorithm <strong>for</strong> generalized <strong>Abel</strong> integral equation 229<br />
y ( x)<br />
= x<br />
2<br />
+<br />
27<br />
40<br />
x<br />
8<br />
3<br />
−<br />
x<br />
∫<br />
0<br />
y ( t )<br />
( x − t )<br />
2<br />
with the exact Solution y ( x)<br />
= x .<br />
1<br />
3<br />
dt ,<br />
Case 1(a) Homotopy perturbation method<br />
27 8<br />
2<br />
3<br />
Taking L 0 ( x)<br />
= x + x , the various iterates are obtained from equation ( 12 ).<br />
40<br />
The first few iterates are as follows:<br />
L ( x ) =<br />
1<br />
−<br />
27<br />
40<br />
x<br />
8<br />
3<br />
4 ⎛ ⎛ 2 ⎞ ⎞<br />
x ⎜ Γ ⎜ ⎟ ⎟<br />
3<br />
2 ( )<br />
⎝ ⎝ ⎠<br />
L x =<br />
⎠<br />
12<br />
x<br />
L ( x)<br />
= −<br />
3<br />
4<br />
−<br />
3<br />
⎛ ⎛ 2 ⎞ ⎞<br />
⎜ Γ⎜<br />
⎟ ⎟<br />
⎝ ⎝ 3 ⎠ ⎠<br />
12<br />
2 x<br />
+<br />
3<br />
10<br />
27<br />
⎛ 3 ⎛<br />
⎜ Γ ⎜<br />
⎝ ⎝<br />
⎛ 13<br />
Γ ⎜<br />
⎝ 3<br />
2<br />
3<br />
⎞<br />
⎟<br />
⎠<br />
⎞<br />
⎟<br />
⎠<br />
⎞<br />
⎟<br />
⎠<br />
10 2 11<br />
3 ⎛ ⎞ ⎛ ⎞<br />
x Γ ⎜ ⎟ Γ ⎜ ⎟<br />
⎝ 3 ⎠ ⎝ 3 ⎠<br />
,<br />
⎛ 13 ⎞<br />
40 Γ ⎜ ⎟<br />
⎝ 3 ⎠<br />
14 ⎛ 2<br />
3 ⎛ ⎞ ⎞<br />
243 x ⎜ Γ⎜<br />
⎟ ⎟<br />
⎝ ⎝ 3 ⎠<br />
−<br />
⎠<br />
6160<br />
The Fig.1 shows the absolute error between the exact solution y ( x )<br />
approximate solution y a ( x ) obtained from (11) by truncating it at level n=13.<br />
4. 10� 7<br />
3. 10� 7<br />
2. 10� 7<br />
1. 10� 7<br />
2<br />
,<br />
3<br />
, ...<br />
0.2 0.4 0.6 0.8 1.0<br />
Figure 1. The absolute error <strong>for</strong> Example1, case 1(a) (n=13).<br />
Case 1(b) Now choosing a different initial guess<br />
L ( x)<br />
= x,<br />
the following iterates of the solution are obtained<br />
0<br />
9 5<br />
2<br />
3<br />
L 1 ( x ) =<br />
− x + x − x +<br />
10<br />
27<br />
40<br />
x<br />
8<br />
3<br />
,<br />
(20)<br />
and the