3 % => Ergebnis RB
' * G * RB = eye;
9 check=RB'*G*RB; %should be eye(NOS);
10 %trace_check=trace(check); %should be NOS
11 %matrix_check=sum(sum(check)); %should be NOS
12 err = max(max(abs(check - eye(NOS))));
function res = RBcheck(RB, G, NOS)
=> Ergebnis RB' * G * RB = eye; check result