@John Fredsted Many thanks.   I have just copied and pasted the example you've just send and myTranspose returned the answer for the variable expr (defined earlier) and not for M.  What am I missing?  

@John Fredsted Many thanks.   Do I always have to define the matrix to be defined as expr?   

@Carl Love Many thanks. a)  I had that and the solution figured out before you sent the reply, but it is good to know that I got it right.  b) I have completely forgotten that.  Many thanks.  

@Carl Love Dear Carl,  thank you very much (again).

If you have the time could you help me with the following issues?

a) OK.   I did as suggested and some of the results I was expecting to see showed up.  Example: eval(CDF(X2,v),[v=0.64,m=3]).  However when I tried eval(Quantile(X2,v),[v=0.95,m=3]), Maple 2017 returns FAIL.  That is a bit discouraging since I need to increase m up to 50 and calculate the quantiles. I hope I am missing something here.

b) Could you explain what does $ 2 do?   

Again many thanks.

@Carl Love Many thanks.   As you know I work with sets of polynomials.  For each of them I need to make a call to Groebner independently. If I got it right, I can send each call to a core (24 core each time) and them save the results, but unfortunately I don't know how. 

@Carl Love models=subsets. I need all the subsets of [f(x,y,z),g(x,y,z),h(x,y,z)] with f(x,y,z) containing that specific monomial. Please note that MyTest will remove some of the subsets. I don't know the package but it sounds worth trying. The functions f,g and h are full quadratic polynomials in x,y and z (V as used in the procedure). Many thanks.

@Carl Love Many thanks.  As usual I am impressed with the solutions.  Many thanks to @Kitonum too.

I believe that you know, by now, what kind of problem I am dealing with, that is, a set of  full quadratic models with coefficients thetas (30 in total), model=[f(x,y,z),g(x,y,z),h(x,y,z)]. Per chance, would you know to create all the  models with a specific monomial in it (I hope I have used the right term)? For instance, the model used in this thread has the monomial x*y*theta[5] in it.  Taking as the starting point the full quadratic model, how many models can one find with the aforementioned monomial?   Bulding all of them does not seem a good idea, 2ˆ30-1.  

@Kitonum Many thanks.  It works.

@Carl Love  Here is the error.  

Error, cyclic dependency detected in parameter V::(depends(And(list(name), satisfies(proc (V) options operator, arrow; nops(V) = nops(F) end proc)))) in procedure MyTest

@Kitonum Many many thanks.  I completely forgot about depends.  

@Carl Love Many thanks for the procedure. Unfortunately it returns a cyclic error on Maple 14.   What is the precisely mathematicsal term to be used, please? (I am not trying to be funny or else)

@acer Many many thanks.   I will try to see if I can understand what Maplets is all about.  

@acer Many many thanks.  Not do I have a progressslider but also to save my own packages. 

I would appreciate if you could find a way to get it working in Maple 14.  I need to know where my proc got stuck after the 30 GB of ram ran out.   

@acer  Many thanks.   It works just fine.   Could that be part of a package?   A similar function that works on Maple 14 would be great. Many thanks.

@emendes I must be dumb! I could not reproduce the example given on examples,ProgrammaticContentGeneration on a clean worksheet. I do not know where the icon just before ProgressBar comes from.  I am completely lost.