qrfa.m

Numerical analysis （Gauss、Bisection、jacobi、Lagrange、Newton、powermethod等）

function s=QRfa(a,N)
for k=1:N
    for i=1:3
        for j=i+1:4
            d=sqrt(a(i,i).*a(i,i)+a(j,i).*a(j,i));
            %disp(d);
            c=a(i,i)./d;
            s=a(j,i)./d;
            Q=eye(4);
            Q(i,i)=c;
            Q(j,j)=c;
            Q(i,j)=s;
            Q(j,i)=-s;
            %disp(Q);
            R=inv(Q);
            a=Q*a*R;
            %disp(a);
        end
    end
end
s=a;
return