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

