add a script that tries to find max likelihood estimates of parameters of a codon model