Perl version of Nilsimsa code