S0 = self.smgauss(X, S0, DX)