[469192798]

设λ=1, i(0)=0.1
>> y=dsolve('Dy=y-y^2','y(0)=0.1')

y =
1/(1+9*exp(-t))

>> x=0:0.01:13;
y=1./(1+9.*exp(-x));
>> plot(x,y)
title('SI模型的i~t曲线');
xlabel('t');
ylabel('i');
axis([0 13 0 1.1]);
将这个算法改为使用龙格库塔方法求解并画出图