计算矛盾方程组MatLab

%Xi=[1,3,4,5,6,7,8];
%Yi=[10,5,4,2,1,1,2];
Xi=[-2,-1,0,1,2];
Yi=[0,1,2,1,0];
N=size(Xi);
i=N([2]); %组数
%y=φ(x)=a0+a1x+a2X^2
for j=1:1:i %计算矛盾方程组
A(j,1)=1;
A(j,2)=Xi([j]);
A(j,3)=(Xi([j]))^2;
X(j,1)=0;
B(j,1)=Yi([j]);
end
a=A\B; %求解
a'