In the present work, the notion of Cubic Spline Super Fractal Interpolation Function (SFIF) is introduced to simulate an object that depicts one structure embedded into another and its approximation properties are investigated. It is shown that, for an equidistant partition points of [x0, xN], the interpolating Cubic Spline SFIF[Formula: see text] and their derivatives [Formula: see text] converge respectively to the data generating function y(x) ≡ y(0)(x) and its derivatives y(j)(x) at the rate of h2-j+ϵ(0 < ϵ < 1), j = 0, 1, 2, as the norm h of the partition of [x0, xN] approaches zero. The convergence results for Cubic Spline SFIF found here show that any desired accuracy can be achieved in the approximation of a regular data generating function and its derivatives by a Cubic Spline SFIF and its corresponding derivatives.