original taxcalc code in SAS from the NBER