Feasible Direction Method
function fdm input Real u0[size(u0, 1), 1]; input Real H[size(u0, 1), size(u0, 1)]; input Real b[size(u0, 1), 1]; input Real A[size(A, 1), size(u0, 1)]; input Real a[size(A, 1), 1]; input Real Aeq[size(Aeq, 1), size(u0, 1)]; input Real aeq[size(Aeq, 1), 1]; constant Integer nu = size(u0, 1); constant Integer na = size(A, 1); constant Integer naeq = size(Aeq, 1); output Real u[size(u0, 1), 1]; end fdm;