function res = T(nn,kk,qq,nb)

res = [cos(nb*qq*3^(nn-kk)),-sin(nb*qq*3^(nn-kk));sin(nb*qq*3^(nn-kk)), cos(nb*qq*3^(nn-kk))];


