// Updated(18-7-07)
// 7.1

H = trfu(1,[1 -1 0],-1);
evans(H)