计算矛盾方程组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'