function Sep16()

function delta = Max3Jump(S)
% S is a pp-represenattion of a spline S(x) with breakpoints
% x(1)<...<x(n). I is the integral of S from x(1) to x(n)
% Let
% alpha(k) = 3rd derivative of the kth local cubic.
% delta = max(abs(alpha(2:n-1)-alpha(1:n-2)))
[x,Coef] = unmkpp(S);
n = length(x);
for k=1:n-1
alpha(k) =
end
delta = max(abs(alpha(2:n-1)-alpha(1:n-2)));