function Y = matvect(L,U,C,X) 

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