Comparative Study of Matlab ODE Solvers for the Korakianitis and Shi Model

Abstract

Changing parameters of the Korakianitis and Shi heart valve model over a cardiac cycle has led to the investigation of appropriate numerical technique(s) for good speed and accuracy. Two sets of parameters were selected for the numerical test. For the seven MATLAB ODE solvers, the computed results, computational cost and execution time were observed for varied error tolerance and initial time steps. The results were evaluated with descriptive statistics; the Pearson correlation and ANOVA at α0.05. The dependence of the computed result, accuracy of the method, computational cost and execution time of all the solvers, on relative tolerance and initial time steps were ascertained. Our findings provide important information that can be useful for selecting a MATLAB ODE solver suitable for differential equation with time varying parameters and changing stiffness properties.