Construct M matrix

construct_M(meioses, mu_dw, mu_up)

Arguments

meioses

number of meioses separating the two individuals

mu_dw

mutation rate for 1-step down-mutation

mu_up

mutation rate for 1-step up-mutation