You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: nc_temporal_energyratio.m
    path of file: ./repos/gift/GroupICAT/icatb/toolbox/noisecloud/scripts
    
    the code of the file until where you have to start completion: function feature = nc_temporal_energyratio(T);
      global TR;
      % Calculate FFT of timeseries;
     freq_data = abs(fft(T,length(T)));
      % FEATURE 112: Ratio of the integral of spectral power below 0.10 Hz to;
     % the integral of power between 0.15 and 0.25 Hz;
     % Integral just means the area under the curve I think :O);
     ntime = length(T);
     nfft = ntime/2;
      % Set spectral power min, middle, and high values;
     spectral_Hz_min = 0.10;
     spectral_Hz_mid = 0.15;

    
function feature = nc_temporal_energyratio(T);
      global TR;
      % Calculate FFT of timeseries;
     freq_data = abs(fft(T,length(T)));
      % FEATURE 112: Ratio of the integral of spectral power below 0.10 Hz to;
     % the integral of power between 0.15 and 0.25 Hz;
     % Integral just means the area under the curve I think :O);
     ntime = length(T);
     nfft = ntime/2;
      % Set spectral power min, middle, and high values;
     spectral_Hz_min = 0.10;
     spectral_Hz_mid = 0.15;
     spectral_Hz_max = 0.25;
      spectral_min_index = round(nfft*spectral_Hz_min)/(1/(2*TR));
     spectral_mid_index = round(nfft*spectral_Hz_mid)/(1/(2*TR));
         spectral_max_index = round(nfft*spectral_Hz_max)/(1/(2*TR));
          % compute energy in the lower range bucket (below 0.10);
     low_bucket_energy = 0;
     for column = 1:spectral_min_index-1;
         low_bucket_energy = low_bucket_energy + freq_data(column)^2;
     end