def METHOD_NAME(self, ts):
    ts = ts * self.x_range
    dxs = self.x_direction[np.newaxis].T
    ys = np.sinh((ts - self.x0) / self.A)
    dys = - ys * self.force[np.newaxis].T
    return (dxs + dys).T