`diagonal_covariance_i` computes the approximate diagonal covariance matrix of the functional data for lag windows defined by i.
diagonal_covariance_i(f_data, i)
A 2-dimensional array, encoding the covariance matrix of the functional data for lag windows defined by i.
the functional data matrix with observed functions in the columns
the index in 1:T that we are computing the covariance for