out ajj;
c:=-mat((2,-1,0),(-1,2,-1),(0,-1,2));
%c:=-mat((2,-1,0,0),(0,-1,2,-1),(0,-1,2,-1),(0,0,-1,2));
%c:=-mat((2,-1,0),(-1,2,-2),(0,-1,2));

%qq:=mat((1,-1,0),(0,1,-1),(0,0,1));
%q:=mat((-1,0,0),(1,-1,0),(0,1,-1));

%c:=mat((2,-1,0,0),(1,2,-2,0),(0,-1,2,-1),(0,0,-1,2));
%ff:=mat((0,1,-1,0),(0,0,1,-1),(0,0,0,1),(1/2,-1/2,-1/2,-1/2));
%f:=mat((0,0,0,-1/2),(-1,0,0,1/2),(1,-1,0,1/2),(0,1,-1,1/2));
%ff*f;

%ff:=mat((1/2,-1/2,-1/2,-1/2),(0,1,-1,0),(0,0,1,-1),(0,0,0,1));


%qq:=mat((-3,0),(1,-2));
%q:=mat((-1,2),(1,-1));

%c:=mat((2,-1),(-1,2));
%dd:=mat((d1,d2),(d2,d1));
%d:=mat((d1,0,0),(0,d1,0),(0,0,d1));
%d:=mat((d1,d2,d3),(d4,d5,d6),(d7,d8,d9));
d:=mat((d1,d2,d3),(d2,d1+d3,d2),(d3,d2,d1));
%d:=mat((d1,d2,d3,d4),(d5,d6,d7,d8),(d9,d10,d11,d12),(d13,d14,d15,d16));
%d:=mat((2*d9,d2,2*d7),(d2,d7+d9,2*d2),(d7,d2,d9));
j:=(1/2)*mat((1,sqrt(2),-1),(sqrt(2),0,sqrt(2)),(-1,sqrt(2),1));
jj:=mat((-3,0,0),(0,-4,0),(0,0,-3));
%j*c*(1/j);
j*jj*(1/j);

shut ajj;
bye;
end;