In this paper, Crank-Nicholson method for solving fractional wave equation is considered. The stability and consistency of the method are discussed by means of Greschgorin theorem and using the stability matrix analysis. Numerical solutions of some wave fractional partial differential equation models are presented. The results obtained are compared to exact solutions.