program define var_nr_irf_bands, rclass
version 16.0

	syntax , VARname(string) OPTname(string) OUTname(string) [statamatrix(string)]

	mata: `outname' = irf_bands_funct(`varname',`optname')

end