%U S. Department of Commerce
%Bureau of the Census
%Subject Reports
%Employment Status and Work Experience 1950,1960,1970,1990


%-----------------------------------------------------------------------------
%                     Male Hours of work - 1950
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
mhrs50=[  ...
 14 99   40062870   38945970     767430    1803960    1281180    1088580
 14 17     953760     919620     190680     239010      64980      24510
 18 19    1117830    1087230      59760     116760      52950      35460
 20 24    3893190    3800310      79830     199350     129690     107490
 25 29    4900350    4799400      54060     166560     134520     127980
 30 34    4893210    4799190      39750     136740     122460     117180
 35 44    9387780    9183720      81570     258450     242790     241470
 45 54    7471860    7260870      79920     244470     217170     201270
 55 64    5190570    4983900      86610     237000     189450     154740
 65 74    1942080    1826250      71850     159420     104670      65910
 75 99     312240     285480      23400      46200      22500      12570 ]; 

%-----------------------------------------------------------------------------
%                      At Work (cont )               Hours Not    W/ a Job
% Age       40        41-48      49-59       60+      Reported    Not at work
%-----------------------------------------------------------------------------
tem=[  ....
 14 99   15751590    7329780    4063980    5919390     940080    1116900
 14 17     130890      65010      60990      83250      60300      34140
 18 19     391500     169020     100980     129600      31200      30600
 20 24    1626900     727980     381030     458820      89220      92880
 25 29    2119470     954810     501390     640260     100350     100950
 30 34    2062620     963450     526560     737580      92850      94020
 35 44    3766560    1820520    1030440    1553250     188670     204060
 45 54    2974290    1386630     781950    1206870     168300     210990
 55 64    2017770     888540     484380     792570     132840     206670
 65 74     595320     313440     172320     281190      62130     115830
 75 99      66270      40380      23940      36000      14220      26760 ]; 
mhrs50=[mhrs50,tem(:,3:8)];
%-----------------------------------------------------------------------------



%-----------------------------------------------------------------------------
%                     Female Hours of work - 1950
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
fhrs50=[  ...
 14 99   15753000   15119910     672390    1463100     871500    1101210
 14 17     413970     388920     108960      90960      20820      15600
 18 19     887160     863610      46920      69510      33870      60090
 20 24    2378370    2312220      58830     140310      93000     182850
 25 29    1947060    1877460      54150     146520      94500     138780
 30 34    1737300    1671210      59940     159150     102240     116190
 35 44    3652860    3509640     127200     357330     228720     261660
 45 54    2751720    2632980     103740     272040     174870     202980
 55 64    1494510    1414260      72780     164850      93720      98370
 65 74     436770     402600      33840      55050      27120      22350
 75 99      53280      47010       6030       7380       2640       2340 ];
%-----------------------------------------------------------------------------
%                      At Work (cont )               Hours Not    W/ a Job
% Age       40        41-48      49-59       60+      Reported    Not at work
%-----------------------------------------------------------------------------
tem=[  ...
 14 99    6652150    2518010     671430     691920     478200     633090
 14 17      72330      32010      11940       8730      27570      25050
 18 19     423420     158880      26700      16020      28200      23550
 20 24    1236060     412230      83460      42450      63030      66150
 25 29     961200     312090      69390      51840      48990      69600
 30 34     766140     282450      69420      68130      47550      66090
 35 44    1515750     584640     157530     176850      99960     143220
 45 54    1056450     432240     138630     168210      83820     118740
 55 64     500950     234590      85020     113970      50010      80250
 65 74     112560      63390      25050      40470      22770      34170
 75 99       7290       5490       4290       5250       6300       6270 ]; 
fhrs50=[fhrs50,tem(:,3:8)];
%-----------------------------------------------------------------------------


%-----------------------------------------------------------------------------
%                     Male Hours of work - 1960
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
mhrs60=[  ...
 14 99   43490888   42558076    1882999    1959555    1314864    1905776
 14 17    1317583    1282756     633564     339071      63822      34672
 18 19    1154406    1132428     190283     177412      58455      51602
 20 24    3619470    3560882     194497     218705     123005     162554
 25 29    4520650    4459948      91310     140123     116486     206058
 30 34    5193966    5123329      67250     114604     120362     220706
 35 39    5460782    5376884      68644     115716     123527     214226
 40 44    5090571    5002031      69131     113959     122166     191807
 45 49    4835906    4740311      75314     123636     124660     197414
 50 54    4193914    4094715      80926     125635     120637     196020
 55 59    3495840    3389594      81475     123674     117136     177145
 60 64    2495279    2400249      81181     109361      94470     136642
 65 69    1192720    1130678     113459     122316      65186      64672
 70 74     579470     546348      79532      84308      39296      31899
 75 99     340331     317923      56433      51035      25656      20359 ];
%-----------------------------------------------------------------------------
%                      At Work (cont )               Hours Not    W/ a Job
% Age       40        41-48      49-59       60+      Reported    Not at work
%-----------------------------------------------------------------------------
tem=[   ...
 14 99   17700401    8335928    4407080    5051473        0.0     932812
 14 17     113830      49020      21824      26953        0.0      34827
 18 19     368552     153363      69741      63020        0.0      21978
 20 24    1512442     694176     349069     306434        0.0      58588
 25 29    1959348     933168     515705     497750        0.0      60702
 30 34    2201583    1101360     640284     657180        0.0      70637
 35 39    2324531    1151027     659743     719470        0.0      83898
 40 44    2174593    1055460     587015     687900        0.0      88540
 45 49    2056209     982850     521402     658826        0.0      95595
 50 54    1804135     821018     411555     534789        0.0      99199
 55 59    1501253     656125     312620     420166        0.0     106246
 60 64    1071585     442908     192233     271869        0.0      95030
 65 69     390433     178405      76374     119833        0.0      62042
 70 74     146941      76182      32154      56036        0.0      33122
 75 99      74966      40866      17361      31247        0.0      22408 ];
mhrs60=[mhrs60,tem(:,3:8)];
%-----------------------------------------------------------------------------



%-----------------------------------------------------------------------------
%                     Female Hours of work - 1960
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
fhrs60=[   ...
 14 99   21155675   20454329    2011026    2321741    1324866    2362804
 14 17     692464     662261     356946     149286      22212      23327
 18 19    1016640     992350     150070     103942      40215     114279
 20 24    2296516    2232034     161492     167494     105859     303292
 25 29    1816557    1755907     128493     169965     101906     223877
 30 34    2042376    1971469     161750     229620     130325     229405
 35 39    2453641    2375267     183524     286698     163776     263064
 40 44    2544376    2466481     174634     290359     175573     278754
 45 49    2519392    2441607     167776     270165     177883     288963
 50 54    2166184    2096129     141262     218819     153265     259900
 55 59    1679152    1621553     125364     169303     120215     196867
 60 64    1048379    1007154     101152     124693      72413     113252
 65 69     520961     496458      84077      83185      37278      43692
 70 74     231724     218833      47673      38958      15633      15785
 75 99     127313     116826      26813      19254       8313       8347 ];
%-----------------------------------------------------------------------------
%                      At Work (cont )               Hours Not    W/ a Job
% Age       40        41-48      49-59       60+      Reported    Not at work
%-----------------------------------------------------------------------------
tem=[   ...
 14 99    8744142    2408570     611096     670084        0.0     701346
 14 17      77116      21449       6524       5401        0.0      30203
 18 19     459842      97122      17743       9137        0.0      24290
 20 24    1160181     253251      51457      29008        0.0      64482
 25 29     841909     210531      45742      33484        0.0      60650
 30 34     881725     235314      55670      47660        0.0      70907
 35 39    1051849     286370      71568      68418        0.0      78374
 40 44    1079306     306486      78808      82561        0.0      77895
 45 49    1042076     311824      85163      97757        0.0      77785
 50 54     881304     271651      75359      94569        0.0      70055
 55 59     664885     205271      57238      82410        0.0      57599
 60 64     379304     122291      36408      57641        0.0      41225
 65 69     143917      52839      17177      34293        0.0      24503
 70 74      53198      22163       8069      17354        0.0      12891
 75 99      27530      12008       4170      10391        0.0      10487 ];
fhrs60=[fhrs60,tem(:,3:8)];
%-----------------------------------------------------------------------------


%-----------------------------------------------------------------------------
%                     Male Hours of work - 1970
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
mhrs70=[   ...
 14 99   48272031   46986910    2132474    2681302    2361175    2308900
 14 15     541370     528116     288409      95175      45040       6703
 16 99   47730661   46458794    1844065    2586127    2316135    2302197
 16 17    1240478    1213670     429337     447142      97724      37568
 18 19    1743203    1700447     285533     372171     146031      85867
 20 21    1708320    1662758     174485     217561     117456      85573
 22 24    3287921    3213840     121800     212357     184994     175338
 25 29    5659685    5543890      96384     187855     258936     293007
 30 34    5007620    4903858      55937     124574     210096     238402
 35 39    4878103    4767278      54133     106108     195608     224917
 40 44    5301457    5171604      57244     113886     203387     233744
 45 49    5258023    5115912      61279     118006     204324     229204
 50 54    4728191    4584544      65435     116049     201052     211857
 55 59    4041490    3904634      71363     120026     184423     203641
 60 64    2871717    2760795      95082     133491     151514     167008
 65 99    2004453    1915564     276053     316901     160590     116071 ];
%-----------------------------------------------------------------------------
%                      At Work (cont )                 Median      W/ a Job
% Age       40        41-48      49-59       60+        Hours    Not at work
%-----------------------------------------------------------------------------
tem=[    ...
 14 99   20232926    8182313    4692564    4395256       40.7    1285121
 14 15      78552       6662       3497       4078       13.7      13254
 16 99   20154374    8175651    4689067    4391178       40.7    1271867
 16 17     143569      34317      12744      11269       21.0      26808
 18 19     545469     153822      66955      44599       37.7      42756
 20 21     673215     215811     105957      72700       40.4      45562
 22 24    1477697     551737     286446     203471       40.6      74081
 25 29    2489336    1075151     632956     510265       40.8     115795
 30 34    2131770     971498     613656     557925       40.9     103762
 35 39    2080157     930218     594325     581812       40.9     110825
 40 44    2293136    1021479     629599     619129       40.9     129853
 45 49    2366234     987202     577665     571998       40.8     142111
 50 54    2185981     856207     472761     475202       40.8     143647
 55 59    1867835     704725     371468     381153       40.7     136856
 60 64    1301659     459641     217791     234609       40.6     110922
 65 99     598316     213843     106744     127046       40.1      88889 ];
mhrs70=[mhrs70,tem(:,3:8)];
%-----------------------------------------------------------------------------


%-----------------------------------------------------------------------------
%                     Female Hours of work - 1970
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
fhrs70=[    ...
 14 99   29331612   28181527    2610162    3807364    2448274    3272971
 14 15     257102     246766     133563      32022      22635       3982
 16 99   29074510   27934761    2476599    3775342    2425639    3268989
 16 17     778086     754408     329408     264994      44711      20749
 18 19    1555225    1508180     283880     271134     111240     147674
 20 21    1790649    1730201     191356     196970     124820     200574
 22 24    2584842    2483342     126355     235529     199179     340358
 25 29    2941662    2816870     191042     319600     240565     339489
 30 34    2479047    2377781     196967     324688     215194     257686
 35 39    2652506    2547532     191062     368984     236947     284638
 40 44    3083930    2964996     195289     409321     274541     344391
 45 49    3204356    3080232     175617     389143     278363     376599
 50 54    2897005    2786200     153317     319915     250708     358636
 55 59    2401488    2305874     134190     257950     211725     304792
 60 64    1614058    1546518     122801     201270     142861     200451
 65 99    1091656    1032627     185315     215844      94785      92952 ];
%-----------------------------------------------------------------------------
%                      At Work (cont )                 Median      W/ a Job
% Age       40        41-48      49-59       60+        Hours    Not at work
%-----------------------------------------------------------------------------
tem=[    ...
 14 99   12634691    2169476     649057     589532       40.2    1150085
 14 15      48649       3489       1361       1065       13.9      10336
 16 99   12586042    2165987     647696     588467       40.2    1139749
 16 17      78183      10820       2851       2692       17.7      23678
 18 19     589868      78054      17433       8897       38.0      47045
 20 21     857632     118941      26192      13716       40.2      60448
 22 24    1305274     201043      48793      26811       40.3     101500
 25 29    1390905     229916      63480      41873       40.2     124792
 30 34    1097046     183779      57548      44873       40.2     101266
 35 39    1142568     203710      64178      55445       40.2     104974
 40 44    1351254     242354      76375      71471       40.2     118934
 45 49    1440576     260099      79355      80480       40.2     124124
 50 54    1306812     240752      77689      78371       40.2     110805
 55 59    1054903     205037      65154      72123       40.2      95614
 60 64     663428     125097      41979      48631       40.2      67540
 65 99     307593      66385      26669      43084       36.1      59029 ];
fhrs70=[fhrs70,tem(:,3:8)];
%-----------------------------------------------------------------------------


% 1980 coming.


%-----------------------------------------------------------------------------
%                     Male Hours of work - 1990
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
mhrs90=[   ...
 16 99   62114157   61067187    1939471    4130358    2162755    2200269
 16 17     638947     627975     180412     324142      44987      16564
 18 19    1891774    1849178     287328     602118     172983      95598
 20 21    2359886    2309441     205873     435159     170585     121492
 22 24    4028070    3960815     150071     415948     210388     184388
 25 29    8678620    8550913     159933     412218     293622     318138
 30 34    9364408    9228727     131543     306530     251578     288950
 35 39    8488887    8362897     112764     239332     201284     258079
 40 44    7583438    7473909      96064     184238     164377     235579
 45 49    5950910    5856002      75333     147270     129968     177560
 50 54    4675652    4591533      68768     133055     106460     150159
 55 59    3744413    3664907      72295     155605     110680     136877
 60 64    2607397    2542342     111470     236469     119212     111892
 65 99    2101755    2048548     287617     538274     186631     104993];
%-----------------------------------------------------------------------------
%                      At Work (cont )                 Median      W/ a Job
% Age       40        41-48      49-59       60+        Hours    Not at work
%-----------------------------------------------------------------------------
tem=[    ...
 16 99   26527473    8085846    8864029    7156986        NaN    1046970
 16 17      44900       7066       5073       4831        NaN      10972
 18 19     488047      93704      66636      42764        NaN      42596
 20 21     907901     206548     157900     103983        NaN      50445
 22 24    1772805     499797     432638     294780        NaN      67255
 25 29    3838314    1273478    1278433     976777        NaN     127707
 30 34    4028803    1401319    1571249    1248755        NaN     135681
 35 39    3660200    1237317    1446914    1207007        NaN     125990
 40 44    3378138    1060212    1284315    1070986        NaN     109529
 45 49    2702000     806585     976545     840741        NaN      94908
 50 54    2212046     606660     708761     605624        NaN      84119
 55 59    1803540     470358     503710     411842        NaN      79506
 60 64    1145450     291456     296867     229526        NaN      65055
 65 99     545329     131346     134988     119370        NaN      53207 ];
mhrs90=[mhrs90,tem(:,3:8)];
%-----------------------------------------------------------------------------


%-----------------------------------------------------------------------------
%                     Female Hours of work - 1990
%-----------------------------------------------------------------------------
%          Total      Total                        Hours at Work
% Age    Employed    At Work      1-14      15-29      30-34      35-39
%-----------------------------------------------------------------------------
fhrs90=[    ...
 16 99   52454695   51241869    3299706    7884275    3744988    4768380
 16 17     617218     604488     208676     316389      33379      12892
 18 19    1878051    1827625     381773     681920     174554     118131
 20 21    2262373    2195723     297471     550678     191560     178438
 22 24    3673892    3585156     187496     510008     256683     333070
 25 29    7315161    7121236     278862     798263     460536     639278
 30 34    7632749    7442117     378979     981035     519506     653145
 35 39    7199218    7056535     352410     946983     504182     658833
 40 44    6613130    6494424     270757     793721     461285     631641
 45 49    5083572    4983545     201377     592303     356597     498217
 50 54    3822449    3740175     161834     463658     274058     390989
 55 59    2903815    2835397     150573     401850     220380     319974
 60 64    1966406    1912064     162598     372678     159496     219394
 65 99    1486661    1443384     266900     474789     132772     114378 ];
%-----------------------------------------------------------------------------
%                      At Work (cont )                 Median      W/ a Job
% Age       40        41-48      49-59       60+        Hours    Not at work
%-----------------------------------------------------------------------------
tem=[    ...
 16 99   21701717    4738687    3175874    1928242        NaN    1212826
 16 17      24854       3603       2513       2182        NaN      12730
 18 19     365940      59717      28089      17501        NaN      50426
 20 21     738961     134683      66733      37199        NaN      66650
 22 24    1611778     373232     210603     102286        NaN      88736
 25 29    3326505     823783     519777     274232        NaN     193925
 30 34    3295227     790902     528926     294397        NaN     190632
 35 39    3102711     700312     497187     293917        NaN     142683
 40 44    2921960     648961     472001     294098        NaN     118706
 45 49    2259915     478336     354486     242314        NaN     100027
 50 54    1713022     327327     239236     170051        NaN      82274
 55 59    1266345     221163     147691     107421        NaN      68418
 60 64     744788     122949      72465      57696        NaN      54342
 65 99     329711      53719      36167      34948        NaN      43277 ];
fhrs90=[fhrs90,tem(:,3:8)];
%-----------------------------------------------------------------------------

%-----------------------------------------------------------------------------
%           1950 Population
% Age      Male       Female
%-----------------------------------------------------------------------------
pop50 = [   ...
 14 99   54610050   57083190 
 14 19    6443040    6369330 
 14 14    1098030    1042320 
 15 15    1078800    1049460 
 16 16    1057290    1025940 
 17 17    1059720    1044600 
 18 18    1093050    1098240 
 19 19    1056150    1108770 
 20 24    5540520    5869860 
 25 29    5900970    6283110 
 30 34    5553210    5874930 
 35 39    5453430    5720010 
 40 44    4945290    5127480 
 45 49    4463310    4562760 
 50 54    4032810    4125060 
 55 59    3577470    3614720 
 60 64    2969790    3021430 
 65 69    2396430    2605600 
 70 74    1599450    1790120 
 75 99    1734330    2118780 ];
%-----------------------------------------------------------------------------

%-----------------------------------------------------------------------------
%           1960 Population
% Age      Male       Female
%-----------------------------------------------------------------------------
pop60=[   ...
 14 99   61353999   64912621 
 14 17    5729892    5523939 
 14 14    1396738    1350403 
 15 15    1430619    1361949 
 16 16    1449550    1396869 
 17 17    1452985    1414718 
 18 19    2350279    2414634 
 20 24    5276385    5506450 
 25 29    5338872    5537408 
 30 34    5844697    6099867 
 35 39    6094607    6438899 
 40 44    5663293    5897534 
 45 49    5377304    5559569 
 50 54    4783543    4927163 
 55 59    4194542    4405945 
 60 64    3385060    3718904 
 65 69    2883723    3294987 
 70 74    2136352    2528914 
 75 99    2295450    3058408 ];
%-----------------------------------------------------------------------------


%-----------------------------------------------------------------------------
%           1970 Population
% Age      Male       Female
%-----------------------------------------------------------------------------
pop70=[   ...
 14 99   71485878   77910094 
 14 15    4246149    4055480 
 16 99   67239729   73854614 
 16 21   10797149   10956387 
 16 16    2005360    1930161 
 17 17    1956794    1888548 
 18 18    1899720    1875097 
 19 19    1746984    1784999 
 20 20    1637976    1771502 
 21 21    1550315    1706080 
 22 24    4572918    4876927 
 25 29    6569934    6810076 
 30 34    5607593    5868858 
 35 39    5431809    5710623 
 40 44    5829922    6149692 
 45 49    5808783    6255178 
 50 54    5329398    5741230 
 55 59    4800385    5228025 
 60 64    4058508    4599123 
 65 69    3116348    3897364 
 70 74    2324327    3115737 
 75 99    2992655    4645394 ];
%-----------------------------------------------------------------------------


%-----------------------------------------------------------------------------
%           1990 Population
% Age      Male       Female
%-----------------------------------------------------------------------------
pop90=[   ...
 16 99   90061061   97878799
 16 21    9502802    9091580
 22 24    5531940    5401073
 25 29   10649679   10579184
 30 34   11020469   11098529
 35 39    9847801   10032992
 40 44    8707322    8922533
 45 49    6844044    7110521
 50 54    5550784    5896318
 55 59    4979082    5457766
 60 64    4950952    5660279
 65 69    4536086    5591262
 70 74    3427782    4607635
 75 99    4512318    8429127];
%-----------------------------------------------------------------------------


%Statistics for 1950
%-------------------

% Average hours across categories:

avgh= [7.5,22,32,37,40,44.5,54,62.5];

% Find the fraction of 14 year olds of the 14-17 group in 1950:

frac14 = pop50(3,3:4)./sum(pop50(3:6,3:4));

% Calculate hours worked for 15-17 in 1950: 

tem1   = (1-frac14(1))*mhrs50(2,:);
tem2   = (1-frac14(2))*fhrs50(2,:);

% Calculate employment and hours worked for 15-24,25-34, etc  in 1950: 

Mhrs50 = [sum([tem1;mhrs50(3:4,:)]);
          sum(mhrs50(5:6,:));
          mhrs50(7:11,:)];
Fhrs50 = [sum([tem2;fhrs50(3:4,:)]);
          sum(fhrs50(5:6,:));
          fhrs50(7:11,:)];
label  = [15 24;25 34;35 44;45 54;55 64;65 74;75 99];
Pop50  = [sum(pop50(4:9,3:4));
          sum(pop50(10:11,3:4));
          sum(pop50(12:13,3:4));
          sum(pop50(14:15,3:4));
          sum(pop50(16:17,3:4));
          sum(pop50(18:19,3:4));
          pop50(20,3:4)];
Memp50 = Mhrs50(:,3)./Pop50(:,1);
Femp50 = Fhrs50(:,3)./Pop50(:,2);
Mhrs50 = Mhrs50(:,5:12)./(sum(Mhrs50(:,5:12)')'*ones(1,8));
Fhrs50 = Fhrs50(:,5:12)./(sum(Fhrs50(:,5:12)')'*ones(1,8));
Mavg50 = Mhrs50*avgh';
Favg50 = Fhrs50*avgh';
clc
home
disp('1950')
disp('----')
disp('                Average              Employment/')
disp('                 Hours               Population')
disp('            --------------------------------------')
disp('  Ages       Male     Female       Male     Female') 
disp('--------------------------------------------------')
disp(sprintf(' %g %g %10.2f %10.2f %10.2f %10.2f\n', ...
             [label,Mavg50,Favg50,Memp50,Femp50]'))


%Statistics for 1960
%-------------------

% Find the fraction of 14 year olds of the 14-17 group in 1960:

frac14 = pop60(3,3:4)./sum(pop60(3:6,3:4));

% Calculate hours worked for 15-17 in 1960: 

tem1   = (1-frac14(1))*mhrs60(2,:);
tem2   = (1-frac14(2))*fhrs60(2,:);

% Calculate employment and hours worked for 15-24,25-34, etc  in 1960: 

Mhrs60 = [sum([tem1;mhrs60(3:4,:)]);
          sum(mhrs60(5:6,:));
          sum(mhrs60(7:8,:));
          sum(mhrs60(9:10,:));
          sum(mhrs60(11:12,:));
          sum(mhrs60(13:14,:));
          mhrs60(15,:)];
Fhrs60 = [sum([tem2;fhrs60(3:4,:)]);
          sum(fhrs60(5:6,:));
          sum(fhrs60(7:8,:));
          sum(fhrs60(9:10,:));
          sum(fhrs60(11:12,:));
          sum(fhrs60(13:14,:));
          fhrs60(15,:)];
label  = [15 24;25 34;35 44;45 54;55 64;65 74;75 99];
Pop60  = [sum(pop60(4:8,3:4));
          sum(pop60(9:10,3:4));
          sum(pop60(11:12,3:4));
          sum(pop60(13:14,3:4));
          sum(pop60(15:16,3:4));
          sum(pop60(17:18,3:4));
          pop60(19,3:4)];
Memp60 = Mhrs60(:,3)./Pop60(:,1);
Femp60 = Fhrs60(:,3)./Pop60(:,2);
Mhrs60 = Mhrs60(:,5:12)./(sum(Mhrs60(:,5:12)')'*ones(1,8));
Fhrs60 = Fhrs60(:,5:12)./(sum(Fhrs60(:,5:12)')'*ones(1,8));
Mavg60 = Mhrs60*avgh';
Favg60 = Fhrs60*avgh';
disp(' ')
disp(' ')
disp('1960')
disp('----')
disp('                Average              Employment/')
disp('                 Hours               Population')
disp('            --------------------------------------')
disp('  Ages       Male     Female       Male     Female') 
disp('--------------------------------------------------')
disp(sprintf(' %g %g %10.2f %10.2f %10.2f %10.2f\n', ...
             [label,Mavg60,Favg60,Memp60,Femp60]'))



%Statistics for 1970
%-------------------

% Assume the fraction of 14 year olds of the 14-15 group in 1970 is 0.5.
% Calculate employment and hours worked for 15-24,25-34, etc  in 1970: 

tem1   = 0.5*mhrs70(2,:);
tem2   = 0.5*fhrs70(2,:);

Mhrs70 = [sum([tem1;mhrs70(4:7,:)]);
          sum(mhrs70(8:9,:));
          sum(mhrs70(10:11,:));
          sum(mhrs70(12:13,:));
          sum(mhrs70(14:15,:));
          mhrs70(16,:)];
Fhrs70 = [sum([tem2;fhrs70(4:7,:)]);
          sum(fhrs70(8:9,:));
          sum(fhrs70(10:11,:));
          sum(fhrs70(12:13,:));
          sum(fhrs70(14:15,:));
          fhrs70(16,:)];
label  = [15 24;25 34;35 44;45 54;55 64;65 99];
Pop70  = [round(0.5*pop70(2,3:4))+sum(pop70(5:11,3:4));
          sum(pop70(12:13,3:4));
          sum(pop70(14:15,3:4));
          sum(pop70(16:17,3:4));
          sum(pop70(18:19,3:4));
          sum(pop70(20:22,3:4))];
Memp70 = Mhrs70(:,3)./Pop70(:,1);
Femp70 = Fhrs70(:,3)./Pop70(:,2);
Mhrs70 = Mhrs70(:,5:12)./(sum(Mhrs70(:,5:12)')'*ones(1,8));
Fhrs70 = Fhrs70(:,5:12)./(sum(Fhrs70(:,5:12)')'*ones(1,8));
Mavg70 = Mhrs70*avgh';
Favg70 = Fhrs70*avgh';

disp(' ')
disp(' ')
disp('1970')
disp('----')
disp('                Average              Employment/')
disp('                 Hours               Population')
disp('            --------------------------------------')
disp('  Ages       Male     Female       Male     Female') 
disp('--------------------------------------------------')
disp(sprintf(' %g %g %10.2f %10.2f %10.2f %10.2f\n', ...
             [label,Mavg70,Favg70,Memp70,Femp70]'))


%Statistics for 1990
%-------------------

% Assume the population of 15 year olds is equal to one-sixth of
% the population of 16 - 21 year olds.  Assume that 15 year olds
% don't work. Calculate employment and hours worked for 15-24,25-34, etc  
% in 1990: 

Mhrs90 = [sum(mhrs90(2:5,:));
          sum(mhrs90(6:7,:));
          sum(mhrs90(8:9,:));
          sum(mhrs90(10:11,:));
          sum(mhrs90(12:13,:));
          mhrs90(14,:)];
Fhrs90 = [sum(fhrs90(2:5,:));
          sum(fhrs90(6:7,:));
          sum(fhrs90(8:9,:));
          sum(fhrs90(10:11,:));
          sum(fhrs90(12:13,:));
          fhrs90(14,:)];
label  = [15 24;25 34;35 44;45 54;55 64;65 99];
Pop90  = [round(pop90(2,3:4)/6)+sum(pop90(2:3,3:4));
          sum(pop90(4:5,3:4));
          sum(pop90(6:7,3:4));
          sum(pop90(8:9,3:4));
          sum(pop90(10:11,3:4));
          sum(pop90(12:14,3:4))];
Memp90 = Mhrs90(:,3)./Pop90(:,1);
Femp90 = Fhrs90(:,3)./Pop90(:,2);
Mhrs90 = Mhrs90(:,5:12)./(sum(Mhrs90(:,5:12)')'*ones(1,8));
Fhrs90 = Fhrs90(:,5:12)./(sum(Fhrs90(:,5:12)')'*ones(1,8));
Mavg90 = Mhrs90*avgh';
Favg90 = Fhrs90*avgh';

disp(' ')
disp(' ')
disp('1990')
disp('----')
disp('                Average              Employment/')
disp('                 Hours               Population')
disp('            --------------------------------------')
disp('  Ages       Male     Female       Male     Female') 
disp('--------------------------------------------------')
disp(sprintf(' %g %g %10.2f %10.2f %10.2f %10.2f\n', ...
             [label,Mavg90,Favg90,Memp90,Femp90]'))


disp(' ')
disp(' ')
disp('Average Hours of Males')
disp('----------------------')
disp(' ')
disp('                             Average Hours by Age               ')
disp('            ----------------------------------------------------')
disp('Age in 50    15-24    25-34    35-44    45-54    55-64    65-74')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Mavg50(1),Mavg60(2),Mavg70(3),  Mavg90(5)]'))
disp(sprintf(' %g %g             %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Mavg50(2),Mavg60(3),Mavg70(4),      Mavg90(6)]'))
disp(sprintf(' %g %g                      %8.2f %8.2f %8.2f  ', ...
             [35,44,Mavg50(3),Mavg60(4),Mavg70(5)      ]'))
disp(sprintf(' %g %g                               %8.2f %8.2f %8.2f ', ...
             [45,54,Mavg50(4),Mavg60(5),Mavg70(6)]'))
disp(sprintf(' %g %g                                        %8.2f %8.2f', ...
             [54,64,Mavg50(5),Mavg60(6)]'))
disp(sprintf(' %g %g                                                 %8.2f', ...
             [65,74,Mavg50(6)]'))

disp(' ')
disp(' ')
disp('                        Average Hours by Census Year            ')
disp('            ----------------------------------------------------')
disp('Age in 50     1950     1960     1970     1980     1990          ')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Mavg50(1),Mavg60(2),Mavg70(3),  Mavg90(5)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Mavg50(2),Mavg60(3),Mavg70(4),      Mavg90(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f  ', ...
             [35,44,Mavg50(3),Mavg60(4),Mavg70(5)      ]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f ', ...
             [45,54,Mavg50(4),Mavg60(5),Mavg70(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f', ...
             [54,64,Mavg50(5),Mavg60(6)]'))
disp(sprintf(' %g %g    %8.2f', ...
             [65,74,Mavg50(6)]'))


disp(' ')
disp(' ')
disp('Employment/Population of Males')
disp('------------------------------')
disp(' ')
disp('                       Employment/Population by Age             ')
disp('            ----------------------------------------------------')
disp('Age in 50    15-24    25-34    35-44    45-54    55-64    65-74')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Memp50(1),Memp60(2),Memp70(3),  Memp90(5)]'))
disp(sprintf(' %g %g             %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Memp50(2),Memp60(3),Memp70(4),      Memp90(6)]'))
disp(sprintf(' %g %g                      %8.2f %8.2f %8.2f  ', ...
             [35,44,Memp50(3),Memp60(4),Memp70(5)      ]'))
disp(sprintf(' %g %g                               %8.2f %8.2f %8.2f ', ...
             [45,54,Memp50(4),Memp60(5),Memp70(6)]'))
disp(sprintf(' %g %g                                        %8.2f %8.2f', ...
             [54,64,Memp50(5),Memp60(6)]'))
disp(sprintf(' %g %g                                                 %8.2f', ...
             [65,74,Memp50(6)]'))

disp(' ')
disp(' ')
disp('                Employment/Population by Census Year            ')
disp('            ----------------------------------------------------')
disp('Age in 50     1950     1960     1970     1980     1990          ')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Memp50(1),Memp60(2),Memp70(3),  Memp90(5)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Memp50(2),Memp60(3),Memp70(4),      Memp90(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f  ', ...
             [35,44,Memp50(3),Memp60(4),Memp70(5)      ]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f ', ...
             [45,54,Memp50(4),Memp60(5),Memp70(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f', ...
             [54,64,Memp50(5),Memp60(6)]'))
disp(sprintf(' %g %g    %8.2f', ...
             [65,74,Memp50(6)]'))



disp(' ')
disp(' ')
disp('Average Hours of Females')
disp('------------------------')
disp(' ')
disp('                             Average Hours by Age               ')
disp('            ----------------------------------------------------')
disp('Age in 50    15-24    25-34    35-44    45-54    55-64    65-74')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Favg50(1),Favg60(2),Favg70(3),  Favg90(5)]'))
disp(sprintf(' %g %g             %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Favg50(2),Favg60(3),Favg70(4),      Favg90(6)]'))
disp(sprintf(' %g %g                      %8.2f %8.2f %8.2f  ', ...
             [35,44,Favg50(3),Favg60(4),Favg70(5)      ]'))
disp(sprintf(' %g %g                               %8.2f %8.2f %8.2f ', ...
             [45,54,Favg50(4),Favg60(5),Favg70(6)]'))
disp(sprintf(' %g %g                                        %8.2f %8.2f', ...
             [54,64,Favg50(5),Favg60(6)]'))
disp(sprintf(' %g %g                                                 %8.2f', ...
             [65,74,Favg50(6)]'))

disp(' ')
disp(' ')
disp('                        Average Hours by Census Year            ')
disp('            ----------------------------------------------------')
disp('Age in 50     1950     1960     1970     1980     1990          ')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Favg50(1),Favg60(2),Favg70(3),  Favg90(5)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Favg50(2),Favg60(3),Favg70(4),      Favg90(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f  ', ...
             [35,44,Favg50(3),Favg60(4),Favg70(5)      ]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f ', ...
             [45,54,Favg50(4),Favg60(5),Favg70(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f', ...
             [54,64,Favg50(5),Favg60(6)]'))
disp(sprintf(' %g %g    %8.2f', ...
             [65,74,Favg50(6)]'))


disp(' ')
disp(' ')
disp('Employment/Population of Females')
disp('--------------------------------')
disp(' ')
disp('                       Employment/Population by Age             ')
disp('            ----------------------------------------------------')
disp('Age in 50    15-24    25-34    35-44    45-54    55-64    65-74')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Femp50(1),Femp60(2),Femp70(3),  Femp90(5)]'))
disp(sprintf(' %g %g             %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Femp50(2),Femp60(3),Femp70(4),      Femp90(6)]'))
disp(sprintf(' %g %g                      %8.2f %8.2f %8.2f  ', ...
             [35,44,Femp50(3),Femp60(4),Femp70(5)      ]'))
disp(sprintf(' %g %g                               %8.2f %8.2f %8.2f ', ...
             [45,54,Femp50(4),Femp60(5),Femp70(6)]'))
disp(sprintf(' %g %g                                        %8.2f %8.2f', ...
             [54,64,Femp50(5),Femp60(6)]'))
disp(sprintf(' %g %g                                                 %8.2f', ...
             [65,74,Femp50(6)]'))

disp(' ')
disp(' ')
disp('                Employment/Population by Census Year            ')
disp('            ----------------------------------------------------')
disp('Age in 50     1950     1960     1970     1980     1990          ')
disp('----------------------------------------------------------------')
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [15,24,Femp50(1),Femp60(2),Femp70(3),  Femp90(5)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f          %8.2f ', ...
             [25,34,Femp50(2),Femp60(3),Femp70(4),      Femp90(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f  ', ...
             [35,44,Femp50(3),Femp60(4),Femp70(5)      ]'))
disp(sprintf(' %g %g    %8.2f %8.2f %8.2f ', ...
             [45,54,Femp50(4),Femp60(5),Femp70(6)]'))
disp(sprintf(' %g %g    %8.2f %8.2f', ...
             [54,64,Femp50(5),Femp60(6)]'))
disp(sprintf(' %g %g    %8.2f', ...
             [65,74,Femp50(6)]'))
