We review the recent glueball mass calculations using an efficient met
hod for solving the Schrodinger equation order by order with a scheme
preserving the continuum limit. The reliability of the method is furth
er supported by new accurate results for (1+1)-dimensional a models an
d (2+1)-dimensional non-abelian model. We present first and encouragin
g data for the glueball masses in 3+1 dimensional QCD.