function sf=sassenfeld(a)
% Sassenfeldwert = obere Schranke der Zs-Norm der ESV-Iterationsmatrix
[n m]=size(a);
s=zeros(n,1);
for i=1:n
    u=abs(a(i,1:i-1))*s(1:i-1,1);
    v=sum(abs(a(i,i+1:end)));
    s(i,1)=(u+v)/abs(a(i,i));
end
sf=max(s);


