function [phase,magnitude] = phamag(wave1,wave2,frequency,T,length)
%PHAMAG Summary of this function goes here
%   Detailed explanation goes here
wave = [wave1;wave2];
fft1 = fft(wave);


    function g = phasefft(wave,length,T)
        f = fft(wave);
        g = abs(f(1:length/2+1)/length);
        figure;
        plot(f)
    end
        

end