function y = vco1(x)
t = x(1); % seconds
fi = (x(2)/180)*pi; % radian
fr = x(3); % Hertz
df = x(4);
f = fr+df;
w = 2*pi*f; % rad/sec
y = sin(w*t + fi);