r1= 0.2;
r2= linspace(0,100, 5000);
c= 53*1e-6;
l= 83*1e-6;
Uin = 1;
%w = 1 / sqrt(c*l)  % идеален случай, не наистина
w = sqrt((r2.^2) ./ (r2.^2.*l*c+l^2))
f = w ./ (2*pi);
Zreal = r1 -((r2.*w.^2*l^2)./(r2.^2-w.^2*l^2));
Zimag = -1./(w.*c) + (r2.^2.*w.*l)./(r2.^2-w.^2*l^2);

Irezonans = Uin ./ Zreal;
UL = w*l.*Irezonans;
UC = (1./(w.*c)).*Irezonans;
Q = UL ./ Uin;
Q2 = UC ./ Uin;

%plot(r2, Q), hold on
%plot(r2, Q2)
plot(r2, UC)
% plot(f, Zimag), hold on
% plot(w([1,end])/2/pi,[0,0],'k'), hold on
% plot(f, Q)
% grid on, ylim([-200 200])
% xlabel('Freq,Hz'),ylabel('Im(Z)')