Example 1

% Circuit parameters
R = 100;
C = 2e-6;
L = 0.3e-3;

% Input signal
Vin = 1;
f = 1000;
w = 2*pi*f;

% results
alpha = R/(2*L);
w0 = 1/sqrt(L*C);
f0  = w0/(2*pi);
fprintf('LC frequency %g Hz\n',f0);
Vout = w0^2*Vin/(w0^2-w^2+2j*w*alpha);
[rho phi] = complex2phasor(Vout);
fprintf('output signal: amplitude %g V, phase %g degress\n',rho,phi);
LC frequency 6497.47 Hz
output signal: amplitude 0.628406 V, phase -52.1555 degress