In this present work, a new two-level implicit non-polynomial spline in tension method is proposed for the numerical solution of the 1D unsteady quasi-linear biharmonic problem. Using the continuity of the first-order derivative of the spline in tension function, a fourth-order accurate implicit finite-difference method is developed in this manuscript. By considering the linear model biharmonic problem, the given implicit spline method is unconditionally stable. Since the proposed method is based on half-step grid points, so it can be directly applied to 1D singular biharmonic problems without alteration in the scheme. Finally, the numerical experiments of the various biharmonic equation such as the generalized Kuramoto-Sivashinsky equation, extended Fisher-Kolmogorov equation, and 1D linear singular biharmonic problems are carried out to show the efficacy, accuracy, and reliability of the method. From the computational experiments, improved numerical results obtained as compared to the results obtained in earlier research work.