% ECE 111 Lecture #8
% 3-Stage RLC Filter

V0 = 100;
V1 = 0;
V2 = 0;
V3 = 0;

dV1 = 0;
dV2 = 0;
dV3 = 0;

V = [];

t = 0;
dt = 0.02;

while(t < 10)

   ddV1 = 50*V0 - 100*V1 + 50*V2 - 0.01*dV1;
   ddV2 = 50*V1 - 100*V2 + 50*V3 - 0.01*dV2;
   ddV3 = 50*V2 -  50*V3         - 0.01*dV3;

   dV1 = dV1 + ddV1*dt;
   dV2 = dV2 + ddV2*dt;
   dV3 = dV3 + ddV3*dt;

   V1 = V1 + dV1*dt;
   V2 = V2 + dV2*dt;
   V3 = V3 + dV3*dt;
   t = t + dt;
   
   plot([0,1,2,3],[V0,V1,V2,V3],'.-');
   ylim([-300,300]);
   clc
   disp(t)
   pause(0.01);

   V = [V; V1, V2, V3];

   end

t = [1:length(V)]' * dt;
plot(t,V);