function [ sinusoid ] = makesinusoid(A,freq,phi,fs,dur)
t=[1:1/fs:dur];
sinusoid=A*cos(2*pi*freq*t+phi);