lucro

mi =  1.2:0.2:6
L =   [5.4545 2.9167 2.0513 1.6071 1.3333 1.1458 1.0084 0.9028 0.8187 0.7500 0.6926 0.6439 0.6020 0.5655 0.5333 0.5048 0.4793 0.4563 0.4356 0.4167 0.3994 0.3835 0.3689 0.3554 0.3429];
Lq =  [3.7879 1.4881 0.8013 0.4960 0.3333 0.2367 0.1751 0.1335 0.1044 0.0833 0.0676 0.0557 0.0465 0.0392 0.0333 0.0286 0.0248 0.0216 0.0189 0.0167 0.0148 0.0132 0.0118 0.0106 0.0095];
W =   [2.7273 1.4583 1.0256 0.8036 0.6667 0.5729 0.5042 0.4514 0.4094 0.3750 0.3463 0.3220 0.3010 0.2827 0.2667 0.2524 0.2397 0.2282 0.2178 0.2083 0.1997 0.1918 0.1845 0.1777 0.1714];
Wq =  [1.8939 0.7440 0.4006 0.2480 0.1667 0.1184 0.0875 0.0668 0.0522 0.0417 0.0338 0.0279 0.0232 0.0196 0.0167 0.0143 0.0124 0.0108 0.0095 0.0083 0.0074 0.0066 0.0059 0.0053 0.0048];
rho = [0.5556 0.7143 0.6250 0.5556 0.5000 0.4545 0.4167 0.3846 0.3571 0.3333 0.3125 0.2941 0.2778 0.2632 0.2500 0.2381 0.2273 0.2174 0.2083 0.2000 0.1923 0.1852 0.1786 0.1724 0.1667];


//caso 1 - a=0.5 e b =50
a =   0.5
b =   50
R = b/W;
D = a*mi^2;
L = R-D';
figure(0)
plot2d(mi,[R D' L],[1 2 5], leg="Receita@Despesa@Lucro");

//caso 2 - a=1 e b =50
a =   1
b =   50
R = b/W;
D = a*mi^2;
L = R-D';
figure(1)
plot2d(mi,[R D' L],[1 2 5], leg="Receita@Despesa@Lucro");

//caso 3 - a=1 e b =25
a =   1
b =   25
R = b/W;
D = a*mi^2;
L = R-D';
figure(2)
plot2d(mi,[R D' L],[1 2 5], leg="Receita@Despesa@Lucro");

//Gráfico mi x L
figure(3)
plot2d(mi,L, leg="Nº médio de pessoas no sistema em função de mi");
//Gráfico mi x Lq
figure(4)
plot2d(mi,Lq, leg="Nº médio de pessoas na fila em função de mi");
//Gráfico mi x W
figure(5)
plot2d(mi,W, leg="Tempo médio de espera no sistema em função de mi");
//Gráfico mi x Wq
figure(6)
plot2d(mi,Wq,leg="Tempo médio de espera na fila em função de mi");
//Gráfico mi x rho
figure(7)
plot2d(mi,rho, leg="Taxa de utilização em função de mi");