import nigelzGLPK

ia = nigelzGLPK.new_intArray(10);
ja = nigelzGLPK.new_intArray(10);
ar = nigelzGLPK.new_doubleArray(10);
lp = nigelzGLPK.glp_create_prob();

nigelzGLPK.glp_set_prob_name(lp, "sample");
nigelzGLPK.glp_set_obj_dir(lp, nigelzGLPK.GLP_MAX);
nigelzGLPK.glp_add_rows(lp, 3);
nigelzGLPK.glp_set_row_name(lp, 1, "p");
nigelzGLPK.glp_set_row_bnds(lp, 1, nigelzGLPK.GLP_UP, 0.0, 100.0);
nigelzGLPK.glp_set_row_name(lp, 2, "q");
nigelzGLPK.glp_set_row_bnds(lp, 2, nigelzGLPK.GLP_UP, 0.0, 600.0);
nigelzGLPK.glp_set_row_name(lp, 3, "r");
nigelzGLPK.glp_set_row_bnds(lp, 3, nigelzGLPK.GLP_UP, 0.0, 300.0);
nigelzGLPK.glp_add_cols(lp, 3);
nigelzGLPK.glp_set_col_name(lp, 1, "x1");
nigelzGLPK.glp_set_col_bnds(lp, 1, nigelzGLPK.GLP_LO, 0.0, 0.0);
nigelzGLPK.glp_set_obj_coef(lp, 1, 10.0);
nigelzGLPK.glp_set_col_name(lp, 2, "x2");
nigelzGLPK.glp_set_col_bnds(lp, 2, nigelzGLPK.GLP_LO, 0.0, 0.0);
nigelzGLPK.glp_set_obj_coef(lp, 2, 6.0);
nigelzGLPK.glp_set_col_name(lp, 3, "x3");
nigelzGLPK.glp_set_col_bnds(lp, 3, nigelzGLPK.GLP_LO, 0.0, 0.0);
nigelzGLPK.glp_set_obj_coef(lp, 3, 4.0);

nigelzGLPK.intArray_setitem(ia, 1, 1);
nigelzGLPK.intArray_setitem(ja, 1, 1);
nigelzGLPK.doubleArray_setitem(ar, 1, 1.0); #/* a[1,1] =  1 */

nigelzGLPK.intArray_setitem(ia, 2, 1);
nigelzGLPK.intArray_setitem(ja, 2, 2);
nigelzGLPK.doubleArray_setitem(ar, 2, 1.0); #/* a[1,2] =  1 */

nigelzGLPK.intArray_setitem(ia, 3, 1);
nigelzGLPK.intArray_setitem(ja, 3, 3);
nigelzGLPK.doubleArray_setitem(ar, 3, 1.0); #/* a[1,3] =  1 */

nigelzGLPK.intArray_setitem(ia, 4, 2);
nigelzGLPK.intArray_setitem(ja, 4, 1);
nigelzGLPK.doubleArray_setitem(ar, 4, 10.0); #/* a[2,1] =  1 */

nigelzGLPK.intArray_setitem(ia, 5, 3);
nigelzGLPK.intArray_setitem(ja, 5, 1);
nigelzGLPK.doubleArray_setitem(ar, 5, 2.0); #/* a[3,1] =  2 */

nigelzGLPK.intArray_setitem(ia, 6, 2);
nigelzGLPK.intArray_setitem(ja, 6, 2);
nigelzGLPK.doubleArray_setitem(ar, 6, 4.0); #/* a[2,2] =  4 */

nigelzGLPK.intArray_setitem(ia, 7, 3);
nigelzGLPK.intArray_setitem(ja, 7, 2);
nigelzGLPK.doubleArray_setitem(ar, 7, 2.0); #/* a[3,2] =  2 */

nigelzGLPK.intArray_setitem(ia, 8, 2);
nigelzGLPK.intArray_setitem(ja, 8, 3);
nigelzGLPK.doubleArray_setitem(ar, 8, 5.0); #/* a[2,3] =  5 */

nigelzGLPK.intArray_setitem(ia, 9, 3);
nigelzGLPK.intArray_setitem(ja, 9, 3);
nigelzGLPK.doubleArray_setitem(ar, 9, 6.0); #/* a[3,3] =  6 */

nigelzGLPK.glp_load_matrix(lp, 9, ia, ja, ar);
nigelzGLPK.glp_simplex(lp,None);

Z  = nigelzGLPK.glp_get_obj_val(lp);
x1 = nigelzGLPK.glp_get_col_prim(lp, 1);
x2 = nigelzGLPK.glp_get_col_prim(lp, 2);
x3 = nigelzGLPK.glp_get_col_prim(lp, 3);
print "\nZ = ",Z,"; x1 = ",x1,"; x2 = ",x2,"; x3 = ",x3,"\n"
nigelzGLPK.delete_intArray(ia);
nigelzGLPK.delete_intArray(ja);
nigelzGLPK.delete_doubleArray(ar);
nigelzGLPK.glp_delete_prob(lp);