function [x,fval] = abslinprog(f,A,b,Aeq,beq,LB,UB,X0,OPTIONS)

[xp,fvp] = linprog( f,A,b,Aeq,beq,LB,UB,X0,OPTIONS);
[xn,fvn] = linprog(-f,A,b,Aeq,beq,LB,UB,X0,OPTIONS);

fvp = abs(fvp);
fvn = abs(fvn);
if fvp > fvn
  fval = fvp;
  x = xp;
else
  fval = fvn;
  x = xn;
end

return