removed extra parameters to substitution matrix fn

The extra parameters were usually set to 1.0, or didn't matter since the
mean-rate of the matrices are divided through to equal 1.0. There are some
improvements still on the functions that could avoid the division of the
constant for JC69 and K80 and F81, but currently these models are correct.