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