Takes the alpha matrix and the j-th eigenvalue lambdaj as input and computes
at each interface (i = 1,...,n+1). All other matrix entries (row # j) are set to zero. The resulting output is the matrix teta.
Release Notes: