Second-order linear hyperbolic equations are solved by using a new three level method based on non-polynomial spline in the space direction and Taylor expansion in the time direction. Numerical results reveal that three level method based on non-polynomial spline is implemented and effective.