@Carl Love 

Thanks. Your functions worked very well. Verified on Maple 2019.1

@tomleslie 

Nice. But I do not think you need max there.

I rewrote your function a little to make it easier for me to understand (I do not like the style of using () over many lines, I find it hard for me to read, I've also seen others use the style you are using here).  (ps. I understand code style is a personal thing. All what I am saying is that I find this harder for me to read, that is all).

DETools:-difforder worked without using max in the following

restart;
get_order:=proc(expr,fun,var)
  local the_list,i;
  the_list := convert( indets(expr, function),list);
  PDETools:-difforder(select( i-> has(i, diff(fun,var)),the_list))
end proc:

Now 

f:= 2*y*(diff(x(t), t))^2+3*(diff(x(t), t$3))-3*x*(diff(y(t), t));
get_order(f, y(t),t);
get_order(f, x(t),t);

Returned 1 and 3 also.

@Carl Love 

Thanks. I verified your solution also works very well on Maple 2019.1

@acer 

Thanks. I verified your code works with no problem on 2019.1

@Joe Riel 

Thanks, but it did not work on the following example:

Do I need something else other than type the command you showed? Using Maple 2019.1 on windows.

Download bug_june_6_2019_3_45_pm.mw

@Carl Love 

Thanks for the reply. But the latex command should internally handle these issues.

What if you have 10,000 equations to process? Is a user supposed to check each time for something and make a workaround to make latex work?

Latex can simply use \text{All Solutions} or anything else it needs to make the latex conversion work.

Sorry, I do not understand your objection about me using "internal" error. I used this to mean the error is coming from internal to Maple. What difference does it make if it comes from int() or Int() or any other function? These all are just internal functions for a user.

@Rouben Rostamian  

"What is the domain of the unknown function u(r,z,t)?   It is
0 < r < 1,    0 < z < 1,   t > 0."

Yes. Actually the book uses 0<r<a and 0<z<H.  To help Maple, I used a numerical value for both radius and height of the cylinder. Both give the same Bessel error.

"Thus, the correct boundary conditions are:
bc := u(r,0,t)=0, u(r,1,t)=0, u(1,z,t)=0, D[1](u)(0,z,t)=0;"

I did not make up the BC. This is a problem from well known book. It is at page 335 from Richard Habeman Applied Partial differential equations, 4th edition.

I see nothing wrong with these boundary conditions. They say that the temperature is zero on the all the cylinder surfaces?

I solved this PDE by hand correctly using these boundary conditions. Because the problem says that it is independent of theta, so I used "u(r,z,t)" only and not "u(r,theta,z,t)" as normally would be the case for a cylinder.  Notice the book uses "a" and "H" which I replaced by "1" each to make it is easier for Maple. Not that changes anything.

"A boundary condition on the edge r=0 is missing"

It is not missing. It is always assumed that at r=0 (center of cylinder) the solution is finite. This is an implicit assumption in such problems. This leads to discarding one of the two Bessel eigenfunctions solutions (the second kind) that results from the radial Bessel ODE since that  blows at r=0 leaving one solution  (The Bessel J_0() of order zero)

@tomleslie 

"where f(r,z) is undefined is going to cause all sorts of weird problems - because it isn't really a meaningful inital condition: it provides no useful information"

Sorry, I am not following what you are saying.

It is not meaningful at all. Initial conditions are needed to obtain the final unknowns in the series solution. This is how half the problems in the books are specified. Need an initial condition!

But thanks for the trick to get the Lapacian without theta dependency. I will use that next time instead of manually type it.

@tomleslie 

The specification of PDE itself is valid inside the domain and when time is involved, for t>0.

The solution of the PDE should ofcourse satisfy the PDE itself and also the initial and boundary conditions.  But this is how the PDE itself is always given/specified in text books:

----------------------------------------

--------------------------------

------------------------

And many more. 

I was just trying to help Maple by giving assumptions (sometimes they help and sometimes they are not needed). In this case, If I use no assumptions, or use t>0 or t>=0, it does not matter, the error in Bessel is there.

I need to check now if this error is new or not. A PDE with Laplacian in Cylinderical or Polar coordinates will generate Bessel ODE for the radial part of the solution, so need to check if this affect other PDE's as well. 

@Thomas Richard 

The Lapacian is correct.  https://en.wikipedia.org/wiki/Del_in_cylindrical_and_spherical_coordinates 

For no theta dependency, the term in red above is dropped. Which gives what I wrote because when expanding the first term above gives

expand(1/r*diff(r*diff(f(r,z,t),r),r))

So you see it is the same as what I had (wikipedica uses "f" instead of "u").

I do not like to do as you did, as it can lead to errors. First of all, time "t" is never part of Laplacian, as Laplacian is spatial only, and I could make errors making substitutions afterwords to remove it, so I type the Lapacian by hand, when there is no theta dependency. Same for Laplacian in spherical. I could not figure how to use Maple for this. But I memorized these things by now so no problem.

@Carl Love 

It does not matter if  r<1 or r<a or no assumption at all on "r".

Initially I had the radius as variable "a", then changed it to "1" just to make it simpler, that is all.

The error shows up in all cases.

@Thomas Richard 

Thanks!  Using "u(x,y__0)" to the set the IC is a nice trick which I will try it next when Maple is having hard time with IC.

your hand simplifcation can not be correct. it is easy to see why. You say it has "beta*mu^3" term at the end but if you type "collect(expr,mu)" where "expr" is your expression, Maple gives

And it is clear that the term multiplying "mu^3" is not "beta" and it does not simplify to "beta" either as your hand simplification says. 

@Mariusz Iwaniuk 

"dsolve convert second order ODE to first order ODE(to Abel)."

I am not sure how dsolve will do that, since "y" is not missing in the original ODE. One can convert second order to first order if "y" is missing.

So letting v(t)=y'(t), then what will dsolve replace y(t) by in the ode? integal of v(t)? Then it becomes integral differential equation, not Abel.

Do you have an example of how dsolve does that to look at to help learn how this is done?

@Mariusz Iwaniuk 

"This is Abel diff equation"

But this is second order ODE. Abel differential equation is first order ODE?