mosquito(Freq,Dur,Amp,stereo,taper)

function [] = mosquito(Freq,Dur,Amp,stereo,taper)
% Manipulates a sine function and plays it as .wav audio.  
%
% 1. Inspired by the "mosquito ringtone" for younger ears.
% 2. To see how old you(r ears) are, try some freq's from 15k to 18k.
% 3. See commented section at end for music note frequencies to try.
% 4. For really high freq's, turn the volume up on your cpu.
% 5. If values are set for the input arguments, they are played. Otherwise
%    the default values below are used.
%

%% Set default params if not given (play with these values)
s = who;
if ~sum(ismember(s,'Freq'))
    Freq = 14500;       % Frequency of the sine wave (Hertz)
end
Fs = 48000;         % Sample rate of the playback (windows audio standard)
if ~sum(ismember(s,'Dur'))
    Dur = 2;            % Length of the playback (seconds)
end
if ~sum(ismember(s,'Amp'))
    Amp = 0.7;          % Height of the sine wave (scaled 0 to 1)
end
if ~sum(ismember(s,'stereo'))
    stereo = true;      % Logical true for sterio playback, false for mono
end
if ~sum(ismember(s,'taper'))
    taper = true;       % Ramps 0.1 seconds of lead-in / lead-out
end

%% Wave generation (don't mess with these)
step = 2*pi*Freq/Fs;
xmax = 2*pi*Dur*Freq;
x = 0:step:xmax;
y = Amp*sin(x)';

%% Wave modifications (turn on/off with logical parameters above)
if taper
    taperlen = Fs/10;
    t = ones(size(y));
    for i = 0:taperlen
        t(i+1) = i*(1/taperlen);
        t(end-i) = t(i+1);
    end
    y = y.*t;
end 
if stereo % (must be the last mod so both channels get prior mods!)
    y(:,2) = y(:,1);
end

%% Play it
wavplay(y,Fs,'async');


%% Frequencies of Musical Notes
% Taken from http://www.phy.mtu.edu/~suits/notefreqs.html
% ("Middle C" is C4 )
% 
% Note      Frequency(Hz)
%	C0          16.35
%	C#0/Db0		17.32
%	D0          18.35
%	D#0/Eb0		19.45
%	E0          20.6
%	F0          21.83
%	F#0/Gb0		23.12
%	G0          24.5
%	G#0/Ab0		25.96
%	A0          27.5
%	A#0/Bb0		29.14
%	B0          30.87
%	C1          32.7
%	C#1/Db1		34.65
%	D1          36.71
%	D#1/Eb1		38.89
%	E1          41.2
%	F1          43.65
%	F#1/Gb1		46.25
%	G1          49
%	G#1/Ab1		51.91
%	A1          55
%	A#1/Bb1		58.27
%	B1          61.74
%	C2          65.41
%	C#2/Db2		69.3
%	D2          73.42
%	D#2/Eb2		77.78
%	E2          82.41
%	F2          87.31
%	F#2/Gb2		92.5
%	G2          98
%	G#2/Ab2		103.83
%	A2          110
%	A#2/Bb2		116.54
%	B2          123.47
%	C3          130.81
%	C#3/Db3		138.59
%	D3          146.83
%	D#3/Eb3		155.56
%	E3          164.81
%	F3          174.61
%	F#3/Gb3		185
%	G3          196
%	G#3/Ab3		207.65
%	A3          220
%	A#3/Bb3		233.08
%	B3          246.94
%	C4          261.63
%	C#4/Db4		277.18
%	D4          293.66
%	D#4/Eb4		311.13
%	E4          329.63
%	F4          349.23
%	F#4/Gb4		369.99
%	G4          392
%	G#4/Ab4		415.3
%	A4          440
%	A#4/Bb4		466.16
%	B4          493.88
%	C5          523.25
%	C#5/Db5		554.37
%	D5          587.33
%	D#5/Eb5		622.25
%	E5          659.26
%	F5          698.46
%	F#5/Gb5		739.99
%	G5          783.99
%	G#5/Ab5		830.61
%	A5          880
%	A#5/Bb5		932.33
%	B5          987.77
%	C6          1046.5
%	C#6/Db6		1108.73
%	D6          1174.66
%	D#6/Eb6		1244.51
%	E6          1318.51
%	F6          1396.91
%	F#6/Gb6		1479.98
%	G6          1567.98
%	G#6/Ab6		1661.22
%	A6          1760
%	A#6/Bb6		1864.66
%	B6          1975.53
%	C7          2093
%	C#7/Db7		2217.46
%	D7          2349.32
%	D#7/Eb7		2489.02
%	E7          2637.02
%	F7          2793.83
%	F#7/Gb7		2959.96
%	G7          3135.96
%	G#7/Ab7		3322.44
%	A7          3520
%	A#7/Bb7		3729.31
%	B7          3951.07
%	C8          4186.01
%	C#8/Db8		4434.92
%	D8          4698.64
%	D#8/Eb8		4978.03