Matlab - Potencia Optica con 2 plots

%Potencia interna generada en el dispositivo LED

clear all;

close all;

tuo_r = input('Enter the radiative recombination time in nano seconds: ');

tuo_nr = input('Enter the non-radiative recombination time in nano seconds: ');

C = input('Enter the velocity of light in free space: ');

%ppp

h = input('Enter the value of Planks constant: ');

q = input('Enter the value of charge of one electron: ');

I = input('Enter the value of Drive current in Amperes: ');

Etta_internal = 1/(1+(tuo_r/tuo_nr));

%Pot. Opt. vs. Long. Onda

n=1;

for Lambda=800:50:1600

	Pinternal= Etta_internal*(h*C/q)*(I/(Lambda*(10^(-9))));

	PO= Pinternal*1000;

	Pot(n)=PO;

	n=n+1;

	fprintf('Pot. Ópt. entreg. LED en mW. para %d nm. = %f\n ',Lambda, PO);

end

Lambda=800:50:1600;

figure(1)

plot(Lambda,Pot,'b--','linewidth',3);

title('POTENCIA INTERNA GENERADA EN EL DISPOSITIVO LED');

xlabel('Longitud de Onda (nm.)');

ylabel('Potencia Óptica (mW.)');

grid;

legend('CURVA Long. Onda vs. Pot. Ópt.');

%Pot. Opt. vs. Corr. Eléctr.

Long = input('Ingrese la Long. Onda en nm.: ');

m=1;

for Int=20:5:50

	Pinternal= Etta_internal*(h*C/q)*(Int*(10^(-3))/(Long*(10^(-9))));

	PO= Pinternal*1000;

	Pot(m)=PO;

	m=m+1;

	fprintf('Pot. Ópt. entreg. LED en mW. para %d Amp. = %f\n ',Int, PO);

end

Int=20:5:50;

figure(2)

plot(Int,Pot,'r--','linewidth',3);

title('POTENCIA INTERNA GENERADA EN EL DISPOSITIVO LED');

xlabel('Corriente Eléctrica (Amp.)');

ylabel('Potencia Óptica (mW.)');

grid;

legend('CURVA Corr. Eléctr. vs. Pot. Ópt.');