>>64 >>66の続き
int main(){
double R2 = sqrt(2), R3 = sqrt(3), R6 = sqrt(6);
double A [3][3] = {{2,1,1},{1,2,1},{1,1,2}};
double P [3][3] = {{1/R2, 1/R6, 1/R3},{-1/R2, 1/R6, 1/R3},{0,-2/R6,1/R6}};
double tP [3][3], tPA [3][3], tPAP[3][3];
int i,j;
show(A);
show(P);

inverse (tP ,P);
show(tP);

multiple(tPA ,tP ,A);
show(tPA);

multiple(tPAP,tPA,P);
show(tPAP);

}