find the indexes of constant colums
[xlag]=search_cte(x,n,init)
* x = a (n x k) real matrix
* indcte = a (p x 1) vector of integers
search_cte([grand(15,2,'nor',0,1) ones(15,1)]) // Returns 3.