# What is the solution of $u_t=u_{xx}+\frac{1}{x}u_x$?

What does the solution $u(x,t)$ of $u_t=u_{xx}+\frac{1}{x}u_x$ on $[0,1]$ with the following initial condition look like? $u(\frac{1}{2},0)=\delta(\frac{1}{2})$ (i.e. delta function at $x=\frac{1}{2}$)

I tried the separation of variables method as follows but I couldn't derive the precise solution:

Suppose $u(x,t)=X(x)T(t)$, then $\frac{\dot{T}}{T}=c$ and $\frac{X''}{X}+\frac{X'}{xX}=c$. This implies $T(t)=e^{ct}$ and the ODE $xX''+X'-cxX=0$. I guess the solution of the ODE looks like $p(x)e^{\alpha x}$ for some polynomial $p(x)=a_0+a_1x+\ldots$ and some constant $\alpha$. I plugged it into the ODE and I got $(m+1)^2a_{m+1}+(2m+1)\alpha a_m+(\alpha^2-c)a_{m-1}=0$ for any $m=1,2,\ldots.$ I think solutions like $p(x)e^{\alpha x+ct}$ with the above relation for the coefficients of $p$ could form a basis for the space of solutions. But I don't know how to use the initial condition to find the coefficient!

-
The ODE you give can be solved in terms of Bessel functions. But, perhaps more fundamentally, you did not specify boundary conditions for your problem. –  Michael Renardy Feb 5 '13 at 18:32
Thanks for your comment. We can assume Neumann boundary conditions. –  Soli Feb 5 '13 at 19:36
The same post on MSE math.stackexchange.com/questions/295570/… –  Andrew Feb 7 '13 at 15:14