% 10-stage RC Filter
% Lecture #4

V = rand(10,1);
dV = 0*V;

dt = 0.01;
t = 0;
V0 = 0;

REF = 1;
Z = 0;

DATA = [];

while(t < 2)
   
   dV(1) = 10*V0   - 20.2*V(1) + 10*V(2);
   dV(2) = 10*V(1) - 20.2*V(2) + 10*V(3);
   dV(3) = 10*V(2) - 20.2*V(3) + 10*V(4);
   dV(4) = 10*V(3) - 20.2*V(4) + 10*V(5);
   dV(5) = 10*V(4) - 20.2*V(5) + 10*V(6);
   dV(6) = 10*V(5) - 20.2*V(6) + 10*V(7);
   dV(7) = 10*V(6) - 20.2*V(7) + 10*V(8);
   dV(8) = 10*V(7) - 20.2*V(8) + 10*V(9);
   dV(9) = 10*V(8) - 20.2*V(9) + 10*V(10);
   dV(10) = 10*V(9) - 10.2*V(10);

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

   plot([0:10], [V0;V], 'b.-');
   ylim([0,1.2]);
   xlim([0,10]);
   pause(0.01);
   DATA = [DATA ; V'];
end
   
%pause(2);
%plot(DATA)