
for N from 5 to 10 do   
  n := 2*N;
  X := [seq(x||i,i=1..N)]:

  B := rand(2);
  PM[1]:=7*x2^2+3*x3*x4:
  PM[2]:=x2*x4*x5:
  PM[3]:=6*x3^2:
  PM[4]:=5*x2*x4+2*x3*x4:
  PM[5]:=9*x3*x5+5*x4^2:
  PM[6]:=7*x5^2:
  PM[7]:=4*x3*x6+10*x4*x5:
  PM[8]:=10*x6^3*x2:
  PM[9]:=3*x7*x6:
  PM[10]:=x3*x7^2+2*x2^2:

  QM[1]:=7*x3^2:
  QM[2]:=18*x4*x2*x3 + 31*x4*x3:
  QM[3]:=x2*x3:
  QM[4]:=10*x2*x4 + 2*x3^2:
  QM[5]:=94*x2^2*x5 + 41*x3*x4^2:
  QM[6]:=42*x2^2*x5:
  QM[7]:=14*x5^3*x6 + 93*x3*x6^2:
  QM[8]:=15*x2^3*x6:
  QM[9]:=3*x7^2:
  QM[10]:=5*x2*x4^2+3*x7*x6:

  T1 := LinearAlgebra:-ToeplitzMatrix([X[1],seq(PM[i],i=1..N-1)],symmetric):
  T1 := convert(T1,listlist):
  T1 := Matrix(T1):

  T2 := LinearAlgebra:-ToeplitzMatrix([X[1],seq(QM[i],i=1..N-1)],symmetric):
  T2 := convert(T2,listlist):
  T2 := Matrix(T2):

  Id := Matrix(N,N):
  AA0 := <<T1|Id>,<Id|T2>>;

  P1 := Matrix(2*N,2*N):
  for i to 2*N do P1[i,i] := 1; od:
  for i to 2*N do 
    for j from i+1 to 2*N do
      P1[i,j] := B();
    od;
  od;
  AA:=AA0.P1;    
   
  P2 := Matrix(2*N,2*N):
  for i to 2*N do P2[i,i] := 1; od:
  for i to 2*N do 
    for j from 1 to i-1 do
      P2[i,j] := B();
    od;
  od;
  AA:=AA.P2:
  
  saver := subs(FILENAME=Test||N, proc() save AA0,T1,T2,P1,P2,AA,FILENAME; end):
  saver():

od:
