The discretized Lyapunov functional method for the stability problem of time-delay systems is further refined using a combination of integral inequality and variable elimination technique. As a result, the computational requirement is further reduced for the same discretization mesh. For systems without uncertainty, the convergence to the analytical solution is greatly accelerated. For uncertain systems, the new method is much less conservative. Numerical examples are presented to illustrate the effectiveness of the method.