function [y,phi]=osc(f,phi0)
% osc - sequence of oscillator output
%
% osc :: seq [[N]] ~'persample normalised frequencies' -> seq [[N]] ~'sinusoidal data'.
%
% osc :: 
%    seq([[N]]) ~'persample normalised frequencies', 
%    real       ~'initial phase' 
% -> seq([[N]]) ~'sinusoidal data',
%    seq([[N]]) ~'phase data'.

if nargin<2, phi0=zeros(size(f,1),1); end
y=srdata(sine,phi0,f);
