% Sinus sweep

clear all;

tstart = -1;
tstop  = 1;
fmax   = 10;

t = tstart:.0001:tstop;

y = atan(7*t);

y = (y - y(1));
y = y / y(length(y));

subplot(3,1,1);
plot(t, y);

dy = diff(y)./diff(t);

subplot(3,1,2);
plot(t(1:(length(t)-1)), dy);

s = sin(100*y);

subplot(3,1,3);
plot(t, s);
