function y = vcoref(x)
t=x(1);
fr=1;
f=x(2);
y=sin(2*pi*(fr+f)*t);