|
 |
 |
 |
 |
 |
 |
 |
 |
 |
D.3.1.9 power
Procedure from librarymatrix.lib (see matrix_lib).
- Usage:
- power(A,n); A a square-matrix of type intmat or matrix, n=integer>=0
- Return:
- intmat resp. matrix, the n-th power of A
- Note:
- for A=intmat and big n the result may be wrong because of int overflow
LIB "matrix.lib"; intmat A[3][3]=1,2,3,4,5,6,7,8,9; print(power(A,3));""; ==> 468 576 684 ==> 1062 1305 1548 ==> 1656 2034 2412 ==> ring r=0,(x,y,z),dp; matrix B[3][3]=0,x,y,z,0,0,y,z,0; print(power(B,3));""; ==> yz2, xy2+x2z,y3+xyz, ==> y2z+xz2,yz2, 0, ==> y3+xyz, y2z+xz2,yz2 ==> |