clear all;

a=1.78;b=0.56;c=0; 
F=@(t,y)[-a*y(1).*y(2)+c*(1-y(1)-y(2)); a*y(1).*y(2)-b*y(2)];


figure(1);
%Ir=[2;7;22;71;194;384;492;446;336;231;152;97;62;39];
%plot([0:13], Ir/1526)
%hold on;
[t,y]=ode45(F,[0 30],[1524/1526 2/1526]);
plot(t,y) 
hold off;

figure(2);
[S,I]=meshgrid(0:0.05:1,0:0.05:0.5);
Sp=-a*I.*S+c*(1-I-S);Ip=a*I.*S-b*I;
quiver(S,I,Sp,Ip);
hold on;
axis tight;
plot(y(:,1), y(:,2))
%hold on;
hold off;



