def mu_disc(x, x0= -3.0): integral, _ = quad(p_discontinuous, x0, x) return np.exp(integral)

There are several types of differential equations, including: