This paper investigates the recovery for time-dependent coefficient and free boundary for heat equation. They are considered under mass/energy specification and Stefan conditions. The main issue with this problem is that the solution is unstable and sensitive to small contamination of noise in the input data. The Crank-Nicolson finite difference method (FDM) is utilized to solve the direct problem, whilst the inverse problem is viewed as a nonlinear optimization problem. The latter problem is solved numerically using the routine optimization toolbox lsqnonlin from MATLAB. Consequently, the Tikhonov regularization method is used in order to gain stable solutions. The results were compared with their exact solution and tested via

     Transformation and many other substitution methods have been used to solve non-linear differential fractional equations. In this present work, the homotopy perturbation method to solve the non-linear differential fractional equation with the help of He’s Polynomials is provided as the transformation plays an essential role in solving differential linear and non-linear equations. Here is the α-Sumudu technique to find the relevant results of the gas dynamics equation in fractional order. To calculate the non-linear fractional gas dynamical problem, a consumer method created on the new homotopy perturbation a-Sumudu transformation method (HP TM) is suggested. In the Caputo type, the derivative is evaluated. a-Sumudu homotopy pe

A new method based on the Touchard polynomials (TPs) was presented for the numerical solution of the linear Fredholm integro-differential equation (FIDE) of the first order and second kind with condition. The derivative and integration of the (TPs) were simply obtained. The convergence analysis of the presented method was given and the applicability was proved by some numerical examples. The results obtained in this method are compared with other known results.

This paper considers a new Double Integral transform called Double Sumudu-Elzaki transform DSET. The combining of the DSET with a semi-analytical method, namely the variational iteration method DSETVIM, to arrive numerical solution of nonlinear PDEs of Fractional Order derivatives. The proposed dual method property decreases the number of calculations required, so combining these two methods leads to calculating the solution's speed. The suggested technique is tested on four problems. The results demonstrated that solving these types of equations using the DSETVIM was more advantageous and efficient

The goal of this paper is to expose a new numerical method for solving initial value time-lag of delay differential equations by employing a high order improving formula of Euler method known as third order Euler method. Stability condition is discussed in detail for the proposed technique. Finally some examples are illustrated to verify the validity, efficiency and accuracy of the method.

This paper introduces a non-conventional approach with multi-dimensional random sampling to solve a cocaine abuse model with statistical probability. The mean Latin hypercube finite difference (MLHFD) method is proposed for the first time via hybrid integration of the classical numerical finite difference (FD) formula with Latin hypercube sampling (LHS) technique to create a random distribution for the model parameters which are dependent on time [Formula: see text]. The LHS technique gives advantage to MLHFD method to produce fast variation of the parameters’ values via number of multidimensional simulations (100, 1000 and 5000). The generated Latin hypercube sample which is random or non-deterministic in nature is further integ

The approach given in this paper leads to numerical methods to find the approximate solution of volterra integro –diff. equ.1st kind. First, we reduce it from integro VIDEs to integral VIEs of the 2nd kind by using the reducing theory, then we use two types of Non-polynomial spline function (linear, and quadratic). Finally, programs for each method are written in MATLAB language and a comparison between these two types of Non-polynomial spline function is made depending on the least square errors and running time. Some test examples and the exact solution are also given.