clc
clear

r=1000
c=10*(10^(-6))
tau=r*c
t=0:0.0001:50*10^(-3)
a=1/(1+4*pi*pi)
b=(2*pi)/(1+4*pi*pi)
d=-1/(1+4*pi*pi)

new=a*cos(200*pi*t)+b*sin(200*pi*t)+d*exp(-100*t)
figure(1)
plot(t,new,'k')
xlabel('Time')
ylabel('Volt of source (Vs) ')