optimization with pyomo