@kevinmaple The arrow notation -> is an abbreviation allowed for procedures (usually one-liners) that contain only an expression, an expression sequence, or an if statement. So

CountUnique:= (L::list)-> nops({L[]}):

is essentially equivalent to

CountUnique:= proc(L::list) return nops({L[]}) end proc:

See ?operators,functional

It's possible to do what you want using the parse statement. However, there's almost always a better way. What's your ultimate goal? to generate a Cartesian product of sets or lists? or some other combinatorial object? These things can be done without using loops nested to arbitrary depth. Nonetheless, in the Answer below, I show how to do what you asked for.

@Preben Alsholm Yes, thanks, that's what I meant. I corrected the Answer.

@GPY You need to minimize the 2-norm, which is equivalent to minimizing its square. Try this:

Minimize(LinearAlgebra:-Norm(minmat, 2)^2);

This should resolve the absolute-value error.

@ger89 The second argument to pdsolve needs to be w(x,t), not w(u,t).

@mzaman Nobody wants to retype your code. Instead of a screen shot, please post plaintext or upload a worksheet. You can upload worksheets and a few other file types by using the green uparrow that is the last item on the second row of the toolbar in the MaplePrimes editor.

@Kitonum 

Note, however, the significance of round-off. You had to change the OP's 0.36 to 0.3603 in order to get equality. I recommend that the OP work in exact rational arithmetic for such problems. That can be done by not using any decimal points in the input.

@lt I still don't see the problem that is the subject of this Question. The plot shown in your worksheet clearly shows that your textplot3d commands are working.

There is no need for you to direct my attention to your other Questions. I read and consider everything posted on MaplePrimes. 

@UndergradMaths 

Right clicking on the output of your equation brings up a context menu. Select "Solve DE."

If you want to type a command, type dsolve(%);.

@lt Would you please post a worksheet where the above command doesn't work? Your syntax looks fine to me.

What do you mean by "a label coordinate"?

Are the circles in the same plane with one inside the other? On first thought, that seems an easier problem than two arbitrary circles. 

@Axel Vogt

I'm just writing this to clarify Axel's advice: If you want to avoid this anomaly, use 1/10 instead of 0.1. The floating-point anomaly that you point out is just a standard example, not an bug.

@Kitonum I agree with Markiyan. What relation is there between your workaround plot and the matter discussed in the paper which the OP referenced?

@Kitonum Here's another improvement. This avoids the once-per-rotation stutter caused by 0 and 2*Pi being actually the same frame. We start the rotation a little higher than 0 based on the frame count. This improvement can be applied to any continuous-play rotating animation:

Frames:= 54:
plots:-animate(
     plots:-tubeplot,
     [[cos(t+phi), sin(t+phi), t], t= 0..6*Pi, radius= 1/3, numpoints=200],
     phi= (2*Pi)/Frames..2*Pi, frames= Frames,
     orientation= [40,70], shading= z, axes= box, style= patchnogrid,
     scaling= constrained
);