function    value=TV(A,V,H)
 n=size(A,1);
 n1=n-1;
 Vi=sparse([(1/n)*ones(n1,1)   ((1/n)*ones(n1)-speye(n1)) ]);
value=sparse(Vi*T(A)*Vi');
