clear;

[U,X]=qr(rand(80,80));
[V,X]=qr(rand(80,80));

S=diag((2.0).^(-1:-1:-80));
A=U*S*V';

exec("mGS.sci")
exec("clGS.sci")

[QC,RC] = clGS(A);
[QM,RM] = mGS(A);
//[Q,R] = qr(A);
ii=(1:1:80)';
xbasc()
plot2d(ii,[diag(RC), diag(RM), diag(S)],...
          style=[-1 2 3],...
          leg="Klasicni GS@Modificirani GS@Singularne vrijednosti",...
          logflag="nl")