#!/usr/bin/perl -w
use strict;
use PDL;

# A 3 by 3 matrix
my $A = pdl [[1,2,3],[-1,0,-3], [5,6,7]];

# An identity matrix of size 3
my $B = identity 3;
print $B;

# Simple addition
my $C = $A + $B;

# Invert the $C matrix
my $D = inv $C;

print $C, $D;
