Double integration method