@C_R  You can easily plot a cube with this method:

plots:-implicitplot3d(max(-x,x-1,-y,y-1,-z,z-1), x=-0.5..1.5, y=-0.5..1.5, z=-0.5..1.5, style=surface, axes=normal, grid=[100,100,100]);

@C_R  You wrote "...Can you give me hint what  max  does with implict geometries?...".

The intersection of 2 of these tori (including the surfaces themselves) is given by the system of inequalities  {T1<=0, T2<=0}  that is equivalent to 1 inequality  max(T1,T2) <=0 . The equality  max(T1,T2)=0 defines the surface of this solid.

@vs140580  Remove the  print  command from the body of your procedures:

trying_new.mw

@vs140580  I don't understand why you can't do this:

op(0,RandomTools:-Generate(choose({A(Graphln, k),B(Graphln, k),C(Graphln, k)})));

@Carl Love  Thanks, I've edited my answer.

@nm  Thanks a lot for the answer!

@acer  Thanks a lot for the detailed answer!

@mmcdara  I agree, but this is not a fundamental difference. We can add the  collect  command:

x^(n-1)*(collect(expand((n*x^n - 2*n*x^(n - 1) + x^n)/x^(n-1)), x));

@lcz If you are interested in arbitrary set partitions, you can also take a look at my  post  https://mapleprimes.com/posts/200776-Partition-Of-A-Set-Into-Parts-Of-Any-Sizes .  It was written long ago, when Maple didn't have the  Iterator package yet.

@Deeshani 

 

Download metriccalculaions_new_(12).mw

@dharr  OK. Vote up.

@ecterrab In fact, the exact solution is easily expressed explicitly (see below). But a natural question arises: why does this solution not continue to the right from  x~2.55 , as shown by the numerical solution?

restart;
ode:=diff(y(x),x,x)+2*diff(y(x),x)^2+8*y(x)=0:
ics:=y(0)=0, D(y)(0)=1:
dsolve({ode, ics}, y(x), implicit);
value(%);
isolate(%, y(x));
sol := simplify(%);
odetest(sol, [ode, ics]);

@Thomas Dean  Using the P procedure, you can plot the body of revolution for any value of the parameter  a  (this is the angle of rotation).  P (2*Pi+Pi/2)  will be the whole body:

P(5*Pi/2);

@Axel Vogt  A good idea! This can be done even without using the plottools package:

restart;
plot([erf(y),y, y=-3..3]);

@Christian Wolinski 

It takes approximately 10 sec in Maple 2018.2.