function fracsound1213(mod,filen,cycles)
%generates a fractal sound which is continually ascending 
%although cyclic in mod steps (try mod=13,12,5) asoundtst(13,0,'')
%saves file if filen~='' 
%mod=13;
div=12; %semitones in an octave
samp=3*mod*1000;
quo=samp/mod;
t=1/2^5;
be=2^(mod/div);
b=zeros(1,samp*cycles);
for g=1:cycles
    for h=(1:mod)
        for i=(1:quo)
            pos=i+(h-1)*quo+(g-1)*samp;
            for j=(0:20)
                b(pos)=b(pos)+cos((be^j)*(t*2^(h/div))*i);
            end	
        end
    end
end
mx=max(b);
mn=min(b);
b=b/(mx-mn);
%plot(b);
sound(b,3*4000); 
%player=audioplayer(b,3*4000);%alternative playing method
%play(player);
if ~isempty(filen)
    wavwrite(b,filen);
end