clear aa;
j=0;
k0=5.4;
for m=exp1:exp2-1;
   jj=inter-1;
   for n=0:jj;
      a=2^(m+n/inter);
      j=j+1;
      aa(j)=a;
   end;
end;
a=2^exp2;
aa(j+1)=a;
omega0=1/2*(k0./aa+sqrt(2+k0*k0)./aa);
period=1./omega0*2*pi;
aa=aa';
period=period';



  



