niu=0.1:0.
01:2;
k=0.18;
Q1=8;
Q2=12;
Q3=16;
A=1-(1./(niu.^2));
B=(1-k)/k;
C1=Q1.*(niu-(1./niu)); % rezolutia aleasa da vector de 20 de elemente
C2=Q2.*(niu-(1./niu));
C3=Q3.*(niu-(1./niu));
D=(1-k)/(2*k);
E1=1./(sqrt((1+B*A).^2 + (C1.*(1+D*A)).^2));
E2=1./(sqrt((1+B*A).^2 + (C2.*(1+D*A)).^2));
E3=1./(sqrt((1+B*A).^2 + (C3.*(1+D*A)).^2));
plot(niu,E1,'b')
hold on
plot(niu,E2,'r')
hold on
plot(niu,E3,'g')
legend('Q=8','Q=12','Q=16')