Don't divide by zero when theta is much smaller than 1.
Add test13.hs

This file is for comparing the output of diploid to haploid analysis of
diploid data.