R1=300; R2=100; C1=5e-6; C2=10e-6; C3=2e-6; L=1;
f=logspace(0,3,300); omega=2*pi*f;
ZL=j*omega*L; ZC1=1./(j*omega*C1); ZC2=1./(j*omega*C2);
ZC3=1./(j*omega*C3); Z1=R1+ZL+ZC1;
Z=R2+ZC2+Z1.*ZC3./(Z1+ZC3);
loglog(f,abs(Z),'LineWidth',2);
set(gca, 'FontSize', 14);
xlabel('f (Hz)', 'FontSize',14); ylabel('|Z| (\Omega)','FontSize', 14);
grid on;
figure;
semilogx(f,angle(Z)*180/pi,'LineWidth',2);
set(gca, 'FontSize', 14);
xlabel('f (Hz)', 'FontSize',14); ylabel('\phi (deg)', 'FontSize', 14);
grid on;