def dfdx(self, x): sech_aR = 1.0 / np.cosh(self.a * (x - self.R0)) return self.V0*self.a*sech_aR*sech_aR