samer@34: function [y,phi]=block(o,phi,n,duty,f)
samer@34: % block - generate square wave
samer@34: 
samer@34: t=ceil(1/f);
samer@34: y=mod(phi+(0:n-1)/t,1)<duty;
samer@34: phi=mod(phi+n/t,1);
samer@34: