In this paper, a cubic Hermite spline interpolating scheme reproducing both linear polynomials and hyperbolic functions is considered. The interpolating scheme is mainly defined by means of integral values over the subintervals of a partition of the function to be approximated, rather than the function and its first derivative values. The scheme provided is C2 everywhere and yields optimal order. We provide some numerical tests to illustrate the good performance of the novel approximation scheme.