function a=lsquad(x,y)
% [SUMx][a]=[B]
n=length(x); x=x(:); y=y(:);
sumx=sum(x);
sumx2=sum(x.^2);
sumx3=sum(x.^3);
sumx4=sum(x.^4);
sumy=sum(y);
sumxy=sum(x.*y);
sumx2y=sum(x.*x.*y);
SUMx=[n      sumx    sumx2;
     sumx   sumx2   sumx3;
     sumx2  sumx3   sumx4];
B=[sumy   sumxy   sumx2y]';

%Cebirsel denklem takımını çözmek için aşağıdaki yöntemlerden birini
%   seçebilirsiniz
%a=uptrbk(SUMx,B)%Gaus eliminasyon ile denklem sisteminin çözümü
%a=gseidel(SUMx,B,0,0.0001,1000)%Gauss Seidel ile çözüm
%a=gseidel(SUMx,B,1.5,1.5,1.5)%Gauss Seidel ile çözüm
%a=inv(SUMx)*B %Ters matris yöntemi ile çözüm
a=SUMx\B %Denklem sistemini direkt çözen MATLAB komutu