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

