In this paper, we considered an inverse problem of recovering the space-dependent source coefficient in the third-order pseudo-parabolic equation from final over-determination condition. This inverse problem appears extensively in the modelling of various phenomena in physics such as the motion of non-Newtonian fluids, thermodynamic processes, filtration in a porous medium, etc. The unique solvability theorem for this inverse problem is supplied. However, since the governing equation is yet ill-posed (very slight errors in the final input may cause relatively significant errors in the output source term), we need to regularize the solution. Therefore, to get a stable solution, a regularized cost function is to be minimized for retrieval of the unknown force term. The third-order pseudo-parabolic problem is discretized using the Cubic B-spline (CB-spline) collocation technique and reshaped as non-linear least-squares optimization of the Tikhonov regularization function. Numerically, this is effectively solved using the lsqnonlin routine from the MATLAB toolbox. Both perturbed data and analytical solutions are inverted. Numerical outcomes are reported and discussed. The computational efficiency of the method is investigated by small values of CPU time. In addition, the von Neumann stability analysis for the proposed numerical approach has also been discussed.