close all
clear all

VT = 26e-3;
REIDCQvec=[0.5 1 3 5 10];
colores = {'red', 'blue', 'green', 'black', 'magenta'};

xh = 20*(((1/20):0.005:1));
xl = 0.01:0.05:0.99;
x = [xl xh];

I0 = besseli(0,x);
I1 = besseli(1,x);

GmSobregm = (2*I1)./(x.*I0);

GmSobregmqvec = REIDCQvec;

figure
hold on
for i=1:length(REIDCQvec)
   GmSobregmqvec = ((2*I1)./(x.*I0)).*(1+((VT*log(I0))./(REIDCQvec(i))));
	plot(x,GmSobregmqvec, colores{i})
end
    legend(num2str(REIDCQvec(:)))
hold off
title('Gm sobre gmQ para RE*IDCQ variando entre 0.5V y 10V')
grid
xlabel('x')
ylabel('Gm/gmQ')