This is certainly a bug. The number following "integer : " is varying with the exponent; try exp(-2*x) or exp(-3*x) in the ic setting.

If you obtain the general solution of pde (w/o ic), you can see that y=0 may be problematic (just like x=0, by symmetry). So let's replace it by a symbol in ic:

ic := u(x,y__0) = exp(-x);
sol := pdsolve([pde,ic]);
pdetest(sol,[pde,ic]);

Yes, MapleSim can do such simulations. Some aspects are covered by examples in our Model Gallery. Please see the Robotics category in particular.

No, they don't require entering complex formulae for drawing. The animations are produced as (part of the) simulation results.

I don't think we have a direct link to 3D printing, but we can import STL files (in the base product MapleSim), and many more formats via our CAD Toolbox.

@kambiz1199 Very old-fashioned, but if you insist...

M:=Matrix(6,2):
for k from 1 to 6 do x:=(k+1.)/2; M[k,1]:=x: M[k,2]:=f(x) end do:
M;

Kitonum's suggestion is much better, though. Make sure you get the function definition right, as he demonstrated.

What other representation would you prefer? Please specify.

The paper ends with this remark:

The computations associated with the examples in this Letter were performed using Maple 10.

So it's probably best if you ask the authors for their worksheet(s).

@michalkvasnicka The installer is designed to be generic. If it came with libraries for specific distributions, users of others would run into problems.

@ahmadtalaei I was unable to obtain any solution - sorry. One of the attempts exhausted my computer's memory of 20 GB RAM.

BTW: the multiplicative separation ansatz T(t)*R(r)*Phi(phi) can be written a bit shorter: HINT=`*`

Please upload your worksheet to MaplePrimes. I followed the SE link and then the Dropbox link - which doesn't open for me.

@Carl Love Ctrl + mousewheel supports some more factors (20 if I counted correctly).

@rahinui Looks like a glitch in odetest of Maple 2018. In Maple 2019, it returns all the expected zeros.

The problem seems to be related to the repesentation in hyperbolic functions, as you can see from the exponential representation obtained via

sol := dsolve(laplace union BCs);
odetest(sol, laplace union BCs);

which succeeds in 2018 as well.. You can run sol := convert(sol,trigh) afterwards.

@a_simsim If you enter ?MapleSim,LinkModel,Simulate, you will see the datapoint option. Also take a look at the last example on that page.

I have long-running calculations all the time, and routinely lock my laptop. Never had any problem with that. Maple always remained responsive.

So it would be nice if you provide more details (Maple version, OS version, sample worksheet(s), configuration details). Please send them to support@maplesoft.com for analysis.

@dharr If you replace pi by Pi, you will get a slightly shorter result.

@J F Ogilvie The 2018.2 update can also be applied to an installed Maple 2018.0, as always.

@tomleslie Glad to hear the problem is solved. I'm not a developer, but AFAIK, the Physics package has very high (or specialized) demands on 2D math (e.g. bra-ket notation), so it's modifying some parts of the typesetting library. I'm not surprised that side effects may occur.