The authors prove the stepwise stability for a finite difference schem
e for the heat equation with an integral constraint. The resulting mat
rix is nonsymmetric and does not have the usual band structure. The pr
oof is based on the method of matrix analysis. The eigenvalues of seve
ral matrices are found explicitly or their location described precisel
y. This method relies upon the relationship of the characteristic poly
nomials of these matrices with orthogonal polynomials.