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