var={5795.2 2814.8 1655.8 4967.3 5298.4,
     2814.8 1367.2  804.2 2412.7 2573.5,
     1655.8  804.2  473.1 1419.2 1513.8,
     4967.3 2412.7 1419.2 4257.7 4541.5,
     5298.4 2573.5 1513.8 4541.5 4844.3};
mu={285.00 184.50 105.70 347.50 204.00};

a={  1.00    1.00    1.00     1.00   1.00,
   225.00   50.00   33.33   333.33  25.00,
     1.00    0.10    0.05     0.75   0.10,
    -1.00    0.00    0.00     0.00   0.00,
     0.00   -1.00    0.00     0.00   0.00,
     0.00    0.00   -1.00     0.00   0.00,
     0.00    0.00    0.00    -1.00   0.00,
     0.00    0.00    0.00     0.00  -1.00};

x={0.00, 436.92382, 0.00, 55.38462, 147.69157};
a*x;

g=mu'-0.00002*var*x;
g;
h=-0.00002*var;
h;

astr=a[1 3 4 6,.];
aa=astr'~g;
rref(aa[1:4,.]);
z=null(astr);

print "Projected Gradient" z'g;
print "Projected Hessian" z'h*z;