Matrix-exponentiation operator: Difference between revisions

=={{header|Scheme}}==

For simplicity, the matrix is represented as a list of lists, and no dimension checking occurs.

<lang scheme>

(define (row*col row col)

(apply + (map * row col)))

(define (multiply-matrix m1 m2)

(map

(lambda (row)

(apply map (lambda col (row*col row col))

m2))

m1))

(define (matrix-expo mat exp)

(cond ((= exp 1) mat)

((even? exp) (square-matrix (matrix-expo mat (halve exp))))

(else (multiply-matrix mat (matrix-expo mat (dec exp))))))

</lang>