|
Abstract
|
In this paper, we present a numerical method to solve a linear fractional differential equations. This new investigation is based on ultraspherical integration matrix to approximate the highest order derivatives to the lower order derivatives. By this approximation the problem is reduced to a constrained optimization problem which can be solved by using the penalty quadratic interpolation method. Numerical examples are included to confirm the efficiency and accuracy of the proposed method.
|