;; Carregamos o pacote "simplex":

(%i1) load("simplex")$

;; Definimos a função objetivo e a lista de restrições:

(%i2) obj: x1+4*x2$
(%i3) res: [x1 +  x3 = 9, x2 + 2*x4 = 3]$

;; Usamos "nonegative_lp=true", de outra forma o Maxima 
;; permitirá que as variáveis sejam negativas:

(%i4) maximize_lp(obj,res),nonegative_lp=true;
(%o4)               [21, [x4 = 0, x2 = 3, x3 = 0, x1 = 9]]

Ali estão o valor objetivo (21) e os valores das variáveis
para a solução ótima.