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

t=ceil(1/f);
y=mod(phi+(0:n-1)/t,1)<duty;
phi=mod(phi+n/t,1);
	
