function xdot=Hahn_eq(t,x)
xdot=[(x(1)^2*(x(2)-x(1))+x(2)^5)/((x(1)^2+x(2)^2)*(1+(x(1)^2+x(2)^2)^2));
      (x(2)^2*(x(2)-2*x(1)))/((x(1)^2+x(2)^2)*(1+(x(1)^2+x(2)^2)^2))];