Quasilinearization method and summation of the WKB series

Solutions obtained by the quasilinearization method (QLM) are compared with the WKB solutions. Expansion of the pth QLM iterate in powers of h reproduces the structure of the WKB series generating an infinite number of the WKB terms with the first 2^p terms reproduced exactly. The QLM quantization condition leads to exact energies for the Pöschl-Teller, Hulthén, Hylleraas, Morse, Eckart potentials, etc. For other, more complicated potentials the first QLM iterate, given by the closed analytic expression, is extremely accurate. The iterates converge very fast. The sixth iterate of the energy for the anharmonic oscillator and for the two-body Coulomb-Dirac equation has an accuracy of 20 significant figures.