test.in

setrand(1);
Alea(n,m)=
{ my(Q);
  Q = matrix(n,n,i,j, random(m)-m\2);
  return((Q+Q~) \ 2);
}
p = default(path);
default(path, Str(p, ":./libre"));
default(parisize, "100M");

\r libre/Qfmodife.gp
\r programme/Qfn.gp
\r programme/Pinceau.gp
\r programme/groupe_orthogonal.gp
\r programme/Changements_de_bases.gp
\r qf13.gp

default(timer,1)

a = Alea(13,10);
b = Alea(13,10);

Qfn(a,b, 5)

a = Alea(13,10^5);
b = Alea(13,10^5);

Qfn(a,b, 5)

{
for(i=1,50,
  print(i,": ");
  a = Alea(13,10^3);
  b = Alea(13,10^3);
  Qfn(a,b, 5)
)
}