Is the following a bug? I am using Maple 2019  64 bit with latest Physics package 357 on windows 10.

Download bug2.mw

This

expr:='1/r^2*diff(r^2*diff(v(r),r),r)'

gives

And when evaluated second time it gives

expr

The question is, how to put expr back to its original form shown at the top?

I tried simplify, combine, etc.. but nothing puts back to same form.

The reason I am asking, is that I need put more complicated expression in form that resembles when I am looking at in textbook. For an example, the Laplacian

VectorCalculus:-Laplacian(u(r,theta,phi),'spherical'[r,theta,phi])

But the above should be the same as

So I figured if I can rewrite each term to look like in the shorter version above (complete derivative version), it will be easier for me to compare what Maple gives and the text book shows,

ps. I found that the following command in Maples gives a little clearer Laplacian, but it is still not as simplified as the book but it is better than using VectorCalculus:-Laplacian

Physics[Vectors]:-Laplacian(u(r,theta,phi));

At least this has each term a little more seprated.

When typing 

z:=exp(I*2*Pi/3);
convert(z,'sincos')

Maple evaluates the intermediate result which is cos(2*Pi/3)+I*sin(2*Pi/3) and  gives

Is there a way to tell it not to do this? I'd like to see the result as when typing

'cos(2*Pi/3)+I*sin(2*Pi/3)'

Is there an option or method to tell Maple not to immediate evaluation in the above? it can do evaluate next time the expression is used.

Should Maple handle this error internally and may be give no solution if it can't solve it instead of this  error?

When setting boundary condition to zero, maple gives error below. Heat PDE in a sphere. No angle dependency. Only the radial part.

unassign('r,u,t');
pde:=diff(u(r,t),t)= 1/r*diff(r*u(r,t),r$2); #Laplacian in spherical
ic:=u(r,0)=1;
bc := u(1,t) =0;
pdsolve([pde,ic,bc],u(r,t),HINT =boundedseries(r=0)) assuming t>0

Error, (in assuming) when calling 'ln'. Received: 'when calling 'ln'. Received: 'numeric exception: division by zero''

Adding assumptions such as 0<r,r<1 did not help.

Physics version 348. Maple 2019 on windows 10. Is there something I am doing wrong? 

Why does Maple display eval(diff(f(r),r),r=C)=0;  differently when C is a number?

unassign('f,C,r');
eval(diff(f(r),r),r=C)=0;
eval(diff(f(r),r),r=5)=0;

I wanted the first one to display the same as the second one. First one looks confusing. it looks like taking derivative w.r.t to C. Is there a way to make the first one look like: