clear
//ＤＣＴ変換行列
n=8;
xd=(1:n);yd=(1:n);

for i = 1:n,
  for j = 1:n, 
      if i==1 then c=1/sqrt(2); else c=1;end;
      a(i,j) = sqrt(2/n)*c*cos((i-1)*(j-1/2)*%pi/n);
  end;
end;
xbasc();plot3d(xd-1,yd-1,a(xd,yd)); 
x=[1:1:n]',
y=a*x,
inv(a)-a', //ユニタリー性のチェック