function [y,phi]=block(o,phi,n,f)
% block - generate sine wave

t=2*pi*f*(0:n-1);
y=sin(repmat(2*pi*phi,[1,n]) + t);
phi=mod(phi+f*n,1);
	
