In this paper, we introduce and discuss an algorithm for the numerical solution of two- dimensional fractional dispersion equation.  The algorithm for the numerical solution of this equation is based on explicit finite difference approximation. Consistency, conditional stability, and convergence of this numerical method are described. Finally, numerical example is presented to show the dispersion behavior according to the order of the fractional derivative and we demonstrate that our explicit finite difference approximation is a computationally efficient method for solving two-dimensional fractional dispersion equation

This paper presents a numerical solution to the inverse problem consisting of recovering time-dependent thermal conductivity and  heat source coefficients  in the one-dimensional  parabolic heat equation.   This  mathematical  formulation  ensures that the inverse problem  has a unique  solution.   However, the problem  is still  ill-posed since small errors  in the input data lead to a drastic  amount  of errors in the output coefficients.  The  finite  difference method  with  the Crank-Nicolson  scheme is adopted  as a direct  solver of the problem in a fixed domain.   The inverse problem is solved sub

 In this paper, we study the convergence theorems of the Modified Ishikawa iterative sequence with mixed errors for the uniformly continuous mappings and solving nonlinear uniformly continuous mappings equation in arbitrary real Banach space.

     The variational iteration method is used to deal with linear and nonlinear differential equations. The main characteristics of the method lie in its flexibility and ability to accurately and easily solve nonlinear equations. In this work, a general framework is presented for a variational iteration method for the analytical treatment of partial differential equations in fluid mechanics. The Caputo sense is used to describe fractional derivatives. The time-fractional Kaup-Kupershmidt (KK) equation is investigated, as it is the solution of the system of partial differential equations via the Boussinesq-Burger equation. By comparing the results that are obtained by the variational iteration method with those obtained by the two-dim

The aim of this article is to solve the Volterra-Fredholm integro-differential equations of fractional order numerically by using the shifted Jacobi polynomial collocation method. The Jacobi polynomial and collocation method properties are presented. This technique is used to convert the problem into the solution of linear algebraic equations. The fractional derivatives are considered in the Caputo sense. Numerical examples are given to show the accuracy and reliability of the proposed technique.

In this article, the backstepping control scheme is proposed to stabilize the fractional order Riccati matrix differential equation with retarded arguments in which the fractional derivative is presented using Caputo's definition of fractional derivative. The results are established using Mittag-Leffler stability. The fractional Lyapunov function is defined at each stage and the negativity of an overall fractional Lyapunov function is ensured by the proper selection of the control law. Numerical simulation has been used to demonstrate the effectiveness of the proposed control scheme for stabilizing such type of Riccati matrix differential equations.

     Fractional calculus has paid much attention in recent years, because it plays an essential role in many fields of science and  engineering, where the study of stability theory of fractional differential equations emerges to be very important. In this paper, the stability of fractional order ordinary differential equations will be studied and introduced the backstepping method. The Lyapunov function  is easily found by this method. This method also gives a guarantee of stable solutions for the fractional order differential equations. Furthermore it gives asymptotically stable.

An evaluation was achieved by designing a matlab program to solve Kepler’s equation of an elliptical orbit for methods (Newton-Raphson, Danby, Halley and Mikkola). This involves calculating the Eccentric anomaly (E) from mean anomaly (M=0°-360°) for each step and for different values of eccentricities (e=0.1, 0.3, 0.5, 0.7 and 0.9). The results of E were demonstrated that Newton’s- Raphson Danby’s, Halley’s can be used for e between (0-1). Mikkola’s method can be used for e between (0-0.6).The term  that added to Danby’s method to obtain the solution of Kepler’s equation is not influence too much on the value of E. The most appropriate initial Gauss value was also determined to