with(DEtools): _Envdiffopdomain:=[D,x]: L:=x^2*D^2+x*D+(x^2-nu^2); gen_exp(L,[D,x],t,x=0); gen_exp(L,[D,x],t,x=infinity); L:=x^2*D^2+x*D-(x^2+nu^2); gen_exp(L,[D,x],t,x=0); gen_exp(L,[D,x],t,x=infinity); L:=x^2*D^2+x*D-(x^2+0^2); formal_sol(L,t,x=0); formal_sol(L,`has logarithm?`,x=0);