@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?   

@acer Many thanks for the explanation.  The main idea was to have a nicer set of plots without using the options under the right click, that is, not a transient change as you have mentioned. 

@Kitonum Many thanks.   Definitely much better, but somehow maple shows a tabulate?  that counts how many times I issue the command (See figure below).  To get rid of it I had to use a colon at the end of the issued command.

@tomleslie Many thanks for the explanation.   

@Carl Love Thank you ever so much. As usual, I will have to go through the whole script to understand what is going on (For instance, from Types to NLMZR). The figures show exactly what I need, that is, no matter what precision one uses the resulting trajectory is very likely to be chaotic except for the case of the exact arithmetic result. In this case, the concept of a transient is dubious.