Two classes of linearly implicit numerical methods for stiff problems: analysis and MATLAB software

The purpose of this work lies in the writing of efficient and optimized Matlab codes to implement two classes of promising linearly implicit numerical schemes that can be used to accurately and stably solve stiff Ordinary Differential Equations (ODEs), and also Partial Differential Equations (PDEs) through the Method Of Lines (MOL). Such classes of methods are the Runge-Kutta (RK) [28] and the Peer [17], and have been constructed using a variant of the Exponential-Fitting (EF) technique [27]. We carry out numerical tests to compare the two methods with each other, and also with the well known and very used Gaussian RK method, by the point of view of stability, accuracy and computational cost, in order to show their convenience.