#!/usr/bin/perl
# [id://1229234]
use warnings;
use strict;
use PDL;

sub matrix_power {
    my ($M,$n) = @_;
    if($n<0) {
        return matrix_power(1/$M, -$n);
    } elsif ($n==0) {
        return 1;
    } elsif ($n==1) {
        return $M;
    } elsif (0 == $n % 2) {
        return matrix_power($M x $M, $n/2);
    } else {
        return $M x matrix_power($M x $M, ($n-1)/2);
    }
}

my $x = pdl([0,1], [2,3]);
for(0 .. 3) {
    print "M ** $_ = " . matrix_power( $x, $_ ) . "\n";
}