function y = vcomat(x)
t=x(1);
phi=0.1;
fr = 100;
f=x(2);
y=sin((2*pi*(fr+f)*t) + phi);