public IOperator identity(int size) { TridiagonalOperator I = new TridiagonalOperator(new Vector(size - 1, 0.0), // lower diagonal new Vector(size, 1.0), // diagonal new Vector(size - 1, 0.0)); // upper diagonal return I; }