Summary: We propose a three level implicit unconditionally stable difference scheme of \(O(k^2+ h^2)\) for the difference solution of second-order linear hyperbolic equation \(u_{tt}+ 2\alpha(x, t)u_t+ \beta^2(x, t)u= A(x, t)u_{xx}+ f(x,t)\), \(0< x< 1\), \(t> 0\) subject to appropriate initial and Dirichlet boundary conditions, where \(A(x, t)> 0\), \(\alpha(x, t)>\beta(x, t)\geq 0\). The proposed formula is applicable to the problems having singularity at \(x= 0\). The resulting tri-diagonal linear system of equations is solved by using the Gauss-elimination method. Numerical examples are provided to illustrate the unconditionally stable character of the proposed method.