fac = 2^(1/12);
t = 1:10000;
c5 = sin(fac^3*440*2*pi*t/10000);
f5 = sin(fac^8*440*2*pi*t/10000);
g5 = sin(fac^10*440*2*pi*t/10000);
a5 = sin(2*440*2*pi*t/10000);
sound ([c5(1:1500) f5(1:2500) c5(1:500) f5(1:1000) c5(1:500) f5(1:1000) g5(1:500) a5(1:3000) f5(1:3000)], 4000)