sin i(t)
w = 2 * pi * 50;
R = 5;
C1 = 0.0005;
C2 = 0.0008;
C3 = 0.0010;
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,it,'b');
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 : 8 nmu F' ,'C3 : 10 nmu F');
xlabel('t');
ylabel('i(t)');
hold off





sin Ur(t)
w = 2 * pi * 50;
R = 5;
C1 = 0.0005;
C2 = 0.0008;
C3 = 0.0010;
m = 50;
n = 100;
U = 1;

T = R*C1;
t = 0:T/n:m*2*C1;
A = asin(w/(sqrt((1/T).^2 + w.^2)));
ur = U*(1/(sqrt(1+(w*T).^2))) * (-exp(-t/T) * sin(A) + cos(w*t - A)*w*T);
plot(t,ur,'b');
hold on

T = R*C2;
t = 0:T/n:m*2*C2;
A = asin(w/(sqrt((1/T).^2 + w.^2)));
ur2 = U * (1/(sqrt(1+(w*T).^2))) * (-exp(-t/T) * sin(A) + cos(w*t - A) *w*T);
plot (t,ur2,'r--')
hold on

T = R*C3;
t = 0:T/n:m*2*C3;
A = asin(w/(sqrt((1/T).^2 + w.^2)));
ur3 = U * (1/(sqrt(1+(w*T).^2))) * (-exp(-t/T) * sin(A) + cos(w*t - A)*w*T);
plot(t, ur3, 'g: ')
legend('C1 : 5 nmu F' ,'C2 : 8 nmu F' ,'C3 : 10 nmu F');
xlabel('t');
ylabel('U_{R}(t)');
hold off




sin Uc(t)
w = 2 * pi * 50;
R = 5;
C1 = 0.0005;
C2 = 0.0008;
C3 = 0.0010;
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))) ;
uc = U *(1/(sqrt(1+(w*T).^2))) * (-exp(-t/T) * sin(A) + sin(w*t-A)*w*T);
plot(t,uc,'b');
hold on

T = R*C2;
t = 0:T/n:m*2*C3;
A = asin(w/(sqrt((1/T).^2 + w.^2)));
uc2 = U * (1/(sqrt(1+(w*T).^2))) * (-exp(-t/T) * sin(A) + sin(w*t - A) *w*T);
plot (t,uc2,'r--')
hold on
T = R*C3;

t = 0:T/n:m*2*C3;
A = asin(w/(sqrt((1/T).^2 + w.^2)));
uc3 = U * (1/(sqrt(1+(w*T).^2)))*(-exp(-t/T) * sin(A) + sin(w*t - A)*w*T);
plot(t, uc3, 'g: ')
legend('C1 : 5 nmu F' ,'C2 : 8 nmu F' ,'C3 : 10 nmu F');
xlabel('t');
ylabel('U_{C}(t)');
hold off