Cómo resolver el circuito de segundo orden mediante el lenguaje de programación matlab

Como no diste un valor determinado para w, supuse que era la frecuencia eléctrica. El código es el siguiente

R1=22;L=50e-3;R2=. 33;C=100e- 6;w=2*pi*50;Uc=200*exp(1i*90*pi/180)/(R2*1/(1i*w*C)/(R2+1/( 1i*w*C) )+R1+1i*w*L)*R2*1/(1i*w*C)/(R2+1/(1i*w*C));

fplot(@(t) (200*cos(w*t+90*pi/180)),[0,0.08],'r')

espera

fplot (@(t)( abs(Uc)*cos(w*t+angle(Uc))),[0,0.08],'b')

legend('Us(t)', 'Uc(t)' )

cuadrícula activada