function A=interpoladd(x,y)
A=[x(1) y(1)];
n=length(x);
for k=2:n
    A=addpoint(x(k),y(k),A);
end
