Added Zangwill optimization method