double integration method