A REDUCE program based on the two step approach is developed to solve the Schrodinger equation with polynomial potential. The program is very simple and suitable to be carried out on microcomputer, The transformation for Hamiltonian, the obtaining of matrix elements, the derivation of the equation which must be satisfied by the parameter in variation method and even the Fortran program for numerical calculation are all completed or generated by computer automatically, As an example the pure quartic oscillator is calcuated and the results are compared with those available in literatures, Agreements are excellent.