Below is the syntax highlighted version of LPDemo.java
from §9.6 Numerical Linear Algebra.
/****************************************************************************** * Compilation: javac -cp .:or124.jar LPDemo.java * Exception: java -cp .:or124.jar LPDemo * Dependencies: or124.jar * * max 13x1 + 23x2 * 5x1 + 15x2 <= 480 * 4x1 + 4x2 <= 160 * 35x1 + 20x2 <= 1190 * x1 , x2 >= 0 * * Copied from http://opsresearch.com/OR-Objects/faq.html * % java LPDemo * 800.0 * Vector[2] = * [0] 11.999999999999998 * [1] 28.0 * * The optimal value is 800 with x1 = 12 and x2 = 28. * ******************************************************************************/ import drasys.or.mp.Problem; import drasys.or.mp.lp.DenseSimplex; public class LPDemo { public static void main(String[] args) throws Exception { // final byte LTE = Constraint.LESS; // final byte GTE = Constraint.GREATER; // prob.newConstraint("c2").setType(LTE).setRightHandSide( 100.0); // prob.newConstraint("c3").setType(LTE).setRightHandSide(6000.0); Problem prob = new Problem(3, 2); prob.getMetadata().put("lp.isMaximize", "true"); prob.newVariable("x1").setObjectiveCoefficient(13.0); prob.newVariable("x2").setObjectiveCoefficient(23.0); prob.newConstraint("corn").setRightHandSide( 480.0); prob.newConstraint("hops").setRightHandSide( 160.0); prob.newConstraint("malt").setRightHandSide(1190.0); prob.setCoefficientAt("corn", "x1", 5.0); prob.setCoefficientAt("corn", "x2", 15.0); prob.setCoefficientAt("hops", "x1", 4.0); prob.setCoefficientAt("hops", "x2", 4.0); prob.setCoefficientAt("malt", "x1", 35.0); prob.setCoefficientAt("malt", "x2", 20.0); StdOut.println(prob); DenseSimplex lp = new DenseSimplex(prob); StdOut.println(lp.solve()); StdOut.println(lp.getSolution()); } }