Pk,n(X)=((Xk-X)Pn-1,n-1-(Xn-1-X)Pk,n-1)/(Xn-Xn-1)
dim x(22),y(255);
labei li;
if h>22 then print"too many datas" go l1;
for i:=1 to k;
input x(i);
input y(i*(i-1)/2+1);
next i;
i1:=2;
for i:=i1 to k;
for j:=2 to i;
i0:=i*(i-1)/2;
j0:=j-1;
y(i0+j)=((x(i)-x)*y(j0*(j0+1)/2)-(x(j0)-x)*y(i0+j0)))/(x(i)-x(j0));
next j;
next i;
end