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

w=2*pi*f;
y=sin(2*pi*phi + w*(0:n-1));
phi=mod(phi+f*n,1);
	
