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

N = 30;   % number of nodes

V = zeros(N,1);
dV = zeros(N,1);

t = 0;
dt = 0.02;

while(t < 100)

   if (t < 2) V0 = 100 * ( ( sin(0.5*pi*t) )^2 ); 
      else V0 = 0;
      end

   ddV(1) = 40*V0 - 80*V(1) + 40*V(2) - 0.01*dV(1);

   for i=2:N-1
      ddV(i) = 40*V(i-1) - 80*V(i) + 40*V(i+1) - 0.01*dV(i);
      end

   ddV(N) = 40*V(N-1) - 40*V(N) - 0.01*dV(N);
%                                 ^^^ 
%                              change this term


   for i=1:N
      dV(i) = dV(i) + ddV(i)*dt;
      V(i) =  V(i) +  dV(i)*dt;
      end
   t = t + dt;
 
   plot([0:N],[V0;V],'.-');
   ylim([-150,250]);
   clc
   disp(t)
   pause(0.01);

   end