This paper mainly presents numerical solutions to an initial-boundary value problem of the time-fractional Klein-Gordon equations. We developed a numerical scheme with the help of the finite difference methods and the predictor-corrector methods to find numerical solutions of the considered problems. The proposed scheme is based on discretizing the considered problems with respect to spatial and temporal domains. Numerical results are derived for some illustrative problems, and the outputs are compared with the exact solution in the integer order case. The solution behavior and 3D graphics of the discussed problems are demonstrated using the proposed scheme. Finally, the proposed scheme, which does not require solving large systems of linear equations, can be extended and modified to handle other classes of time-fractional PDEs.