tested integrator - pretty slow for high precision (need adaptive functions to improve further)