@Carl Love It was fixed between Maple 2015.2 and Maple 2016.2.

@Arastas The trailing tilde makes the command map across entries of a Matrix, or list, etc, in an elementwise manner.

Some commands automatically get applied across elements of such structures, but most do not.

@etian2 It works for me. In 2D Input mode I type in,

  a := 5

but before I press <Enter> I select the "a", make it bold, then right-click to make it an atomic identifier. Then I press <Enter>.

Then in a new statement I type "a", repeat the process to make is bold and an atomic identifier. And after I press <Enter> the evaluated output shows as 5.

Download bold_assignment.mw

Why do you think that your procedure is thread-safe, given that it calls signum?

I don't see signum listed on the Help page for Topic threadsafe.

I don't know whether that is part cause of your problem.

It does not make sense to say that you cannot upload a Maple 11 worksheet because "Too complex to insert the Maple 11 version because I have  a complicated function  and further on".

The integral you claimed to work in Maple 11 is short. It either works, by itself and simple, or it doesn't.

@Preben Alsholm I get the following:

restart;

kernelopts(version);

    Maple 11.02, X86 64 LINUX, Nov 9 2007 Build ID 330022

with(IntegrationTools):
with(DEtools):

fx:=exp(-sqrt(x^2 + c)):

integral:=int(fx,x) assuming positive;

      int(exp(-(x^2+c)^(1/2)),x)

value(integral);

      int(exp(-(x^2+c)^(1/2)),x)

int(fx,x=xl..xu,method = _DEFAULT) assuming (x^2+c)>0,x>0,c>0;

              2     1/2
  int(exp(-(x  + c)   ), x = xl .. xu, method = _DEFAULT)

int(fx,x=xl..xu) assuming (x^2+c)>0,x>0,c>0;

   int(exp(-(x^2+c)^(1/2)),x = xl .. xu)

It is poor etiquette and unhelpful to show only an image, instead of actual code to reproduce.

You could lprint the expression, and show us that as plaintext, or you could upload and attach a worksheet using the green up-arrow in the Mapleprimes editor.

Do you work with interface(typesetting) as extended, or standard? It can make a difference, for some workarounds.

@weidade37211 I've added a few faster approaches to my Question, now edited.

One key problem in your original was coding it so as to perform the symbolic integration  int(f(x),x=0..y)  where y is a nonnumeric symbolic name. Your 2D Input of the integral is pretty-printed in black in your original, denoting the active int command.

Using a restricted range can also help the numeric root-finder (ie. fsolve).

@nmacsai I also added an alternative workaround, using a wrapping call to plots:-display.

The approach of altering the Histogram procedure is easier only if you put the preamble code in an initialization file (which would be done just once). After that the Histogram command handles your example directly.

But if you don't mind the little extra bit of typing then the plots:-display approach is also an easy way to handle your example.

@mmcdara Your comment doesn't pertain to the OP's example. And so I removed the additional, concocted example to which you objected -- which was only a sidebar here and didn't demonstrate the OP's problem.

You might discuss the new format=stacked functionality elsewhere. That's not the main subject of this Question thread.

@Guimzo Have you considered describing exactly what you want?

@David Sycamore 

Here is a worksheet. Does it run ok for you (without your making any edits to it)?

If not then you should tell us what version of Maple you are using. You should also upload and attach you'r problematic worksheet using the green up-arrow in the Mapleprimes editor.

So far you've only told us that you cannot get what you want, without actually providing us with exactly whatever you've been trying.

Download ptplot.mw

@Preben Alsholm I wasn't paying close attention, sorry -- it could have been as short as 10 minutes, and I doubt it took longer than 3/4 of an hour. I did it on Maple 2022.0 on Linux.

@bobtom Do you mean output like this?

[seq([seq([a[i,j],b[(9/2-j/2)*j-i+1]],
          i=5-j..1,-1)],
     j=1..4)];

   [[[a[4, 1], b[1]], [a[3, 1], b[2]], [a[2, 1], b[3]], [a[1, 1], b[4]]],

    [[a[3, 2], b[5]], [a[2, 2], b[6]], [a[1, 2], b[7]]],

    [[a[2, 3], b[8]], [a[1, 3], b[9]]], [[a[1, 4], b[10]]]]

Another way to consider that is as the following, when N=4,

[seq([seq([a[i,j],b[j*(2*N-j+1)/2-i+1]],
          i=N-j+1..1,-1)],
     j=1..N)];