% Problema de ejemplo de la clase de INEL 4102 Feb 26, 2001 (p1.m)
% estas lineas deben aparecer si ejecuto la instruccion
% >> help p1

% de la linea en blanco en adelante no debe aparecer nada al ejecutar help

Va  = 15;               % voltaje de bateria
R1  = 1e3;              % primera resistencia
R2  = 4e3;              % segunda resistencia
C   = 6.25e-9;          % capacitancia
L   = 0.625e-3;         % inductancia
iL0 = 3e-3;             % corriente incial por el inductor

% constantes del numerador
n2 = iL0;
n1 = iL0 / (R2 * C);
n0 = Va / (R2 * L * C);

% constantes del denominador
d1 = 1 / (R2 * C);
d0 = 1 / (L * C);

% raiz con parte imaginaria positiva
rtpos = -d1 / 2 + 1 / 2 * sqrt(d1^2 - 4 * d0);

% como se que es complejo, entonces (rtpos = ro + j * wo)
ro = real(rtpos);
wo = imag(rtpos);

% para determinar intervalo de tiempo para la figura
tau = abs(1 / ro);
t = 4 * tau * (0:.0025:1);    % vector de tiempo

% los coeficientes de una expansion en fracciones parciales
k1 = - n0 / d0;
k2 = (rtpos^2 * n2 + rtpos * n1 - n0) / (rtpos * 2 * j * wo);  % j = sqrt(-1)

% la solucion en tiempo
f = k1 + 2 * abs(k2) * exp(ro * t) .* cos(wo * t + angle(k2));

% dibujar la figura
figure(1)           % abre una pantalla de dibujo nueva o usa la primera
clf
plot(t,f)           % hace una grafica de f vs t
xlabel('t (s)')     % texto para eje horizontal
ylabel('iL(t) (A)') % texto para eje vertical
title('Example exercise for INEL 4102 - Feb 23, 2001')
grid('on')          % dibuja lineas entrecortadas dentro de la grafica