Replicates Saez 2001 Mirrlees analysis using Matlab and AMPL.