has been extensively treated in the literature. Consider the pure quartic oscillator, in which the quadratic term is missing: . For simplicity, take . The Schrödinger equation thus reduces to

.

No analytic solution has been found, but accurate numerical computations have been carried out [1]. This Demonstration applies the operator method [2], a generalization of the canonical operator formulation for the harmonic oscillator, for potential energies that are functions of even powers of . For larger values of the matrix dimension, the results obtained here are comparable with the published results.

The computation results in a secular equation for the eigenvalues, which are plotted as red lines on the graph superposed on the potential energy curve. The number of eigenvalues shown is equal to the selected matrix dimension. For comparison, the first eight eigenvalues according to the WKB method are shown as thin gray lines.