%
% �¶‘ÌƒVƒXƒeƒ€�HŠw‚ÌŠî‘b
% �}4.71
%
close all
clear all
dt=0.01;
step=10000;
x(1)=2;
y(1)=2;
for i=2:step
    x(i) = x(i-1) + dt * (1 - x(i-1) * y(i-1) * y(i-1));
    y(i) = y(i-1) + dt * (x(i-1) * y(i-1) * y(i-1) - y(i-1));
end
plot(x, y);
xlabel('x');
ylabel('y');
xlim([0 4]);
ylim([0 4])
