Select Code to Download

PAM1 := MutationMatrix:
lnPAM1 := ln(MutationMatrix):
...
    printf('After %d iteration(s), k=%f, RateMutation=%f
', iter, k, RateMutation);
od:

After 1 iteration(s), k=44.530219, RateMutation=0.014324
After 2 iteration(s), k=63.785069, RateMutation=0.010027
...
After 7 iteration(s), k=63.960992, RateMutation=0.010000
After 8 iteration(s), k=63.960992, RateMutation=0.010000

my @matrix =(
[0.54, 0.03, 0.03, 0.03, 0.08, 0.04, 0.05, 0.05, 0.02, 0.02, 0.03, 0.0
+4, 0.03, 0.01, 0.06, 0.11, 0.06, 0.01, 0.01,0.06],
...
    
}