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