Je n'ai trouvé que ceci qui fonctionnât à peu près, avec un ajustement manuel:
manuel :
\documentclass{article}
\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{lmodern}
\usepackage{amsmath, amsfonts}
\usepackage{systeme}
\begin{document}
\sysdelim|.
\sysaddeqsign{~}
\sysequivsign{~}{\in}
\begin{align*}
\max z& = x_1-x_2\\
\mathrm{s. c.} & \phantom{ =}
\raisebox{-0.77\height}{$ \,\systeme{\frac{1}{3}x_1 + x_2 = 2,
-2x_1 + 5x_2 <= 7,
x_1 + x_2 <= 4,
x_1 >= 0,
x_2 ~ \mathbb{R}}$}
\end{align*}
\end{document}
![alt text][1]
[1]: http://texnique.fr:80/osqa/upfiles/optimisationlineaire.png<img src="/upfiles/optimisationlineaire.png" style="width: 100%;"/>