Anfonwch hwn fel neges destun: A boundary preserving numerical algorithm for the Wright-Fisher model with mutation