@acer Many thanks.   As you said it is pretty straightforward but I still had to pass the result of map(int,eq,t) to Parts to get the result.  

@Christian Wolinski Would you know how to use Rule from Student[Calculus1] to do the same calculations as above? I could not figure out how to do it.  Many thanks.  

@Christian Wolinski Many thanks.  I did not know about the student[intparts].  

@acer Many thanks.  I am not sure what form I need. The only thing I have in mind is to get all indets that have solution sorted out and the rest of indets in a form that is as simple as possible.  Any special advice on that?   I know that I need to carefully go through the commands in the worksheets you sent in order to understand what they do.   

@acer Thank you ever so much. I will have to go through the worksheets step by step to understand what you did.  Meanwhile, would you mind looking at a much more complicated example, please?  example1_eqns.mw

@acer Many thanks.  I have uploaded a worksheet as you requested.
 

Download example_eqns.mw
 

@Christian Wolinski thanks for pointing out eliminate and a couple of new alternatives that I wasn't aware of.  The downside of the solution is that you entered with the values of theta[2,1] and theta[3,0],  

@acer Many thanks

eqns := {-theta[1, 2]*theta[2, 6]-theta[1, 3]*theta[3, 6]+1 = 0, theta[1, 2]*theta[2, 6]+theta[1, 3]*theta[3, 6]-1 = 0, -theta[1, 2]*theta[2, 6]*theta[3, 6]-theta[1, 3]*theta[3, 6]^2+1 = 0, -theta[1, 1]*theta[1, 2]*theta[2, 2]*theta[2, 6]*theta[3, 6]-theta[1, 1]*theta[1, 3]*theta[2, 2]*theta[3, 6]^2+1 = 0, -theta[1, 3]*theta[2, 2]+theta[1, 3]*theta[3, 3]+alpha+rho = 0, theta[1, 0]*theta[1, 3]*theta[2, 2]^2*theta[3, 3]-theta[1, 0]*theta[1, 3]*theta[2, 2]*theta[3, 3]^2-alpha^2*beta-alpha*beta*rho = 0, theta[1, 1]*theta[1, 2]*theta[2, 6]+theta[1, 1]*theta[1, 3]*theta[3, 6]+theta[1, 2]*theta[2, 2]*theta[2, 6]+theta[1, 3]*theta[2, 2]*theta[3, 6]-alpha = 0, theta[1, 1]*theta[1, 2]*theta[2, 6]*theta[3, 6]+theta[1, 1]*theta[1, 3]*theta[3, 6]^2+theta[1, 2]*theta[2, 2]*theta[2, 6]*theta[3, 6]+theta[1, 3]*theta[2, 2]*theta[3, 6]^2-alpha = 0, theta[1, 1]*theta[1, 3]*theta[2, 2]-theta[1, 1]*theta[1, 3]*theta[3, 3]+theta[1, 3]*theta[2, 2]^2-theta[1, 3]*theta[3, 3]^2-alpha^2+rho^2 = 0, -theta[1, 1]*theta[1, 2]*theta[2, 6]-theta[1, 1]*theta[1, 3]*theta[3, 6]-theta[1, 2]*theta[2, 2]*theta[2, 6]-theta[1, 2]*theta[2, 6]*theta[3, 3]-2*theta[1, 3]*theta[3, 3]*theta[3, 6]-2*rho = 0, -theta[1, 0]*theta[1, 2]*theta[2, 2]*theta[2, 6]*theta[3, 6]-theta[1, 0]*theta[1, 3]*theta[2, 2]*theta[3, 6]^2-theta[1, 1]*theta[1, 2]*theta[2, 2]*theta[2, 6]*theta[3, 3]+theta[1, 1]*theta[1, 3]*theta[2, 2]^2*theta[3, 6]-2*theta[1, 1]*theta[1, 3]*theta[2, 2]*theta[3, 3]*theta[3, 6]-alpha-2*rho = 0, theta[1, 1]*theta[1, 2]*theta[2, 6]*theta[3, 3]-theta[1, 1]*theta[1, 3]*theta[2, 2]*theta[3, 6]+2*theta[1, 1]*theta[1, 3]*theta[3, 3]*theta[3, 6]+theta[1, 2]*theta[2, 2]*theta[2, 6]*theta[3, 3]-theta[1, 3]*theta[2, 2]^2*theta[3, 6]+2*theta[1, 3]*theta[2, 2]*theta[3, 3]*theta[3, 6]+alpha^2+2*alpha*rho = 0, -theta[1, 0]*theta[1, 2]*theta[2, 2]*theta[2, 6]*theta[3, 3]+theta[1, 0]*theta[1, 3]*theta[2, 2]^2*theta[3, 6]-2*theta[1, 0]*theta[1, 3]*theta[2, 2]*theta[3, 3]*theta[3, 6]+theta[1, 1]*theta[1, 3]*theta[2, 2]^2*theta[3, 3]-theta[1, 1]*theta[1, 3]*theta[2, 2]*theta[3, 3]^2+alpha*beta+alpha*rho+rho^2 = 0, -theta[1, 0]*theta[1, 2]*theta[2, 2]*theta[2, 6]-theta[1, 0]*theta[1, 3]*theta[2, 2]*theta[3, 6]-theta[1, 1]*theta[1, 3]*theta[2, 2]^2+theta[1, 1]*theta[1, 3]*theta[3, 3]^2-theta[1, 3]*theta[2, 2]^2*theta[3, 3]+theta[1, 3]*theta[2, 2]*theta[3, 3]^2-alpha^2*rho-alpha*rho^2+alpha+beta+rho = 0}

and

fc := {theta[1, 0], theta[1, 1], theta[1, 2], theta[1, 3], theta[2, 2], theta[2, 6], theta[3, 3], theta[3, 6]}

I found out how to add the second module. Save allows me to add more modules.

@acer Many thanks. It worked as you said (including loading automatically).   I have followed all the steps in the help file but to no avail. 

As I said I have two modules and if I got it right Save only saves one module.  How can I append the second module to the same mla file?

Yes, it is OSX. The files were created using the line codes of the example in the help files.  Maple put .mapleinit in /Users/eduardo . What bugs me is that Maple help uses around 15 lines of code to do that same as you did in 5 lines and does not work.  

Regarding adding help files,  what you recommend?  

Many thanks.

@Christian Wolinski  Many thanks.   Now it returns 14 lines.  I will apply the function to other examples to see how it goes.

Thanks but something is wrong.   Line 7 reads [0, 0, 0, 0, 1, 0, 0, 0, 0] but it should be [0, 0, 0, 0, 1, 0, 0, 0, 1]. There is one product of two indeterminates in the fourth element of eq and no summation which means that only one line should be added to the matrix with two ones (since we have a product of two indeterminates).  I guess when there is just a term and the term is a product, the procedure fails.   That seems the case since the number of lines of resulting matrix is 15 instead of 14.  Lines 7 and 8 should be one instead. 

@Kitonum Many many thanks.  1) That seems to be the case and 2) It worked as a sequence that can be easily be converted to a list.  

@Kitonum Many thanks.   Why do I need to use zip in this case?   I have used a similar command in one of my functions, ff:=(l::list)-> eval(f,vars~=l)), that works in Maple 14 without a problem.  (f and vars are inputs to the function).

As for the second option, it does not work when I try f([x,y,z]) and then f(f([x,y,z])).  

@Carl Love Many thanks.  I retyped the equation twice, but I didn't pay attention to the other command line.  

There is a singularity at 26.3 (p[w] tends to zero)  Any suggestions to check if there is a solution beyond that point?