^{1}

^{2}

In this paper, we have studied and developed some analytical methods for solving some linear and fuzzy differential equations of the first degree with their fuzzy elementary conditions using the homotopy analytical method, using the approximate Padè method, and comparing the results with the exact solution of the equation, which indicates that the error decreases exponentially with the N degree of approximation. Also, we made some improvements to the solution results by using Padè enhancements as shown in the two examples with tables and figures.

In this thesis, we solve the equation linear through a systematic method called the homotopy trepidation method (HPM). The Chinese mathematician He first proposed this method in 1999 [

x ′ ( t ) = f ( t , x ( t ) ) , t 0 ≤ t ≤ t 0 + a (1)

Subject to the initial conditions

x ( t 0 ) = x 0 (2)

In this paper, we used the homotopy, analytical, blurry method with Padè Approximants to solve elementary values’ problems.

In this paper, show some basic definitions.

Definition of Fuzzy number 2.1 [

[ u ] r = { s ∈ R : u ( s ) ≥ r } (3)

And also

[ u ] 0 = { s ∈ R : u ( s ) > 0 } (4)

It is easy to prove that u is a fuzzy number if and only if [ u ] r is a convex compressed a subset of R for each r ∈ [ 0 , 1 ] and [ u ] ≠ ∅ and so on if u is a fuzzy number then

[ u ] r = [ u _ ( r ) , u ¯ ( r ) ] (5)

u _ ( r ) = min { s : s ∈ [ u ] r } (6)

u ¯ ( r ) = max { s : s ∈ [ u ] r } (7)

For each r ∈ [ 0 , 1 ] , [ u ] r is called an r-cut representation or a parameter form of a fuzzy number.

u : R → [ 0 , 1 ] defined by

u ( s ) = sup { r : u _ ( r ) ≤ s ≤ u ¯ ( r ) } (8)

Theorem 2.2 [

(i) [ u + v ] r = [ u ] r + [ v ] r = [ u _ r + v _ r , u ¯ r + v ¯ r ] (9)

(ii) [ λ u ] r = λ [ u ] r = { [ λ u ¯ r , λ u _ r ] : λ ≥ 0 [ λ u ¯ r , λ u _ r ] : λ < 0 (10)

(iii) [ u v ] r = [ u ] r [ v ] r = [ min S r , max S r ] , where

S r = { u _ r v _ r , u _ r v ¯ r , u ¯ r v _ r , u ¯ r v ¯ r } (11)

All the fuzzy numbers that are collected by adding and multiplying are defined by part (i) and (ii), is a convex cone [

Definition 2.3 [

1) If μ A ( X ) = 0 if it is not an element in A.

2) If μ A ( X ) ≃ 1 if the degree of belonging is high of A.

3) If μ A ( X ) ≃ 0 if the degree of belonging is weak.

4) If μ A ( X ) = 1 if zero is an element in A.

μ A ( X ) : E → [ 0 , 1 ]

μ A ( X ) = { 1 x ∈ A 0 x ∉ A (12)

Definition 2.4 [

R N , M ( x ) = P A [ N / M ] = P N ( x ) Q M ( x ) , a ≤ x ≤ b (13)

Q M ( x ) , P N ( x ) Polynomials are solved so that the function f ( x ) and R N , M ( x ) plus their derivatives to the N + M rank are identical when x = 0, Assuming that f ( x ) is analytical and has a (Maclaurin Series), we get:

f ( x ) = a 0 + a 1 x + a 2 x 2 + ⋯ + a k x k + ⋯ (14)

I will discuss some of basic concepts of the homotopy method of analysis.

If we have the following non-linear equation

N [ y ( t ) ] = 0 , t ≥ 0 (15)

Since N represents a non-linear operator, y is an unknown function, and t is an independent variable. The zero-order distortion equation was derived by generalizing and formulating the homotopy method:

( 1 − q ) L [ ∅ ( t , q ) − y 0 ( t ) ] = q h ¯ H ( t ) N [ ∅ ( t , q ) ] (16)

Since q ∈ [ 0 , 1 ] the inclusion parameter is called the homotopy parameter, L is the auxiliary linear bound that achieves the property of L ( 0 ) = 0 when L ( y ) = 0 , ∅ ( t , q ) solve Equation (15), y 0 ( t ) initial estimation of the exact solution y ( t ) , h ¯ ≠ 0 is called the convergence control parameter, H ( t ) ≠ 0 auxiliary function, respectively, within the homotopy we can have great freedom to choose the auxiliary linear term L, the initial estimate y 0 ( t ) , And the parameters h ¯ , H ( t ) when q = 0 , due to the property L ( 0 ) = 0 , we get from Equation (16) the solution:

∅ ( t , 0 ) = y 0 ( t ) (17)

When q = 1 , and since H ( t ) ≠ 0 , h ¯ ≠ 0 , and Equation (15) is equivalent to Equation (16) then:

∅ ( t , 1 ) = y ( t ) (18)

As y ( t ) is a solution to the original Equation (15) and so whenever the homotopy parameter q increases from zero to one, the solution ∅ ( t , q ) changes continuously from the initial guess y 0 ( t ) to the exact solution y ( t ) , this A type of constant change in solution is called homotypic deformation. By expanding ∅ ( t , q ) using the Taylor series with respect to q, we get the following:

∅ ( t , q ) = y 0 ( t ) + ∑ m = 1 + ∞ y m ( t ) q m (19)

y m ( t ) = 1 m ! ∂ m ∅ ( t , q ) ∂ q m | q = 0 (20)

If h, H ( t ) , y 0 ( t ) , L are correctly chosen then Equation (19) approaches when q = 1, therefore based on this hypothesis, the solution becomes as follows:

y ( t ) = ∅ ( t , 1 ) = y 0 ( t ) + ∑ m = 1 + ∞ y m ( t ) (21)

Now we know the vector y n ( t ) as follows:

y n ( t ) = { y 0 ( t ) , y 1 ( t ) , y 2 ( t ) , ⋯ , y n ( t ) } (22)

By deriving the zero deformation Equation (16), m times for parameter q, then substituting q for zero, and finally, we divide the equation by m! (the deformation equations) of the order m are obtained as follows:

L [ y m ( t ) − x m y m − 1 ( t ) ] = h ¯ H ( t ) R m − 1 ( y m − 1 ( t ) ) (23)

According to the conditions

y m ( 0 ) = 0 (24)

As that

R m − 1 ( y m − 1 ( t ) ) = 1 ( m − 1 ) ! ∂ m − 1 N [ ∅ ( t , q ) ] ∂ q m − 1 | q = 0 (25)

x m = { 0 m ≤ 1 1 m > 1 (26)

We will demonstrate the fuzzy homotopy analytic method with Padè approximants to solve elementary value problems:

Case I. If x ( t ) is differentiable,

[ D 1 x ( t ) ] r = [ x _ ′ r ( t ) , x ¯ ′ r ( t ) ] (27)

To solve (FIVP) we follow the following steps:

Step (i): to solve the system of (ODEs) subsequent for x _ r ( t ) , x ¯ r ( t ) :

x _ ′ r ( t ) = f 1 , r ( t , x _ r ( t ) , x ¯ r ( t ) ) x ¯ ′ r ( t ) = f 2 , r ( t , x _ r ( t ) , x ¯ r ( t ) ) (28)

Subject to the initial conditions

x _ r ( t 0 ) = x _ r 0 x ¯ r ( t 0 ) = x ¯ r 0 (29)

Step (ii): ensure that the solution [ x _ r ( t ) , x ¯ r ( t ) ] and its derivative [ x _ ′ r ( t ) , x ¯ ′ r ( t ) ] are functional sets for each r ∈ [ 0 , 1 ] ,

Step (iii): using Equation (8) to construct a solution x ( t ) such that

[ x ( t ) ] r = [ x _ r ( t ) , x ¯ r ( t ) ] , r ∈ [ 0 , 1 ] (30)

Case II. If x ( t ) is differentiable, then

[ D 2 x ( t ) ] r = [ x ¯ ′ r ( t ) , x _ ′ r ( t ) ]

And to solve (FIVP) we follow the following steps:

Step (i): to solve the system of (ODEs) subsequent for x _ r ( t ) , x ¯ r ( t ) :

x _ ′ r ( t ) = f 2 , r ( t , x _ r ( t ) , x ¯ r ( t ) ) x ¯ ′ r ( t ) = f 1 , r ( t , x _ r ( t ) , x ¯ r ( t ) ) (31)

Subject to the initial conditions

x _ r ( t 0 ) = x _ r 0 ( t ) x ¯ r ( t 0 ) = x ¯ r 0 ( t ) (32)

Step (ii): ensure that the solution [ x _ r ( t ) , x ¯ r ( t ) ] and its derivative [ x ¯ ′ r ( t ) , x _ ′ r ( t ) ] are functional sets for each r ∈ [ 0 , 1 ] ,

Step (iii): using Equation (8) to construct a solution x ( t ) such that

[ x ( t ) ] r = [ x _ r ( t ) , x ¯ r ( t ) ] , r ∈ [ 0 , 1 ] (33)

Case III. Gather approximate iterations x 0 ( t ) , x 1 ( t ) , ⋯ , x n ( t )

ψ x _ r , k ( t ) = ∑ m = 0 k − 1 x _ r , m ( t ) and ψ x ¯ r , k ( t ) = ∑ m = 0 k − 1 x ¯ r , m ( t ) (34)

its value determined by the number of limits after summing the approximate iterations k, representing the number of iterations, when collecting the frequencies we substitute for the best values of h ¯ , h ¯ = − 1 so that the sequence (34) is convergent we get the sequenced solution to the problem of elementary fuzzy values in an analytical homotopy method.

Case IIII. Using the Padè approximation [N/M] where M , N ∈ ℕ ∪ { 0 } we obtain:

P A [ N / M ] = P N ( x ) Q M ( x ) = ∑ v = 0 N b v x v 1 + ∑ z = 1 M c z x z (35)

By taking different values of t and substituting them into the Padè sequence, (35), we obtain approximate solutions in an analytical homotopy method with the upper and lower Padè approximations.

Example 4.1. Consider the following linear (FIVP) [

x ′ r ( t ) = 2 t x ( t ) + t u , 0 ≤ t ≤ 1 (36)

Subject to the fuzzy initial condition

x ( 0 ) = u , where u = [ r − 1 , 1 − r ] , x ( 0 ) = [ r − 1 , 1 − r ] (37)

The fuzzy (1)―differentiable exact solution is x _ ( t ) = 1 2 ( 3 e t 2 − 1 ) ( r − 1 ) .

When the fuzzy (2), differentiable exact of the solution is x ¯ ( t ) = 1 2 ( 3 e t 2 − 1 ) ( 1 − r ) .

Solution:

The first Case 1: When the system of the (ODEs) matching to (1) and (2), then differentiability is

The first step: If x(t) is differentiable (1) then

x _ ′ r ( t ) = 2 t x _ r ( t ) + t ( r − 1 ) x ¯ ′ r ( t ) = 2 t x ¯ r ( t ) + t ( 1 − r ) (38)

Subject to the initial condition

x _ r ( 0 ) = ( r − 1 ) x ¯ r ( 0 ) = ( 1 − r ) (39)

The second step:

To find the iterations x _ r , 1 ( t ) , x ¯ r , 1 ( t ) , ⋯ lower and upper we will use the equation to find the homotopic equation in the formula:

x _ r , m ( t ) = X m x _ r , m − 1 ( t ) + h ¯ ∫ t 0 t R x _ r , m ( x _ r , m − 1 ( τ ) , x ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ x ¯ r , m ( t ) = X m x ¯ r , m − 1 ( t ) + h ¯ ∫ t 0 t R x ¯ r , m ( x _ r , m − 1 ( τ ) , x ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ (40)

Through the equation for the deformation of order m, (40) m ≥ 1, we get:

x _ r , m ( t ) = X m x _ r , m − 1 ( t ) + h ¯ ∫ t 0 t x _ ′ r , m − 1 ( t ) − 2 t x _ r , m − 1 ( t ) − t ( r − 1 ) ( 1 − x m ) d τ x ¯ r , m ( t ) = X m x ¯ r , m − 1 ( t ) + h ¯ ∫ t 0 t x ¯ ′ r , m − 1 ( t ) − 2 t x ¯ r , m − 1 ( t ) − t ( 1 − r ) ( 1 − x m ) d τ (41)

When m = 1, 2, 3, 4 by substituting into Equation (41) we get:

x _ r , 1 ( t ) = − 3 2 h ¯ t 2 r + 3 2 h ¯ t 2

x ¯ r , 1 ( t ) = − 3 2 h ¯ t 2 + 3 2 h ¯ t 2 r

x _ r , 2 ( t ) = − 3 2 h ¯ t 2 r + 3 2 h ¯ t 2 + 3 4 h ¯ 2 t 4 r − 3 4 h ¯ 2 t 4 − 3 2 h ¯ 2 t 2 r + 3 2 h ¯ 2 t 2

x ¯ r , 2 ( t ) = − 3 2 h ¯ t 2 + 3 2 h ¯ t 2 r + 3 4 h ¯ 2 t 4 − 3 4 h ¯ 2 t 4 r − 3 2 h ¯ 2 t 2 + 3 2 h ¯ 2 t 2 r

The third case:

Sum the approximate iterations lower and upper x _ r , 0 ( t ) , x ¯ r , 0 ( t ) , x _ r , 1 ( t ) , x ¯ r , 1 ( t ) , ⋯ .

After substituting in the best values of, h ¯ , h ¯ = − 1 , r = 0.9 , we get an approximate solution as follows:

ψ x _ r , k ( t ) = ∑ m = 0 k − 1 x _ r , m ( t ) , ψ x ¯ r , k ( t ) = ∑ m = 0 k − 1 x ¯ r , m ( t )

ψ x _ r , 5 ( t ) = − 0.00625 t 8 − 0.02500 t 6 − 0.07500 t 4 − 0.15000 t 2 − 0.1 ψ x ¯ r , 5 ( t ) = 0.00625 t 8 + 0.02500 t 6 + 0.07500 t 4 + 0.15000 t 2 + 0.1 (42)

Fourth case: We connect the solution series (42) to Padè approximation P A [ N / M ] and we get:

P A [ N / M ] = P N ( x ) Q M ( x ) = ∑ v = 0 N b v x v 1 + ∑ z = 1 M c z x z (43)

lower : P A [ 4 / 4 ] = − 0.09999 − 0.09999 t 2 − 0.00833 t 4 0.99999 − 0.49999 t 2 + 0.08333 t 4 upper : P A [ 4 / 4 ] = 0.09999 + 0.09999 t 2 + 0.00833 t 4 0.99999 − 0.49999 t 2 + 0.08333 t 4 (44)

By taking different values of (t) and substituting them into the Padè series (44), we obtain approximate solutions by homotopic analytical method with Padè approximation, as shown in

The second case II:

The first step:

For find the upper and lower values of x _ r , m ( t ) , x ¯ r , m ( t ) , if x ( t ) is differentiable then:

x _ ′ r ( t ) = 2 t x ¯ r ( t ) + t ( r − 1 ) x ¯ ′ r ( t ) = 2 t x _ r ( t ) + t ( 1 − r ) (45)

subject to the initial condition

x _ r ( 0 ) = ( r − 1 ) x ¯ r ( 0 ) = ( 1 − r ) (46)

The second step:

To find the iterations x r , 1 ( t ) , x r , 2 ( t ) , x r , 3 ( t ) , ⋯ upper and lower by using the equations

x _ r , m ( t ) = X m x ¯ r , m − 1 ( t ) + h ¯ ∫ t 0 t R x ¯ r , m ( x _ r , m − 1 ( τ ) , x ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ x ¯ r , m ( t ) = X m x _ r , m − 1 ( t ) + h ¯ ∫ t 0 t R x _ r , m ( x _ r , m − 1 ( τ ) , x ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ (47)

When m = 1, by substituting in the upper and lower Equations (47) we get:

x _ r , 1 ( t ) = − 3 2 h ¯ t 2 r + 3 2 h ¯ t 2

x ¯ r , 1 ( t ) = − 3 2 h ¯ t 2 + 3 2 h ¯ t 2 r

When m = 2 , 3 , ⋯ in Equation (54) and using the Riemann-Liouville integral properties, we get:

x _ r , 2 ( t ) = − 3 2 h ¯ t 2 r + 3 2 h ¯ t 2 + 3 4 h ¯ 2 t 4 r − 3 4 h ¯ 2 t 4 − 3 2 h ¯ 2 t 2 r + 3 2 h ¯ 2 t 2

x ¯ r , 2 ( t ) = − 3 2 h ¯ t 2 + 3 2 h ¯ t 2 r + 3 4 h ¯ 2 t 4 − 3 4 h ¯ 2 t 4 r − 3 2 h ¯ 2 t 2 + 3 2 h ¯ 2 t 2 r

The third case III:

We assume the upper and lower approximate iterations x r , 0 ( t ) , x r , 1 ( t ) x r , 2 ( t ) , x r , 3 ( t ) , x r , 4 ( t ) .

After substituting in the best values of h ¯ , h ¯ = − 1 , we get an approximate solution of the upper and lower values

ψ x _ r , 5 ( t ) = − 0.00625 t 8 − 0.02500 t 6 − 0.07500 t 4 − 0.15000 t 2 − 0.1 ψ x ¯ r , 5 ( t ) = 0.00625 t 8 + 0.02500 t 6 + 0.07500 t 4 + 0.15000 t 2 + 0.1 (48)

The fourth step:

We connect the solution series (48) to the approximations of the upper and lower pads P A [ N / M ] and we get:

P A [ N / M ] = P N ( x ) Q M ( x ) = ∑ v = 0 N b v x v 1 + ∑ z = 1 M c z x z (49)

By taking different values of (t) and substituting them into the Padè series (49), we get approximate solutions in an analytical homotopic method with upper and lower Padè approximations as shown in

Example 4.2. Consider the following linear (FIVP) [

y ′ r ( t ) = y ( t ) , 0 ≤ t ≤ 1 (50)

Subject to the fuzzy initial condition

y _ ( 0 ) = ( r − 1 ) , y ¯ ( 0 ) = ( 1 − r ) (51)

The fuzzy (1)―differentiable exact solution is y _ ( t ) = ( r − 1 ) e t .

When the fuzzy (2)―differentiable exact of the solution is y ¯ ( t ) = ( 1 − r ) e t .

Solution: The first Case 1: When the system of the (ODEs) equivalent to (1) and (2), then differentiability is

The first step: If y ( t ) is differentiable (1) and (2) then:

Absolute Error | Approximation PAD x _ | Absolute Error | Approximation HAM x _ | Exact Solution | t |
---|---|---|---|---|---|

0 | −0.10 | 0 | −0.10 | −0.10 | 0 |

2.10428 × 10^{−14} | −0.10150 | 1.25208 × 10^{−13} | −0.10150 | −0.10150 | 0.1 |

2.22060 × 10^{−11} | −0.10612 | 1.28858 × 10^{−10} | −0.10612 | −0.10612 | 0.2 |

1.34668 × 10^{−9} | −0.11412 | 7.49328 × 10^{−9} | −0.11412 | −0.11412 | 0.3 |

2.56747 × 10^{−8} | −0.12602 | 1.34648 × 10^{−7} | −0.12602 | −0.12602 | 0.4 |

2.62207 × 10^{−7} | −0.14260 | 1.27344 × 10^{−6} | −0.14260 | −0.14260 | 0.5 |

1.81950 × 10^{−6} | −0.16499 | 8.03618 × 10^{−6} | −0.16499 | −0.16499 | 0.6 |

9.74307 × 10^{−6} | −0.19483 | 3.84079 × 10^{−5} | −0.19480 | −0.16499 | 0.7 |

4.34537 × 10^{−5} | −0.23442 | 1.49955 × 10^{−4} | −0.23432 | −0.23447 | 0.8 |

1.69612 × 10^{−4} | −0.28701 | 5.02252 × 10^{−4} | −0.28668 | −0.28718 | 0.9 |

5.99417 × 10^{−4} | −0.35714 | 1.49227 × 10^{−3} | −0.35625 | −0.35774 | 1.0 |

Absolute Error | Approximation PAD x ¯ | Absolute Error | Approximation HAM x ¯ | Exact Solution | t |
---|---|---|---|---|---|

0 | 0.10 | 0 | 0.1 | 0.10 | 0 |

2.10428 × 10^{−14} | 0.10150 | 1.25208 × 10^{−13} | 0.10150 | 0.10150 | 0.1 |

2.22060 × 10^{−11} | 0.10612 | 1.28858 × 10^{−10} | 0.10612 | 0.10612 | 0.2 |

1.34668 × 10^{−9} | 0.11412 | 7.49328 × 10^{−6} | 0.11412 | 0.11412 | 0.3 |

2.56747 × 10^{−8} | 0.12602 | 1.34648 × 10^{−7} | 0.12602 | 0.12602 | 0.4 |

2.62207 × 10^{−7} | 0.14260 | 1.27344 × 10^{−6} | 0.14260 | 0.14260 | 0.5 |

1.81950 × 10^{−6} | 0.16499 | 8.03618 × 10^{−6} | 0.16499 | 0.16499 | 0.6 |

9.74307 × 10^{−6} | 0.19483 | 3.84079 × 10^{−5} | 0.19480 | 0.16499 | 0.7 |

4.34537 × 10^{−5} | 0.23442 | 1.49955 × 10^{−4} | 0.23432 | 0.23447 | 0.8 |

1.69612 × 10^{−4} | 0.28701 | 5.02252 × 10^{−4} | 0.28668 | 0.28718 | 0.9 |

5.99417 × 10^{−4} | 0.35714 | 1.49227 × 10^{−3} | 0.35625 | 0.35774 | 1.0 |

y _ ′ r ( t ) = y _ r ( t ) y ¯ ′ r ( t ) = y ¯ r ( t ) (52)

subject to the initial condition

y _ r ( 0 ) = ( r − 1 ) y ¯ r ( 0 ) = ( 1 − r ) (53)

The second step: To find the iterations y _ r , 1 ( t ) , y ¯ r , 1 ( t ) , ⋯ lower and upper we will use the equation to find the homotopic equation:

y _ r , m ( t ) = X m y _ r , m − 1 ( t ) + h ¯ ∫ t 0 t R y _ r , m ( y _ r , m − 1 ( τ ) , y ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ y ¯ r , m ( t ) = x m y ¯ r , m − 1 ( t ) + h ¯ ∫ t 0 t R y ¯ r , m ( y _ r , m − 1 ( τ ) , y ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ (54)

When m = 1, 2, 3, 4 by substituting into the Equation (54), we get:

y _ r , 1 ( t ) = h ¯ t − h ¯ r t , y ¯ r , 1 ( t ) = − h ¯ t + h ¯ r t

y _ r , 2 ( t ) = h ¯ t − h ¯ r t − 1 2 h ¯ 2 t 2 + 1 2 h ¯ 2 t 2 r + h ¯ 2 t − h ¯ 2 r t ,

y ¯ r , 2 ( t ) = − h ¯ t + h ¯ r t + 1 2 h ¯ 2 t 2 − 1 2 h ¯ 2 t 2 r − h ¯ 2 t + h ¯ 2 r t

The third case: Sum the approximate iterations lower and upper y r , 0 ( t ) , y r , 1 ( t ) , ⋯ .

After substituting in the best values of h ¯ , h ¯ = − 1 , r = 0.9 , we get an approximate solution as follows:

ψ y _ r , k ( t ) = ∑ m = 0 k − 1 y _ r , m , ψ y ¯ r , k ( t ) = ∑ m = 0 k − 1 y ¯ r , m

ψ y _ r , 5 ( t ) = − 0.00416 t 4 − 0.01666 t 3 − 0.05000 t 2 − 0.1 t − 0.1 ψ y ¯ r , 5 ( t ) = 0.00416 t 4 + 0.01666 t 3 + 0.05000 t 2 + 0.1 t + 0.1 (55)

Fourth case: We connect the solution series (55) to Padè approximation P A [ N / M ] and we get:

P A [ N / M ] = P N ( x ) Q M ( x ) = ∑ v = 0 N b v x v 1 + ∑ z = 1 M c z x z (56)

lower : P A [ 1 / 3 ] = − 0.1 − 0.25000 t 1.0 − 0.74999 t + 0.24999 t 2 − 0.04166 t 3 upper : P A [ 1 / 3 ] = 0.1 + 0.02500 t 1.0 − 0.74999 t + 0.24999 t 2 − 0.04166 t 3 (57)

By taking different values of (t) and substituting them into the Padè series (57), we obtain approximate solutions by homotopic analytical method with Padè approximation, as shown in

The second case II: The first step:

For find the upper and lower values of x ¯ r , m ( t ) , x _ r , m ( t ) , if x ( t ) is differentiable then:

y _ ′ r ( t ) = y _ r ( t ) y ¯ ′ r ( t ) = y ¯ r ( t ) (58)

Subject to the initial condition

y _ r ( 0 ) = ( r − 1 ) y ¯ r ( 0 ) = ( 1 − r ) (59)

The second step: To find the iterations y r , 1 ( t ) , y r , 2 ( t ) , ⋯ upper and lower by using the equations

y _ r , m ( t ) = X m y ¯ r , m − 1 ( t ) + h ¯ ∫ t 0 t R y ¯ r , m ( y _ r , m − 1 ( τ ) , y ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ y ¯ r , m ( t ) = X m y _ r , m − 1 ( t ) + h ¯ ∫ t 0 t R y _ r , m ( y _ r , m − 1 ( τ ) , y ¯ r , m − 1 ( τ ) ) ( 1 − x m ) d τ (60)

The third case: We assume the upper and lower approximate iterations y r , 0 ( t ) , y r , 1 ( t ) , ⋯ .

After substituting in the best values of h ¯ , h ¯ = − 1 , we get an approximate solution of the upper and lower values

ψ y _ r , 5 ( t ) = − 0.00416 t 4 − 0.01666 t 3 − 0.05000 t 2 − 0.1 t − 0.1 ψ y ¯ r , 5 ( t ) = 0.00416 t 4 + 0.01666 t 3 + 0.05000 t 2 + 0.1 t + 0.1 (61)

The fourth step:

We connect the solution series (61) to the approximations of the upper and lower pads P A [ N / M ] and we get:

P A [ N / M ] = P N ( x ) Q M ( x ) = ∑ v = 0 N b v x v 1 + ∑ z = 1 M c z x z (62)

By taking different values of (t) and substituting them into the Padè series (62), we get approximate solutions in an analytical homotopy method with upper and lower Padè approximations as shown in

Absolute Error | Approximation PAD x _ | Absolute Error | Approximation HAM x _ | Exact Solution | t |
---|---|---|---|---|---|

0 | −0.1 | 0 | −0.1 | −0.1 | 0 |

2.40151 × 10^{−9} | −0.11051 | 8.47423 × 10^{−9} | 0.11051 | 0.11051 | 0.1 |

8.86663 × 10^{−8} | −0.12214 | 2.75816 × 10^{−7} | −0.12214 | −0.12214 | 0.2 |

7.77537 × 10^{−7} | −0.13498 | 2.13075 × 10^{−6} | −0.13498 | 0.13498 | 0.3 |

3.78701 × 10^{−6} | −0.14918 | 9.13643 × 10^{−6} | −0.14917 | −0.14918 | 0.4 |

1.33691 × 10^{−5} | −0.16488 | 2.83770 × 10^{−5} | −0.16484 | −0.16487 | 0.5 |

3.85161 × 10^{−5} | −0.18225 | 7.18800 × 10^{−5} | −0.18214 | −0.18221 | 0.6 |

9.64731 × 10^{−5} | −0.20147 | 1.58187 × 10^{−4} | −0.20121 | −0.20137 | 0.7 |

2.18184 × 10^{−4} | −0.22277 | 3.14092 × 10^{−4} | −0.22224 | −0.22255 | 0.8 |

4.56586 × 10^{−4} | −0.24641 | 5.76561 × 10^{−4} | −0.24538 | −0.24596 | 0.9 |

8.99089 × 10^{−4} | −0.27272 | 9.94849 × 10^{−4} | −0.27083 | −0.27182 | 1.0 |

Absolute Error | Approximation PAD x ¯ | Absolute Error | Approximation HAM x ¯ | Exact Solution | t |
---|---|---|---|---|---|

0 | 0.1 | 0 | 0.1 | 0.1 | 0 |

2.40151 × 10^{−9} | 0.11051 | 8.47423 × 10^{−9} | 0.11051 | 0.11051 | 0.1 |

8.86663 × 10^{−8} | 0.12214 | 2.75816 × 10^{−7} | 0.12214 | 0.12214 | 0.2 |

7.77537 × 10^{−7} | 0.13498 | 2.13075 × 10^{−6} | 0.13498 | 0.13498 | 0.3 |

3.78701 × 10^{−6} | 0.14918 | 9.13643 × 10^{−3} | 0.14917 | 0.14918 | 0.4 |

1.33691 × 10^{−5} | 0.16488 | 2.83770 × 10^{−5} | 0.16484 | 0.16487 | 0.5 |

3.85161 × 10^{−5} | 0.18225 | 7.18800 × 10^{−5} | 0.18214 | 0.18221 | 0.6 |

9.64731 × 10^{−5} | 0.20147 | 1.58187 × 10^{−4} | 0.20121 | 0.20137 | 0.7 |

2.18184 × 10^{−4} | 0.22277 | 3.14092 × 10^{−4} | 0.22224 | 0.22255 | 0.8 |

4.56586 × 10^{−4} | 0.24641 | 5.76561 × 10^{−4} | 0.24538 | 0.24596 | 0.9 |

8.99089 × 10^{−4} | 0.27272 | 9.94849 × 10^{−4} | 0.27083 | 0.27182 | 1.0 |

In the case of fuzzy differential equations, it is found that the results of numerical solutions of the homotopy method with Padè approximations are very close when compared with the exact solutions (accurate). There is no problem with this sentence and it has been validated. The proposed method is very powerful, straightforward, and promising algorithm in finding a solution to the analytical approximate form of broad categories of Fuzzy Initial Values Problems (FIVPS).

For the University of Al-Mosul College of Education for Pure Sciences Department of Mathematics for the assistance they provided to improve the quality of this work.

The authors declare no conflicts of interest regarding the publication of this paper.

Ali1, R.H. and Ibraheem, K.I. (2020) Solution of Fuzzy Initial Value Problems Using Homotopy Analysis Method and Padè Approximate. Open Access Library Journal, 7: e6803. https://doi.org/10.4236/oalib.1106803