Simplex Algorithm

[Albert Chan]

(11-12-2023 05:40 PM)ftneek Wrote:  Reminder: to avoid incorrect answers for the Dual solution, the 15/7 should really be a 0. This likely just requires either a small modification or new sol() command.

Code:

sol(a,v):=
BEGIN
LOCAL x,b,c, j,k;
x:=makelist(v); /* all zeroes */
b:=col(a,0);    /* last column */
FOR j FROM 1 TO v DO
 IF abs(trn(c:=col(a,j)))!=1.0 THEN continue END
 IF (k:=contains(c,1.0))==0 THEN continue END
 IF b[k] THEN x[j]:=b[k]; b[k]:=0 END
END;
RETURN x;
END;