We propose an approach to calculate asymptotic series for low lying eigenvalues of Schrödinger operator based on normal forms in the formal graded Weyl–Heisenberg algebra. The difference from a traditional scheme is that we don’t use any symbol map (Weyl,
,
, etc.). We show that our method may be useful for different reasons. Firstly, it enables to estimate the growth of the eigenvalues expansion coefficients, and secondly it may be efficient for practical calculations, e.g. for treating inverse problems. In particular, we prove that under some restrictions in the one-dimensional case the knowledge of asymptotic series for any pair of low lying eigenvalues is enough to recover the potential.