function y=mymod(x,Fc,Fs)
if length(Fs)<2;ip=0;else;ip=Fs(2);Fs=Fs(1);end;
[r,c]=size(x);
if r==1;x=x(:);end;
t=(0:1/Fs:((size(x,1)-1)/Fs))';
y=x.*cos(2*pi*Fc*t+ip);
if r==1;y=y.';end;