1. definicja ograniczeń

OG:[13/6*x1+x2<=13,x1+4*x2<=20,x1>=0,x2>=0];

2. wyrysowanie ograniczeń

wxdraw2d(
    implicit(OG[1],x1,-1,7,x2,-1,7),
    implicit(OG[2],x1,-1,7,x2,-1,7),
    implicit(OG[3],x1,-1,7,x2,-1,7),
    implicit(OG[4],x1,-1,7,x2,-1,7)        
)$

3. wyrysowanie funkcji celu

wxplot3d(x1+x2,[x1,-1,7],[x2,-1,7]);

4. albo alternatywnie połączenie w jednej scenie  ograniczeń i konturowego wykresu funkcji celu 

z0:-5;

z0:-10;
wxdraw3d(    
    contour_levels=10,
    contour=both,
    surface_hide=false,wired_surface=false,
              user_preamble="set for [i=1:8] linetype i dashtype i linecolor 0",
    explicit(x1+x2,x1,-1,7,x2,-1,7),
    parametric(1*t+0,-13/6*t+13,z0,t,-1,7),
    parametric(1*t+0,-1/4*t+5,z0,t,-1,7),
    parametric(1*t+0,0*t+0,z0,t,-1,7),
    parametric(0*t+0,1*t,z0,t,-1,7),
    xlabel="x1",
    ylabel="x2",
    xrange=[-1,7],
    yrange=[-1,7]
);

solve([13/6*x1+x2=13,x1+4*x2=20], [x1,x2]);