function[vec] = M(n,k,q)

na = 3;
nb=1;

    
    
C = [1,0 ; 0,na/nb];

D = [1,0;0,nb/na];
if n == 1
    vec = T(0,k,q,nb)*C*T(0,k,q,nb)*D*T(0,k,q,nb);
    return
else
    vec = M(n-1,k,q)*T(n-1,k,q,nb)*M(n-1,k,q);
end
end





