Select Code to Download

    for my $aref ( @lnPAM ) {
        $_ /= ( $mutRate * 100 ) for @$aref;
    }

@a = map[ map int(rand$_), 1 .. 10],1..10;;
cmpthese -1,{
...
     Rate    b    a
b 52726/s   -- -25%
a 70388/s  33%   --

while( abs( $mutRate - 0.01) > 1e-16 ) {
    ++$iter;
...
    $mutRate = sum map{ $freq[ $_ ] * ( 1 - $matrix[ $_ ][ $_ ] ) } 0.
+.19;
    printf "After %d iteration(s), k=%f, RateMutation=%f\n", $iter, $k
+, $mutRate;
}

while abs(RateMutation - 0.01) > DBL_EPSILON do
    iter := iter + 1;
...
    RateMutation := sum(Freq[i] * (1 - PAM1[i,i]), i= 1..20);
    printf('After %d iteration(s), k=%f, RateMutation=%f', iter, k, Ra
+teMutation);
od: