function [bdiag, odiag] = bdiag_odiag(RT,n);

% 14.02.2002, rsotirov, klagenfurt
% testing wether b^0diag & o^0diag is fulfiled
% call: [bdiag, odiag] = bdiag_odiag(Z);

  
joj = 1:n;
bdiag = zeros(n);

for i=1:n,
  bdiag = bdiag + RT(n*(i-1) + 1+joj, n*(i-1) + 1+joj );
end;


odiag =  zeros(n);
for i=1:n,
odiag = odiag + RT((joj-1)*n +1 +i ,(joj-1)*n +1 + i);
end;