function [res] = INTE(J,JJ) % This function is used by Cubic Spline Rayleigh Ritz Method % (Algorithm 11.6). res = JJ-J+3;