Added Oct 6, 2019
Problem 2.4.18 from Peter Olver, Intoduction to Partial differential equations, 4th edition.
Solve for \(u(x,t)\)
With \(u(x,0)=0,u_t(x,0)=g(x)\). Note, in the book, it says to assume \(g(x)\) is even function. In the code below, this assumption is not used. When I find the correct way to implement this assumption in CAS, will have to re-run these.
Mathematica ✗
ClearAll["Global`*"]; pde = D[u[x, t], {t, 2}] - D[u[x, t], {x, 2}] - 2/x*D[u[x,t],x] == 0; ic = {u[x,0]==0,Derivative[0,1][u][x,0]==g[x]}; sol = AbsoluteTiming[TimeConstrained[DSolve[{pde, ic}, u[x, t], {x, t}], 60*10]];
Failed
Maple ✓
restart; pde := diff(u(x,t),t$2)- diff(u(x,t),x$2) - 2/x*diff(u(x,t),x)=0; ic:= u(x,0)=0, D[2](u)(x,0)=g(x); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve([pde,ic],u(x,t)) ),output='realtime'));