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

u=cumsum([phi,f],2);
y=sin(2*pi*u(1:end-1));
phi=mod(u(end),1);
	
