function result=ste(name)
 [x,f,n] = wavread(name); 
 blocks = floor(length(x)/64);
 result = zeros(1,blocks);
 for i=0:blocks-1
  for j=1:64
   result(i+1) = result(i+1) + x(i*64+j)*x(i*64+j);
  endfor
 endfor
endfunction
