%otimizaçõ de lucro
%criar malha de -5 a 25 com um incremento de 0.5 para as variaveis x1 e x2
[x1,x2]=meshgrid(-5:0.5:25.0,-5:0.5:25.0);
%entra funcoes para o problema de maximização de lucro
f=400*x1+600*x2;
g1=x1+x2-16;
g2=x1/28+x2/14-1;
g3=x1/14+x2/24-1;
g4=-x1;
g5=-x2;
%inicializacao, sem ponto e virgula no final
    cla reset
    axis auto %minimo e maximo para os eixos determ. automaticamente
              %limites para eixos x e y podem tambem ser dados por
              %axis([xmin xmax ymin y max])
    xlabel('x1'),ylabel('x2') %da nome aos eixos
    title('Maximizacao de lucro') %titulo do problema
    hold on %mantem o grafico atual e propriedades dos eixos daqui p/frente
%uso do comando "contour" para plotar as restriçoes e funcao custo
    pause
    cv1=[0 0]; %especifica  valores de contorno, 0
    const1=contour(x1,x2,g1,cv1,'k'); %plota contorno p/ g1; k=preto
    clabel(const1) %coloca o valor do contorno no grafico
    text(1,16,'g1') %escreve g1 no ponto (1,16)
    pause
    cv2=[0 0]; 
    const2=contour(x1,x2,g2,cv2,'k');
    clabel(const2) 
    text(23,3,'g2')
    pause
    const3=contour(x1,x2,g3,cv2,'k'); 
    clabel(const3) 
    text(1,23,'g3')
    pause
    cv3=[0 0];
    const4=contour(x1,x2,g4,cv3,'k');
    clabel(const4) 
    text(.25,20,'g4')
    pause
    const5=contour(x1,x2,g5,cv3,'k');
    clabel(const5)
    text(19,.5,'g5')
    text(1.5,7,'regiao viável')
    pause
    fv=[2400, 4800, 7200, 8800]; %define 4 contornos para a funcao lucro
    fs=contour(x1,x2,f,fv, '.'); %linhas coloridas
    clabel(fs)
    hold off %fim da sequencia de plotagens
             %plotagens seguintes aperecerao em outras janelas