lambda=2:2:22;
lambda=lambda';
ts=l/c;
rho=lambda.*ts;
traficofrecido=lambda.*ts./nserv;
pperdida=((1-rho)./(1-rho.^(K+1))).*rho.^(K);
lambdacursada=lambda.*(1-pperdida);
traficursado=lambdacursada.*ts./nserv;

figure(1)
plot([traficofrecido, traficursado, pperdida]);
title("Teoricos");
%Simulacion
pp=load('report_MM1K3.txt');
pp=pp(:,10);
tofre=lambda.*ts;
tcurs=tofre.*(1-pp);
figure(2)
plot([tofre, tcurs, pp]);
title("Simulado");