samer@34: function [y,phi]=osc(f,phi0)
samer@34: % osc - sequence of oscillator output
samer@34: %
samer@34: % osc :: seq [[N]] ~'persample normalised frequencies' -> seq [[N]] ~'sinusoidal data'.
samer@34: %
samer@34: % osc :: 
samer@34: %    seq([[N]]) ~'persample normalised frequencies', 
samer@34: %    real       ~'initial phase' 
samer@34: % -> seq([[N]]) ~'sinusoidal data',
samer@34: %    seq([[N]]) ~'phase data'.
samer@34: 
samer@34: if nargin<2, phi0=zeros(size(f,1),1); end
samer@34: y=srdata(sine,phi0,f);