Select Code to Download

use PDL:
sub p{$t=pdl+pop;$t=$t*(pdl$_)for@_;[list$t]}

print join " ", @{p([1..5],[1..5])}; 
1 4 9 16 25

sub p{my $t=pdl+pop;$t=$t*(pdl$_)for@_;[list$t]}

A = p1'*p2;
product = crossdiags(A);

sub p{$t=pdl+pop;for(@_){
  my $a=outer$t,pdl+pop;
...
  $t->slice("$_:".($_+$n))+=$a->slice("$_,:")for 0..$n;
}
[list$t]}