clear;
clf;
function dX = LC(t,X) // X(1) = x, X(2) = y
dX=zeros(2,1);
dX(1) = 10 - X(1) - 4 * X(1) * X(2) / (1 + X(1)^2) ;
dX(2) = X(1) * (1 - X(2) / (1 + X(1)^2)) ;
endfunction
xc=-1:1:11;
yc=-10:10:110;
fchamp(LC,0,xc,yc,1)
T = linspace(0,100,100000);
x01 = [0 ; 0];
x1 = ode(x01,0,T,LC);
x02 = [1.0143 ; 4.06];
x2 = ode(x02,0,T,LC);
x03 = [10 ; 0];
x3 = ode(x03,0,T,LC);
x04 = [10 ; 101];
x4 = ode(x04,0,T,LC);
x05 = [10 ; 50];
x5 = ode(x05,0,T,LC);
x06 = [10 ; 20];
x6 = ode(x06,0,T,LC);
plot(x1(1,:),x1(2,:),x2(1,:),x2(2,:),x3(1,:),x3(2,:),x4(1,:),x4(2,:),x5(1,:),x5(2,:),x6(1,:),x6(2,:));
zoom_rect([-1,-10,11,110])
xlabel('x');
ylabel('y');
gh=gca();
gh.font_style=2;
gh.font_size=5;
gh.x_label.font_style=3;
gh.x_label.font_size=6;
gh.y_label.font_style=3;
gh.y_label.font_size=6;