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: Wrapper_BFR_iHARPERELLA.m
    path of file: ./repos/sepia/wrapper/WrapperBFR
    
    the code of the file until where you have to start completion: %% [RDF] = Wrapper_BFR_iHARPERELLA(totalField,mask,matrixSize,voxelSize,algorParam, headerAndExtraData)
%
% Input
% --------------
% totalField    : total field map (background + tissue fields), in Hz
% mask          : signal mask
% matrixSize    : size of the input image
% voxelSize     : spatial resolution of each dimension of the data, in mm
% algorParam    : structure cont
    
function to access iHARPERELLA for SEPIA
%
% Kwok-shing Chan @ DCCN
% k.chan@donders.ru.nl
% Date created: 8 March 2020 (v0.8.0)
% Date last modified:
%
%
function [RDF] = Wrapper_BFR_iHARPERELLA(totalField,mask,matrixSize,voxelSize,algorParam, headerAndExtraData)
sepia_universal_variables;

% get algorithm parameters
algorParam  = check_and_set_algorithm_default(algorParam);
method      = algorParam.bfr.method;
iteration  	= algorParam.bfr.iteration;

% add path
sepia_addpath('STISuite');

%% Display algorithm parameters
disp('The following parameter is being used...');
disp(['Maximum iterations = ' num2str(iteration)]);

%% main
RDF = iHARPERELLA(totalField, mask,'voxelsize',voxelSize,'niter',iteration);
       
end