clear all
close all
clc

x = 0:11;
K = 50000;
a = 5;
b = -.8;
yR = [0.0485;
    0.0142;
    0.1084;
    0.3163;
    0.7735;
    1.2478;
    2.1594;
    3.3621;
    3.6407;
    4.5476;
    4.7924;
    4.8808]';
plot(x + 1981,yR,'r*','MarkerSize',10,'LineWidth',2)
hold on
grid on
xlabel('Years')
ylabel('Number of New Cases (in 10K)')
title('Number of New AIDS Cases in the US Since 1981')
axis([1980 1993 0 5.5])
display('Press Any Key To Continue')
pause


sol = ode45(@odeFun,[0 20], .0485*10000);
x = 0:.01:20;
y = deval(sol,x)/10000;

plot(x + 1981,y)

legend('Measured Data','Model Predicted','Location','SouthEast')
axis([1980 2004 0 5.5])