a0=1;
a2=(1-(2*a0))/4
a1=(a0+(2*a2))/2
T=0.001;
f=1:1:1/T;

y= a0+2*a1*cos(2*pi*f*T)+2*a2*cos(2*pi*f*2*T);
z= 2*sin(f);

#plot(f,z);

k=z+y;

#plot(f,k);
plot(f,z, f,k, f,y);

