Create two_locus_model_multinomial.R

initial version