A first approach in solving initial-value problems in ODEs by elliptic fitting methods

Exponentially-fitted and trigonometrically-fitted methods have a long successful history in the solution of initial-value problems, but other functions might be considered in adapted methods. Specifically, this paper aims at the derivation of a new numerical scheme for approximating initial value problems of ordinary differential equations using elliptic functions. The example considered is the undamped Duffing equation where the forcing term is of autonomous type affected by a perturbation parameter. The new scheme is constructed by considering a suitable approximation to the theoretical solution based on elliptic functions. The proposed elliptic fitting procedure has been tested on a variety of problems, showing its good performance.