function Y = matvec(L,U,C,X); 

Y = -U\(L\(C*X));
