# do1
ceps = fft(log(abs(fft(my1))));
plot(log(abs(fft(my1))));
cutout1 = ceps;

for i = 4:(length(my1)-2)
    cutout1(i)=0;
end

result1 = (cutout1(1)^2 + cutout1(2)^2)^(1/2);

# do2
ceps = fft(log(abs(fft(my2))));
plot(log(abs(fft(my2))));
cutout2 = ceps;

for i = 4:(length(my2)-2)
    cutout2(i)=0;
end

result2 = (cutout2(1)^2 + cutout2(2)^2)^(1/2);


// do 3
ceps = fft(log(abs(fft(test))));
plot(log(abs(fft(test))));
cutout3 = ceps;

for i = 4:(length(test)-2)
    cutout3(i)=0;
end

result3 = (cutout3(1)^2 + cutout3(2)^2)^(1/2);