function Sep2()

% Generate some data and the interpolant
n = 5;
x = linspace(0,pi,n);
y = sin(x);
c = polyfit(x,y,4);
%   p(x) = c(1)x^4 + c(2)x^3 + c(3)x^2 + c(4)x + c(5)
%  dp(x) = 4c(1)x^3 + 3c(2)x^2 + 2c(3)x + c(4)
%        = b(1)x^3 + b(2)x^2 + b(3)x + b(4)
b = zeros(n-1,1)
for k=1:n-1
    b(k) =   
end
xVals = linspace(x(1),x(n))
dpVals = polyval(b,xVals)
plot(xVals,dpVals)
shg


function a = InterpV(x,y,u,v)
n = length(x);
w = (x-u)/v;
V = zeros(n,n);
for k=1:n
    V(:,k) =
end
a = V\y;
