function ret = sinusoid(fs, f, phase, n)
    ret=zeros(1,n);
    for i=1 : n
        ret(1,i) =  sin(phase + ((2*pi*f)/fs)*i);
    end
end


sinus = sinusoid(4000, 10, 0.5*pi, 800);
plot(sinus);
figure;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

n=800;
sinus1 = sinusoid(16000, 144, 0, n);
sinus2 = sinusoid(16000, 56, 0, n);

sinus3=sinus1+sinus2;

csinus3=zeros(1, 10 );
for i=1 : size(csinus3,2)
    csinus3(1,i)=sinus3(i*80);
end

hold all;
plot(sinus1);
plot(sinus2);
plot(sinus3);
plot([0:80:800](2:end),csinus3);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


cv4