V0=325;
R=10;
C=0.47e-6;
L=6;
w=2.*pi*[1:1:200];
w0=1./(2.*pi.*sqrt(L.*C))

t=[0:0.0001:1]-1.01;

for i = 1:1:length(w)
Imax(i)=V0./(sqrt(R.^2+(w(i).*L-1./(w(i).*C)).^2));
fi(i) = atan((w(i).*L-1./(w(i).*C))./R);
t=t+1.01; 
I=Imax(i).*cos(w(i).*t-fi(i));
figure(1), plot(t, I), hold on, 
set(gca,'FontSize',16,'FontWeight','bold')
xlabel('Czas (s)')
ylabel('I (A)')
end
