Extrapolating for attaining high precision solutions for fractional partial differential equations.

Abstract

[EN]This paper aims at obtaining a high precision numerical approximation for fractional partial differential equations. This is achieved through appropriate discretizations: firstly we consider the application of shifted Legendre or Chebyshev polynomials to get a spatial discretization, followed by a temporal discretization where we use the Implicit Euler method (although any other temporal integrator could be used). Finally, the use of an extrapolation technique is considered for improving the numerical results. In this way a very accurate solution is obtained. An algorithm is presented, and numerical results are shown to demonstrate the validity of the present technique.