We analyze the optimal nonlinear income tax schedule when taxpayers earn multiple incomes and differ along many unobserved dimensions. We derive the necessary conditions for the government’s optimum using both a tax perturbation and a mechanism design approach, and show that both methods produce the same results. Our main contribution is to propose a numerical method to find the optimal tax schedule. Applied to the optimal taxation of couples, we find that optimal isotax curves are very close to linear and parallel. The slope of isotax curves is strongly affected by the relative tax-elasticity of male and female income. We make several additional contributions, including a test for Pareto efficiency and a condition on primitives that ensures the government’s necessary conditions are sufficient and the solution to the problem is unique.