In this paper, the auxiliary equation method is used to study the Weierstrass elliptic function solutions and degenerate solutions of the variable coefficient higher order Schrödinger equation, including Jacobian elliptic function solutions, trigonometric function solutions and hyperbolic function solutions. The types of solutions of the variable coefficient higher-order Schrödinger equation are enriched, and the method of seeking precise and accurate solutions is extended. It is concluded that the types of degenerate solutions are related to the coefficients of the equation itself when the degenerate solutions are obtained from the solutions of the Weierstrass elliptic functions. In addition, the solutions form of the equation is extended from the power series expansion form to the Laurent series expansion form, and the corresponding solutions are obtained. After the conversion formula between the Weierstrass elliptic function solutions and the Jacobian elliptic function solutions is constructed, the Jacobian elliptic function solutions of the higher order Schrödinger equation with variable coefficients are also obtained. These have not been previously studied.