w = 2 * pi * 50
R = 5
C1 = 0.0005
C2 = 0.0009
C3 = 0.0011
m = 50
n = 100
U = 1
T = R*C1
t = 0:T/n:m*2*C3
A = asin (w / (sqrt ((1/T)ˆ2 + wˆ2))) ;
it = U/R * (1/(sqrt(1+(w*T)ˆ2))) * (-exp(-t/T) * sin(A) + cos(w*t-A) * ...
w * T);
plot(t,i t)
hold on
T = R*C2
t = 0:T/n:m*2*C3
A = asin (w/(sqrt((1/T)ˆ2 + wˆ2)));
it2 = U/R * (1/(sqrt(1+(w*T)ˆ2))) * (-exp(-t/T) * sin(A) + cos (w * t - A) ...
* w * T);
plot (t,it2, 'r - - ')
 hold on
 T = R*C3
 t = 0:T/n:m*2*C3
 A = asin (w/(sqrt((1/T)ˆ2 + wˆ2)));
 it3 = U/R * (1/( sqrt(1+(w*T)ˆ2)))*(-exp(-t/T) * sin(A) + cos (w * t - A) ...
* w * T);
 plot ( t, it3, 'g : ')
 legend('C1 : 5 nmu F' ,'C2 : 9 nmu F' ,'C3 : 11 nmu F' ,1);
 xlabel('t');
ylabel('i(t)');
hold off