function [y,phi]=block(o,phi,n,f)
% block - generate block of sine wave samples
%
% sine/block :: 
%    sine      ~'sine generator object',
%    real      ~'phase at start of block',
%    N:natural ~'number of samples to generate',
%    real      ~'frequence in cycles per sample'
% -> [[1,N]]   ~'block of output samples',
%    real      ~'phase at end of block'.

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