def polder_function(self, x, y): s = .5 * np.exp(2 * x) * erfc(x / y + y) + \ .5 * np.exp(-2 * x) * erfc(x / y - y) return s