t = [0,300];
sir0 = [1000;0.001];
[t,sir] = ode45('sir',t,sir0);
plot(t,sir);
xlabel('日');
ylabel('千人');
legend('S','I');

