In this paper, author’s study sub diffusion bio heat transfer model and developed explicit finite difference scheme for time fractional sub diffusion bio heat transfer equation by using caputo fabrizio fractional derivative. Also discussed conditional stability and convergence of developed scheme. Furthermore numerical solution of time fractional sub diffusion bio heat transfer equation is obtained and it is represented graphically by Python.

In this paper Volterra Runge-Kutta methods which include: method of order two and four will be applied to general nonlinear Volterra integral equations of the second kind. Moreover we study the convergent of the algorithms of Volterra Runge-Kutta methods. Finally, programs for each method are written in MATLAB language and a comparison between the two types has been made depending on the least square errors.

This study is concerned with the derivation of differential equation of motion for the free coupled vertical – torsional and lateral vibration of opened thin-walled curved beams. The curved beam to be considered in this study is of isotropic opened thin – walled (I) section with equal top and bottom flanges. The derivation depends on Hamilton's principle which required finding the potential and kinetic energy of the curved beam section due to internal stresses and all types of movements (Vertical,Torsional and Lateral) .The effect of restrained warping displacement is also considered in this study. Three differential equations are derived for vertical, torsional and lateral movement .and approximate solutions are developed by using the

In this paper we use Bernstein polynomials for deriving the modified Simpson's 3/8 , and the composite modified Simpson's 3/8 to solve one dimensional linear Volterra integral equations of the second kind , and we find that the solution computed by this procedure is very close to exact solution.

     The method of operational matrices is based on the Bernoulli and Shifted Legendre polynomials which is used to solve the Falkner-Skan equation. The nonlinear differential equation converting to a system of nonlinear equations is solved using Mathematica®12, and the approximate solutions are obtained. The efficiency of these methods was studied by calculating the maximum error remainder ( ), and it was found that their efficiency increases as  increases. Moreover, the obtained approximate solutions are compared with the numerical solution obtained by the fourth-order Runge-Kutta method (RK4), which gives  a good agreement.

Lagrange series and the Bessel function are two classical methods that were created by series expanding from Taylor series. In this paper, the purpose of those two methods was to find the values of the eccentric anomaly for one period (0–360)°. The Matlab program is used to apply the results, the input parameters were eccentricity (0–1), mean anomaly (0–360)°, and finally the parameter W (1–13). The program does not need a tolerance to obtain a precise value for eccentric anomaly like other iterative and non-iterative methods to stop the program; it will stop after completing the required period from 0° to 360° for a body that is determined by the solver. The output will be the final value of the eccentric anomaly. Furthermore,