Added gauss_magnus.c This is a standalone (almost, still relies on cubature)
that calculated the third order Magnus correction for the special case of gaussian
functions for th phase matching and pump function.