I'm okay with having the answer returned as a set, as you say Axel it is quite natural, but I do think op(%) is not a very friendly command for beginners (I'm not very comfortable with it myself at times, I do sometimes have to fiddle to get it right).

I see, that looks simpler

... simpler than having to explicitly state "UseAssumptions" before using "assume" -- default behavior being, I infer, "AbsolutelyDoNotUseAssumptions".

... but it returns a set, and since users would often want to assign the solution for later use, it is perhaps less simple than I initially thought:

solve({f(x)=0,x>0},x); whattype(%);
                               {x = 1}


                                 set

I wonder why it is thought useful to explicitly require UseAssumptions

I see, that looks simpler

... simpler than having to explicitly state "UseAssumptions" before using "assume" -- default behavior being, I infer, "AbsolutelyDoNotUseAssumptions".

... but it returns a set, and since users would often want to assign the solution for later use, it is perhaps less simple than I initially thought:

solve({f(x)=0,x>0},x); whattype(%);
                               {x = 1}


                                 set

I wonder why it is thought useful to explicitly require UseAssumptions

I was manipulating symbolic expressions, so digits would not be at fault this time. Maybe something similar to your example with cos(x)^2+sin(x)^2 being simplified to 1.

Thanks for your explanations.

The actual expressions were:

newsys := [-u*(-1+u)*(u*b+u*v-s*v-s*w), u^2*A*a*(a-1)*(A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*d^2+2*A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*d*b+2*A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*d*v+A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*b^2+2*A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*b*v+A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*v^2-A^(-(a-1)^2)*a^(-(a-1)^2)*d^3*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)-2*A^(-(a-1)^2)*a^(-(a-1)^2)*d^2*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*b-2*A^(-(a-1)^2)*a^(-(a-1)^2)*d^2*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*v-A^(-(a-1)^2)*a^(-(a-1)^2)*d*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*b^2-2*A^(-(a-1)^2)*a^(-(a-1)^2)*d*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*b*v-A^(-(a-1)^2)*a^(-(a-1)^2)*d*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*v^2), -u*w^2-u*w*b2+2*u*w*b-u*b1*w-u*b1*b2+u*b1*b+u*b*b2-u*b^2+s*w*v-s*w*u*v+s*w^2-s*w^2*u]:
newsystaylor := [b*u^2-s*v*u-s*w*u, C*u^2, (-b1*b2+b1*b+b*b2-b^2)*u+(-b2+2*b-b1)*u*w+s*w*v+s*w^2]:
 

I was manipulating symbolic expressions, so digits would not be at fault this time. Maybe something similar to your example with cos(x)^2+sin(x)^2 being simplified to 1.

Thanks for your explanations.

The actual expressions were:

newsys := [-u*(-1+u)*(u*b+u*v-s*v-s*w), u^2*A*a*(a-1)*(A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*d^2+2*A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*d*b+2*A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*d*v+A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*b^2+2*A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*b*v+A^(-2*a^2+4*a-1)*a^(-2*(a-1)^2)*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(a*(a-2))*v^2-A^(-(a-1)^2)*a^(-(a-1)^2)*d^3*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)-2*A^(-(a-1)^2)*a^(-(a-1)^2)*d^2*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*b-2*A^(-(a-1)^2)*a^(-(a-1)^2)*d^2*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*v-A^(-(a-1)^2)*a^(-(a-1)^2)*d*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*b^2-2*A^(-(a-1)^2)*a^(-(a-1)^2)*d*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*b*v-A^(-(a-1)^2)*a^(-(a-1)^2)*d*(d^2+2*d*b+2*d*v+b^2+2*b*v+v^2)^(-a)*(d+b+v)^(-1+a^2)*v^2), -u*w^2-u*w*b2+2*u*w*b-u*b1*w-u*b1*b2+u*b1*b+u*b*b2-u*b^2+s*w*v-s*w*u*v+s*w^2-s*w^2*u]:
newsystaylor := [b*u^2-s*v*u-s*w*u, C*u^2, (-b1*b2+b1*b+b*b2-b^2)*u+(-b2+2*b-b1)*u*w+s*w*v+s*w^2]:
 

cat - concatenating expressions

The cat function is commonly used to concatenate strings and names together.

maybe the thing called Hglower0det is inside one of your datafiles. I don't have the datafiles to look inside.

I'm not familiar with this, I'm afraid I can't help any further with this. I would suggest you try to break down the code into smaller pieces and check it step by step.

cat - concatenating expressions

The cat function is commonly used to concatenate strings and names together.

maybe the thing called Hglower0det is inside one of your datafiles. I don't have the datafiles to look inside.

I'm not familiar with this, I'm afraid I can't help any further with this. I would suggest you try to break down the code into smaller pieces and check it step by step.

for this lucid explanation, will keep it in mind.

for this lucid explanation, will keep it in mind.

do you have the data file that readlib is wanting to read set up properly?

do you have the data file that readlib is wanting to read set up properly?

useful to know about the command-completion facility. It's for 2D input, right?

another way is to simply copy-paste from the code posted on mapleprimes, which is what I had done.

useful to know about the command-completion facility. It's for 2D input, right?

another way is to simply copy-paste from the code posted on mapleprimes, which is what I had done.

I ran the code above without any problem, did you refresh your session with restart; at the top? what is your version of Maple? are you using classic or standard? someone else may be able to help -- I'm afraid I don't know what's wrong.

I ran the code above without any problem, did you refresh your session with restart; at the top? what is your version of Maple? are you using classic or standard? someone else may be able to help -- I'm afraid I don't know what's wrong.