One-step inverse methods based on deformation theory causes some amount of error. The amount of error is generally increased as the deformation path is more complex. As a remedy, a new three dimensional multi-step inverse method is introduced for optimum design of blank shapes and strain distributions from desired final shapes. The approach extends a one-step inverse method to a multi-step inverse method in order to reduce the amount of error. The algorithm developed is applied to square cup drawing to confirm its validity by demonstrating reasonably accurate numerical results.