simpson_integrator -> c++