Adding SIA4 Integration method

The algorithm is based on Franz J. Vesely, Computational Physics: An Introduction
The method was tested with the simple pendulum, the energy is conserved and the
trajectory is similar to the other methods trajectories