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