The stability problem for systems with distributed delay is considered using discretized Lyapunov functional. The coefficients associated with the distributed delay are assumed to be piecewise constant, and the discretization mesh may be non-uniform. The resulting stability criteria are written in the form of linear matrix inequality. Numerical examples are also provided to illustrate the effectiveness of the method. The basic idea can be extended to a more general setting with more involved formulation.