Here is what part C looks like in mathematica:
First we tell mathematica that $u(r,\theta)=3r^3\cos^2(\theta)\sin(\theta)-r^3\sin^3(\theta))$
u[r_, t_] = 3 r^3 Cos[t]^2 Sin[t] - r^3 Sin[t]^3;
Then we take the derivatives $u_r$,$u_{rr}$, and $u_{\theta\theta}$
ur = D[u[r, t], r]
urr = D[u[r, t], r, r]
utt = D[u[r, t], t, t]
the output:
9 r^2 Cos[t]^2 Sin[t] - 3 r^2 Sin[t]^3
18 r Cos[t]^2 Sin[t] - 6 r Sin[t]^3
-27 r^3 Cos[t]^2 Sin[t] + 9 r^3 Sin[t]^3
So we have
$$\begin{aligned}
u_r(r,\theta)=r^2(9\cos^2(\theta)\sin(\theta)-&3\sin^3(\theta))
\\u_{rr}(r,\theta)=r(18\cos^2(\theta)\sin(\theta)-&6\sin^3(\theta))
\\u_{\theta\theta}=r^3(-27\cos^2(\theta)\sin(\theta)+&9\sin(\theta))
\end{aligned}$$
Its actually pretty easy to do this by hand at this point since the trig terms are all the same, but Mathematica isn't as error prone as we are.
urr + ur/r + utt/r^2 // FullSimplify
Which gives
0
What a surpirise, the Laplacian is satisfied! The worst part in this whole thing is by far the angular derivative, but its pretty easy to see how all this should go together to zero once we have the derivatives.