Add fourth example.

This example is based on issue sloisel/numeric#43.

Author: albertosantini

test/fourth-test.R

@@ -0,0 +1,8 @@
+require(quadprog)
+
+Dmat <- matrix(c(13,18,-6, 18,27,-9, -6,-9,4),3,3)
+dvec <- c(-4,0,-100)
+Amat <- matrix(c(0,0,-1),3,1)
+bvec <- c(-25)
+
+solve.QP(Dmat,dvec,Amat,bvec=bvec,meq=1)

test/fourth-test.js

@@ -0,0 +1,27 @@
+"use strict";
+
+var vows = require("vows"),
+    assert = require("assert"),
+    t = require("./fourth");
+
+vows.describe("Test 4").addBatch({
+    "with the fourth test": {
+        topic: t.fourth,
+
+        "we get the result": function (topic) {
+            assert.equal("", topic.message);
+            assert.equal(-4.000000000000016, topic.solution[1]);
+            assert.equal(11.000000000000007, topic.solution[2]);
+            assert.equal(25, topic.solution[3]);
+            assert.equal(2804.500000000001, topic.value[1]);
+            assert.equal(-3.999999999999967, topic.unconstrained_solution[1]);
+            assert.equal(-30.666666666666707, topic.unconstrained_solution[2]);
+            assert.equal(-100.00000000000006, topic.unconstrained_solution[3]);
+            assert.equal(2, topic.iterations[1]);
+            assert.equal(0, topic.iterations[2]);
+            assert.equal(1, topic.iact[1]);
+            assert.equal(125, topic.Lagrangian[1]);
+        }
+    }
+
+}).export(module);

test/fourth.js

@@ -0,0 +1,47 @@
+"use strict";
+
+var qp = require("../lib/quadprog");
+
+function fourth() {
+    var dmat = [], dvec = [], amat = [], bvec = [], meq;
+
+    dmat[1] = [];
+    dmat[2] = [];
+    dmat[3] = [];
+    dmat[1][1] = 13;
+    dmat[2][1] = 18;
+    dmat[3][1] = -6;
+    dmat[1][2] = 18;
+    dmat[2][2] = 27;
+    dmat[3][2] = -9;
+    dmat[1][3] = -6;
+    dmat[2][3] = -9;
+    dmat[3][3] = 4;
+
+    dvec[1] = -4;
+    dvec[2] = 0;
+    dvec[3] = -100;
+
+    amat[1] = [];
+    amat[2] = [];
+    amat[3] = [];
+    amat[1][1] = 0;
+    amat[2][1] = 0;
+    amat[3][1] = -1;
+
+    bvec[1] = -25;
+
+    meq = 1;
+
+    return qp.solveQP(dmat, dvec, amat, bvec, meq);
+}
+
+module.exports = {
+    name: "fourth test",
+    fourth: fourth,
+    tests: {
+        "fourth": function () {
+            fourth();
+        }
+    }
+};