Optimal taxation with multiple dimensions of heterogeneity

•We characterize optimal income tax schedules with multidimensional heterogeneity.•We account for the possibility that agents have multiple optimal incomes.•We quantify jumping effects that occur when individuals have multiple optima.•We provide a partial characterization of when jumping effects will occur.•We devise a novel simulation method to solve for optimal tax schedules.