function f = Kmult(T,x,K)
K = toeplitz([2 -1 zeros(1,199)]);
f = -K*x;