%PLOTGAINS
%
 
  load trantfWPHG.dat
  curpolt = trantfWPHG;
  load trantfWPHGhighercopay.dat;
  newpolt = trantfWPHGhighercopay;
[WG,Ujf1,Ujc1] = gains(curpolt,newpolt);

figure,plot(WG,'LineWidth',2)
line([0,200],[0,0],'Color','Black')
xlabel('Birth-Year Cohort','Fontsize',14)
axis([0,200,-8,2])
set(gca,'XTickLabel',[100,50,0,-50,-100],'Fontsize',14)
legend('Low','Medium','High','Top 1%','Location','Southwest')
legend('boxoff')
title({''; ...
 '% Welfare changes (VAT vs Higher Copay)'}, ...
 'Fontsize',14)

