Test whether a sample mean differs from a hypothesised when the population standard deviation is known:
Implement z_test(x, mu0, sigma, alternative="two-sided") returning (z, p_value), with alternative in "two-sided", "greater", "less". Use math.erf for the normal CDF.