%% Chapter 15 Problem 4
% Thrasyvoulos Karydis

Cc = 5e-12;
Cm = 20e-15;
Lm = 3e-3;
Rm = 6;

G = tf([Cm*Lm Cm*Rm 1],[Cc*Cm*Lm Cc*Cm*Rm Cm+Cc 0]);

[mag, phase, wout] = bode(G);

mag = squeeze(mag);
phase = squeeze(phase);

Z = mag.*exp(1i*phase);

figure;
subplot(3,1,1);
semilogx(wout, real(Z))
subplot(3,1,2);
semilogx(wout, imag(Z))
subplot(3,1,3);
semilogx(wout, phase)