% 10-stage RC Filter
% Spring 2018 Model
% Rs = 1,  Rp = 80, C = 0.25uF

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

dt = 0.01;
t = 0;

while(t < 100)

   Ref = 1 * (sin(0.5*t) > 0);
    
   V0 = Ref;
   
   dV(1)  = 4*V0   - 8.05*V(1) + 4*V(2);
   dV(2)  = 4*V(1) - 8.05*V(2) + 4*V(3);
   dV(3)  = 4*V(2) - 8.05*V(3) + 4*V(4);
   dV(4)  = 4*V(3) - 8.05*V(4) + 4*V(5);
   dV(5)  = 4*V(4) - 8.05*V(5) + 4*V(6);
   dV(6)  = 4*V(5) - 8.05*V(6) + 4*V(7);
   dV(7)  = 4*V(6) - 8.05*V(7) + 4*V(8);
   dV(8)  = 4*V(7) - 8.05*V(8) + 4*V(9);
   dV(9)  = 4*V(8) - 8.05*V(9) + 4*V(10);
   dV(10) = 4*V(9) - 4.5*V(10);

   V = V +  dV*dt;
   t = t + dt;

   hold off
   plot([0,10],[0,1.2],'w.');
   hold on
   plot([0:10], [V0;V], '.-');
   pause(0.01);
   end