added function for computing the saturationa the front using the integral method