We develop an inverse matrix method to solve for resonance masses from a dispersion relation obeyed by correlation function. Given the operator product expansion (OPE) of function in deep Euclidean region, we obtain nonperturbative spectral density, which exhibits structures naturally. The value gluon condensate OPE is fixed producing $\rho$ meson mass formalism, and then input into relations s...