load Zwick.dat
x           = Zwick'/10^9;
sz14        = x(2:14,[4,2]);
szz14       = x(15:27,[4,2]);
ii          = [1,3,6,8,10,12];
jj          = [2,4,7,9,11,13];
szz14w      = szz14(ii,:);
sz14w       = sz14(ii,:);
szz14y      = szz14(jj,:);
sz14y       = sz14(jj,:);
szz14y(2,:) = szz14y(2,:)+szz14(5,:);
sz14y(2,:)  = sz14y(2,:)+sz14(5,:);
gdp14       = 17527;

sz77        = x(2:14,[3,1]);
szz77       = x(15:27,[3,1]);
ii          = [1,3,6,8,10,12];
jj          = [2,4,7,9,11,13];
szz77w      = szz77(ii,:);
sz77w       = sz77(ii,:);
szz77y      = szz77(jj,:);
sz77y       = sz77(jj,:);
szz77y(2,:) = szz77y(2,:)+szz77(5,:);
sz77y(2,:)  = sz77y(2,:)+sz77(5,:);
gdp77       = 2081;

disp('1977 Income (non-housing)')
disp('=========================')

disp('  SZ        SZZ ')
disp('-----------------')
disp('Top All   Top All')
disp([sz77y/gdp77 szz77y/gdp77])
disp('-----------------')
disp(sum([sz77y(1:5,:)/gdp77 szz77y(1:5,:)/gdp77]))

disp(' ')

disp('1977 Wealth')
disp('===========')

disp('  SZ        SZZ ')
disp('-----------------')
disp('Top All   Top All')
disp([sz77w/gdp77 szz77w/gdp77])
disp('-----------------')
disp(sum([sz77w/gdp77 szz77w/gdp77]))

disp(' ')

disp('2014 Income')
disp('===========')

disp('  SZ        SZZ ')
disp('-----------------')
disp('Top All   Top All')
disp([sz14y/gdp14 szz14y/gdp14])
disp('-----------------')
disp(sum([sz14y/gdp14 szz14y/gdp14]))

disp(' ')

disp('2014 Wealth')
disp('===========')

disp('  SZ        SZZ ')
disp('-----------------')
disp('Top All   Top All')
disp([sz14w/gdp14 szz14w/gdp14])
disp('-----------------')
disp(sum([sz14w/gdp14 szz14w/gdp14]))
