Maxima/行列操作
< Maxima
行列のサイズ 編集
size(Z):=block(fn(f):=1,C:matrixmap(fn,col(Z,1)),R:matrixmap(fn,row(Z,1)),
return(matrix([transpose(C).C,R.transpose(R)])))$
a:matrix([1,2],[3,4],[5,6]);size(a);
対角行列の作成 編集
size(Z):=block(fn(f):=1,C:matrixmap(fn,col(Z,1)),R:matrixmap(fn,row(Z,1)),
return(matrix([transpose(C).C,R.transpose(R)])))$
diag(v):=block(n:size(listarray(y))[1,2],
P:zeromatrix(n,n),
FOR i:1 THRU n STEP 1 DO P[i,i]:y[i],
return(P))$
y[1]:1;y[2]:2;y[3]:3;diag(y);