% Figure 1:
% Real compensation per hour of the business sector (1982=100),
% Table B-47, Economic Report of the President 1994

rcomp = [49.4;68.8;91.3;99.5;103.8];
ahrs  = [22.03;20.97;20.55;22;23.62];
t     = 1950:10:1990;

legend off
figure(1)
plot(t,rcomp), xlabel('year'),ylabel('real compensation (1982=100)'), ...
  set(gca,'XTick',t,'XTickLabels',t),grid, ...
print -deps fig1
figure(2)
plot(t,ahrs),  xlabel('year'), ylabel('average hours') ,  ...
  set(gca,'XTick',t,'XTickLabels',t), ...
  axis([1950,1990,0,50]), grid
print -deps fig2


lct= [ ...
   19.13   32.86   36.65   36.27   17.84    0.69    0.00   20.49
   19.64   30.83   34.49   34.24   18.56    1.53    0.00   19.90
   14.05   28.80   32.62   32.75   19.09    2.44    0.00   18.54
   14.15   26.16   29.89   31.47   20.14    3.11    0.00   17.85
   17.47   24.73   28.03   28.16   20.75    4.06    0.05   17.61
   16.75   24.92   27.00   28.27   20.68    5.15    0.58   17.62
   17.65   23.63   27.09   27.63   23.28    5.11    1.18   17.94
   18.81   22.78   26.16   26.31   22.58    6.91    1.39   17.85
   19.26   22.13   25.53   25.69   22.19    8.43    2.17   17.92
   20.09   21.20   25.04   24.83   22.69   12.03    2.97   18.41
   20.91   20.40   24.36   23.90   22.50   14.34    3.93   18.62];
lct=lct(:,1:7);

lcm= [ ...
   20.23   35.27   37.23   36.57   16.95    0.00    0.00   20.89
   21.59   36.00   37.87   36.73   19.39    0.33    0.00   21.70
   15.96   36.15   38.40   37.28   21.64    2.07    0.00   21.64
   17.65   37.89   38.59   37.75   24.44    3.48    0.00   22.83
   22.65   38.20   39.79   37.20   26.73    5.27    0.39   24.32
   23.00   38.60   39.98   38.95   28.38    7.28    1.17   25.34
   25.35   39.42   41.41   39.15   32.84    7.71    2.16   26.86
   27.92   39.93   42.35   40.06   33.60   10.65    2.57   28.15
   29.67   40.50   43.20   41.01   35.34   13.31    3.53   29.51
   31.90   41.13   44.28   41.70   37.66   20.75    5.12   31.79
   34.08   41.71   45.23   42.55   39.27   25.10    7.46   33.63];
lcm=lcm(:,1:7);

lcf= [ ...
   17.99   30.27   35.74   35.87   18.95    2.59    0.00   20.20
   17.68   25.67   31.05   31.85   18.06    2.75    0.00   18.15
   12.18   21.63   26.96   28.48   17.03    2.93    0.00   15.60
   10.73   15.03   21.53   25.48   16.45    2.98    0.00   13.17
   12.43   11.87   16.97   19.78   15.41    3.23    0.14   11.40
   10.68   11.84   14.70   18.38   13.91    3.48    0.43   10.49
   10.18    8.58   13.17   16.48   14.71    3.13    0.66    9.56
   10.00    6.43   10.38   12.58   12.30    3.98    0.74    8.06
    9.19    4.61    8.19   10.18    8.85    4.23    1.29    6.65
    8.69    2.20    6.02    7.45    7.16    3.87    1.36    5.25
    8.19    0.08    3.63    4.44    4.65    4.12    1.01    3.73];
lcf=lcf(:,1:7);

a=15:10:75;
figure(1)
plot(a,lcm), xlabel('age'), ylabel('Average hours')
text(a-.5,lcm(1,:),'10')
text(a-.5,lcm(2,:),'9')
text(a-.5,lcm(3,:),'8')
text(a-.5,lcm(4,:),'7')
text(a-.5,lcm(5,:),'6')
text(a-.5,lcm(6,:),'5')
text(a-.5,lcm(7,:),'4')
text(a-.5,lcm(8,:),'3')
text(a-.5,lcm(9,:),'2')
text(a-.5,lcm(10,:),'1')
text(a-.5,lcm(11,:),'0')
axis([10,80,0,50])
grid
print -deps fig3

figure(2)
plot(a,lcf), xlabel('age'), ylabel('Average hours')
text(a-.5,lcf(1,:),'10')
text(a-.5,lcf(2,:),'9')
text(a-.5,lcf(3,:),'8')
text(a-.5,lcf(4,:),'7')
text(a-.5,lcf(5,:),'6')
text(a-.5,lcf(6,:),'5')
text(a-.5,lcf(7,:),'4')
text(a-.5,lcf(8,:),'3')
text(a-.5,lcf(9,:),'2')
text(a-.5,lcf(10,:),'1')
text(a-.5,lcf(11,:),'0')
axis([10,80,0,50])
grid
print -deps fig4

figure(3)
plot(a,lct), xlabel('age'), ylabel('Average hours')
text(a-.5,lct(1,:),'10')
text(a-.5,lct(2,:),'9')
text(a-.5,lct(3,:),'8')
text(a-.5,lct(4,:),'7')
text(a-.5,lct(5,:),'6')
text(a-.5,lct(6,:),'5')
text(a-.5,lct(7,:),'4')
text(a-.5,lct(8,:),'3')
text(a-.5,lct(9,:),'2')
text(a-.5,lct(10,:),'1')
text(a-.5,lct(11,:),'0')
axis([10,80,0,50])
grid
print -deps fig5

