function [qrz, p]=QRZerlegung(a)
    r=a
    n=size(a,1)
    epsilon=10^-16
    for i=1:n-1
        v=r(i:n,i)+norm(r(i:n,i)).*eye(n-i+1,1)
        qrz(i:n,i)=v
        Qhelp=eye(n-i+1,n-i+1)-(2./(norm(v)^2)).*(v*v')
        Qhelp2=eye(n,n)
        Qhelp2(i:n,i:n)=Qhelp
        r=Qhelp2*r
        p(i,1)=r(i,i)
        qrz(i,i+1:n)=r(i,i+1:n)
    end
    p(n,1)=r(n,n)
endfunction