A Reliable Numerical Method for Computing Homogenized Coeecients