R=[2^0   2^1   2^2   2^3   2^4   2^5;
   2^1   2^0   2^1   2^2   2^3   2^4;
   2^2   2^1   2^0   2^1   2^2   2^3;
   2^3   2^2   2^1   2^0   2^1   2^2;
   2^4   2^3   2^2   2^1   2^0   2^1;
   2^5   2^4   2^3   2^2   2^1   2^0];
X0=ones(5,1);
Y=zeros(5,1);
K=zeros(5,1);
X0(1,1)=-R(1,2)/R(1,1);
Y(1,1)=R(1,1)-R(1,2)^2/R(1,1);
K(1,1)=X0(1,1)
for i=2:5
    sum=0;
    for j=1:i-1
        sum=sum+X0(j,1)*R(1,i-j+1);
    end
    K(i,1)=-(R(1,i+1)+sum)/Y(i-1,1);
    Y(i,1)=Y(i-1,1)*(1-K(i,1)^2);
    for j=1:i-1
        X1(j,1)=X0(j,1)+K(i,1)*X0(i-j,1);
    end
    X1(i,1)=K(i,1);
    for j=1:i
        X0(j,1)=X1(j,1);
    end
end