F matrices utility function.
replaceBand(a, b, k)
The matrix to be replaced
The matrix with the replacement values
The extend of the replacement: 0 (diagonal only), 1 (diagonal and first extra diagonal), in general an entry is replaced if abs(row(a) - col(a)) <= k
A matrix (same size as a)