#!/usr/bin/perl 

use PDL;

my $matrix = pdl [
	       [4, 2, 9],
	       [2, 5, 5],
	       [6, 3, 5],
	       [8, 4, 6],
	       ];

($a, $b, $c) = svd($matrix);

$bt = transpose($b);

$B1 = $bt*$b;

$B2 = $B1^-1;

$pseudoB = $B2*$bt;