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: ft_specest_neuvar.m
    path of file: ./repos/spm12/external/fieldtrip/specest
    
    the code of the file until where you have to start completion: function [spectrum, freqoi] = ft_specest_neuvar(dat, time, varargin)

% FT_SPECEST_NEUVAR computes a time-domain estimation of overall signal 
% power, having compensated for the 1/f distribution of spectral content.
%
% Use as
%   [spectrum,ntaper,freqoi] = ft_specest_neuvar(dat,time...)
% where
%   dat        = matrix of chan*sample
%   time       = vector, containing time in seconds for each sample
%   neuvar     = matrix of chan*neuvar
%
% Optional arguments should be specified in key-value pairs and can include
%   order      = number, the order of differentation for compensating for the 1/f (default: 1)
%   pad        = number, total length of data after zero padding (in seconds)
%   padtype    = string, indicating type of padding to be used (see ft_preproc_padding, default: 0)
%   verbose    = output progress to console (0 or 1, default 1)
%
% See also FT_FREQANALYSIS, FT_SPECEST_MTMFFT, FT_SPECEST_MTMCONVOL, FT_SPECEST_TFR, FT_SPECEST_HILBERT, FT_SPECEST_WAVELET

    
function [spectrum, freqoi] = ft_specest_neuvar(dat, time, varargin)

% FT_SPECEST_NEUVAR computes a time-domain estimation of overall signal 
% power, having compensated for the 1/f distribution of spectral content.
%
% Use as
%   [spectrum,ntaper,freqoi] = ft_specest_neuvar(dat,time...)
% where
%   dat        = matrix of chan*sample
%   time       = vector, containing time in seconds for each sample
%   neuvar     = matrix of chan*neuvar
%
% Optional arguments should be specified in key-value pairs and can include
%   order      = number, the order of differentation for compensating for the 1/f (default: 1)
%   pad        = number, total length of data after zero padding (in seconds)
%   padtype    = string, indicating type of padding to be used (see ft_preproc_padding, default: 0)
%   verbose    = output progress to console (0 or 1, default 1)
%
% See also FT_FREQANALYSIS, FT_SPECEST_MTMFFT, FT_SPECEST_MTMCONVOL, FT_SPECEST_TFR, FT_SPECEST_HILBERT, FT_SPECEST_WAVELET

% Copyright (C) 2016, Arjen Stolk
%
% This file is part of FieldTrip, see http://www.fieldtriptoolbox.org
% for the documentation and details.
%
%    FieldTrip is free software: you can redistribute it and/or modify
%    it under the terms of the GNU General Public License as published by
%    the Free Software Foundation, either version 3 of the License, or
%    (at your option) any later version.
%
%    FieldTrip is distributed in the hope that it will be useful,
%    but WITHOUT ANY WARRANTY; without even the implied warranty of
%    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%    GNU General Public License for more details.
%
%    You should have received a copy of the GNU General Public License
%    along with FieldTrip. If not, see <http://www.gnu.org/licenses/>.
%
% $Id$


% get the optional input arguments
order     = ft_getopt(varargin, 'order', 1);
pad       = ft_getopt(varargin, 'pad');
padtype   = ft_getopt(varargin, 'padtype', 'zero');
fbopt     = ft_getopt(varargin, 'feedback');
verbose   = ft_getopt(varargin, 'verbose', true);

if isempty(fbopt)
  fbopt.i = 1;
  fbopt.n = 1;
end