%
% �¶‘ÌƒVƒXƒeƒ€�HŠw‚ÌŠî‘b
% �}4.74
%
clear all
close all
dt=0.01;
step=10000;
alpha = 0.4;
beta = 1.0;
m(1)=10;
e(1)=0;
for i=2:step
    m(i) = m(i-1) + dt * (e(i-1) * e(i-1)/ (1 + e(i-1) * e(i-1)) - alpha * m(i-1));
    e(i) = e(i-1) + dt * (m(i-1) - beta * e(i-1));
end
plot(m, e);
xlabel('M');
ylabel('E');
xlim([0 10]);
