Some modified techniques are used in this article in order to have approximate solutions for systems of Volterra integro-differential equations. The suggested techniques are the so called Laplace-Adomian decomposition method and Laplace iterative method. The proposed methods are robust and accurate as can be seen from the given illustrative examples and from the comparison that are made with the exact solution.

     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.

 In this paper, we introduce and discuss an algorithm for the numerical solution of two- dimensional fractional partial differential equation with parameter. The algorithm for the numerical solution of this equation is based on implicit and an explicit difference method. Finally, numerical example is provided to illustrate that the numerical method for solving this equation is an effective solution method.

In this paper, third order non-polynomial spline function is used to solve 2nd kind Volterra integral equations. Numerical examples are presented to illustrate the applications of this method, and to compare the computed results with other known methods.

This paper sheds the light on the vital role that fractional ordinary differential equations(FrODEs) play in the mathematical modeling and in real life, particularly in the physical conditions. Furthermore, if the problem is handled directly by using numerical method, it is a far more powerful and efficient numerical method in terms of computational time, number of function evaluations, and precision. In this paper, we concentrate on the derivation of the direct numerical methods for solving fifth-order FrODEs  in one, two, and three stages. Additionally, it is important to note that the RKM-numerical methods with two- and three-stages for solving fifth-order ODEs are convenient, for solving class's fifth-order FrODEs. Numerical exa

A particular solution of the two and three dimensional unsteady state thermal or mass diffusion equation is obtained by introducing a combination of variables of the form,
η = (x+y) / √ct , and η = (x+y+z) / √ct, for two and three dimensional equations
respectively. And the corresponding solutions are,
θ (t,x,y) = θ₀ erfc (x+y)/√8ct and θ( t,x,y,z) =θ₀ erfc (x+y+z/√12ct)