function [t]=CapacitorSec31(vs, R, C, vd)
 v=0;
 t=0;
 dt=0.01;
    while v<=vd
          v=v+(vs-v)/(R*C)*dt;
          t=t+dt;
    end
end