05.01.2013 Views

An Analytic Algorithm for Generalized Abel Integral Equation

An Analytic Algorithm for Generalized Abel Integral Equation

An Analytic Algorithm for Generalized Abel Integral Equation

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!