function xdot=fdot2DKdVB(xi,t)

global  mu nu  xe1 xe2


x=xi(1);y=xi(2);
xdot=zeros(2,1);


 xdot(1)=y;
 xdot(2)=(nu/mu)*y - 1/(2*mu)*(x-xe1)*(x-xe2);
 
