@AmusingYeti 

My remark refers to eigenvectors, not Cholesky which anyway does not apply directly to Av = aBv.

It would be interesting to see an example where eigenvalues works only after the mentioned transform.

@MDD 
{3, 4, 11}  contains the indices of the facets in the list L.

@acer 

Thanks, I forgot that simplex is still table-based, but I wonder why ":-" works in Maple 2017 and not in Maple 18 (or maybe <18 ?). I changed this in the code.

@MDD 

OK, I shall include a direct solution in your initial question https://www.mapleprimes.com/questions/223770-How-Can-Remove-The-Redundant-Inequalities

@MDD 

You may try to use the Convex package http://www-home.math.uwo.ca/~mfranz/convex/

This is just like computing m*n using the cartesian product {1,2,...,m} x {1,2,...,n}.

@carriewong 

@carriewong 

I mean that no prerequisites are needed to obtain the desired one-line procedure. You asked for num1dsubspaces and you have it in my answer.

@mmcdara 

1. Forget the square roots.
You have a 3 digits FP unit and two numbers x = 1.41 and y = 1.73.
What is x*y using this machine?

2. On the same machine compute x = sqrt(2) and y = sqrt(3).
GO TO 1.

@Giulianot 

Plotting procedures has some advantages and it's the simplest fix. It does not work for you?

@_Maxim_ 

This has a single index and is of course very useful (and consistent with rtables). My opinion was about two indices L[u,v].

@Alex Bowden 

You refuse to try to understand the problem and see the general picture.
It is far beyond your 4(1+x) and (x+1)(x-1) and you do not have a constructive solution.
There are other solutions to remove the ambiguities. For example, Mathematica accepts 4(x+1) but requires square brackets for functions ( Sin[x],  f[x] ). I prefer Maple's approach. So, let's stop here.

@acer 

This is a parametrized solution too, the difference is that it is not global (it has 4 maps) and uses cartesian coordinates.
The corresponding solution with 2 maps  and spherical coordinates (curves included):

el:=[x=9/2*sin(phi)*cos(theta),y=6*sin(phi)*sin(theta), z=3*cos(phi)]:
sp:=[x=4*sin(phi)*cos(theta),y=4*sin(phi)*sin(theta), z=4*cos(phi)]:
phi1 := arccos(sqrt((63*cos(theta)^2-80)/(63*cos(theta)^2-108))): phi2:=Pi-phi1:
phi3 := arcsin(21/(2*sqrt(140+49*sin(theta)^2))): phi4:=Pi-phi3:
p1:=plot3d(rhs~(el), theta=0..2*Pi,phi=phi1..phi2,color=blue,style=surface):
p2:=plot3d(rhs~(sp), theta=0..2*Pi,phi=phi3..phi4,color=red,style=surface):
p3:=plots:-spacecurve( eval(1.02*rhs~(el),phi=phi1), theta=0..2*Pi ,color=yellow, thickness=2):
p4:=plots:-spacecurve( eval(1.02*rhs~(el),phi=phi2), theta=0..2*Pi ,color=yellow, thickness=2):
plots:-display(p1,p2,p3,p4,scaling=constrained);

@Alex Bowden 

If you really hate the explicit multiplication operator then stick to Math 2D notation and put a space instead of *.

a (b)  or  a b  or (a)  (b)  will be interpreted as a*b.

Note that the "mathematical notation" is (or can be) ambiguous without the context  and a language (such as Maple, Mathematica, C etc) cannot cope with ambiguity.

Working with "generic" functions is essential for Maple and for the user.
The interpretation of   a(x+1)   as   a*(x+1)  would destroy  this functionality.

@Alex Bowden 

@Earl 

phi1 and phi2 are obtained from the intersection:

a,b,c,  r := 9/2,6,3, 4;
el:=[x=a*sin(phi)*cos(theta),y=b*sin(phi)*sin(theta), z=c*cos(phi)];
sp:=[x=r*sin(phi)*cos(theta),y=r*sin(phi)*sin(theta), z=r*cos(phi)];
eq:=eval(x^2+y^2+z^2=r^2,el);
eval(eq, sin(phi)^2=1-cos(phi)^2);
isolate(%, cos(phi)^2);
cos2ph:=simplify(rhs(%));
# plot(cos2ph,theta=0..2*Pi);
phi1:=arccos(sqrt(cos2ph));
phi2:=Pi-phi1;

Please note that for other values  of a,b,c,r  we may need to swap phi,psi or the axes to use the same formulae.
Note also that the shere was not cropped; it was simply obfuscated by the ellipsoid.
Actually phi1, phi2 are needed only for the intersection curve(s).