OBJECTIVE: The main purpose of this work is to present cubic non-polynomial spline approximation method for solving Robin-type singularly perturbed reaction-diffusion problems. RESULTS: The solution domain is first discretized using a piecewise mesh. The process begins by defining the cubic non-polynomial spline function and calculating its derivatives. These derivatives are then transformed into difference approximations, forming a linear system of algebraic equations in the form of a three-term recurrence relation, which is solved using an elimination algorithm. The stability and consistency of the method are analyzed, ensuring convergence. Numerical model examples are used to validate the proposed method, and the results are compared with those from other methods found in the literature. The maximum absolute error and the order of convergence for each example demonstrate the effectiveness and core contribution of the method.