{OperatorNormTransformed(n,tol) =
    \\ Gewichte und Stützstellen der Summationsformel
    u = vector(n-1,k,k/n);
    c1 = vector(n-1,k,exp(2/(1-u[k])^2-1/u[k]^2/2)); c = vector(n-1,k,k+1) + c1;
    w = vector(n-1,k,1 + c1[k]*(4/(1-u[k])^3+1/u[k]^3)/n);
    c = concat([1],c); w = concat([1],w); w = sqrt(w);
    \\ Vektoriteration
    x = vector(n); x[1]=1;
    lambda=0; lambda0=1; 
    while(abs(lambda-lambda0)>tol, 
        xhat=x/sqrt(sum(k=-n,-1, x[-k]^2));
        y=vector(n,k, sum(l=-n,-1, w[k]*w[-l]*xhat[-l]/((c[k]+c[-l]-1)*(c[k]+c[-l])/2-(c[-l]-1))));
        x=vector(n,j, sum(k=-n,-1, w[j]*w[-k]*   y[-k]/((c[j]+c[-k]-1)*(c[j]+c[-k])/2-( c[j]-1))));
        lambda0=lambda; lambda=sum(k=-n,-1,xhat[-k]*x[-k]);
    ); 
    sqrt(lambda)
}