Xi = self.X*self.dx + dmax * np.cos(angle-np.pi/2)