clear

%x=.1:10:2000;
x = logspace(-1,5);

C1=24./x+6./(1+sqrt(x))+0.4;           % White (1991) <2x10e5
C2=24./x+4./x.^(1/3);                  % Putnam (1961) <1000
C3=24./x.*(1+0.15.*x.^.687);           % Schiller-Naumann (1933) <1000
C4=12*x.^(-.5);                        % Edwards et al. (2000)
C5=C3+.42./(1+42500./x.^1.16);         % Clift-Gavin (1970) <2x10e5
C6=24./x+2.6.*x/5./(1+(x/5).^1.52) ... % Morrison (2010)
    +.411.*(x/263000).^(-7.94)./(1+(x/263000).^(-8))+x.^.8/461000;

figure(1)
plot(x,C1,'-k')
hold on
plot(x,C2,'.k')
plot(x,C3,'--k')
plot(x,C4,'-.k')
plot(x,C5,':k')
plot(x,C6,'^k','MarkerSize',6)
hold off
axis([0,10000,0,1])
%axis([00,1000,0,1])

%hleg1 = legend('White (1991)','Clift-Gavin (1970)','Morrison (2010)');
hleg1=legend('White (1991)','Putnam (1961)','Schiller-Naumann (1933)', ... 
    'Edwards et al. (2000)','Clift-Gavin (1970)','Morrison (2010)');

% Results for x=10 to 1000
% C2 agrees with C3
% C4 close to C3
% C5 closer to C3 than C4
% C6 closer to C3 than C4

figure(2)
loglog(x,C1,'-k')
hold on
loglog(x,C2,'.k')
loglog(x,C3,'--k')
loglog(x,C4,'-.k')
loglog(x,C5,':k')
loglog(x,C6,'^k','MarkerSize',6)
%grid on
hold off
set(gca,'YTick',[.1,1,10,100,1000])
axis([.1,100000,.1,1000])
%hleg2 = legend('White (1991)','Clift-Gavin (1970)','Morrison (2010)');
hleg2=legend('White (1991)','Putnam (1961)','Schiller-Naumann (1933)', ... 
    'Edwards et al. (2000)','Clift-Gavin (1970)','Morrison (2010)');