$\renewcommand{\Re}{\operatorname{Re}}$ $\renewcommand{\Im}{\operatorname{Im}}$ $\newcommand{\erf}{\operatorname{erf}}$ $\newcommand{\dag}{\dagger}$ $\newcommand{\const}{\mathrm{const}}$ $\newcommand{\arcsinh}{\operatorname{arcsinh}}$

**Problem 1.**

- Find the solutions that depend only on $r$ of the equation \begin{equation*} \Delta u:=u_{xx}+u_{yy}=0. \end{equation*}
- Find the solutions that depend only on $\rho$ of the equation \begin{equation*} \Delta u:=u_{xx}+u_{yy}+u_{zz}=0. \end{equation*}
- (bonus) In $n$-dimensional case prove that if $u=u(r)$ with $r=(x_1^2+x_2^2+\ldots+x_n^2)^{\frac{1}{2}}$ then \begin{equation} \Delta u = u_{rr}+ \frac{n-1}{r}u_r=0. \label{equ-H8.1} \end{equation}
- (bonus) In $n$-dimensional case prove ($n\ne 2$) that $u=u(r)$ satisfies Laplace equation as $x\ne 0$ iff $u=Ar^{2-n}+B$.

**Problem 2.**
Using the proof of mean value theorem (see Subsection 7.2.4) prove that if $\Delta u\ge 0$ in $B(y,r)$ then

- $u(y)$ does not exceed the mean value of $u$ over the sphere $S(y,r)$ bounding this ball: \begin{equation} u(y)\le \frac{1}{\sigma_n r^{n-1}}\int_{S(y,r)} u\,dS. \label{equ-H8.2} \end{equation}
$u(y)$ does not exceed the mean value of $u$ over this ball $B(y,r)$: \begin{equation} u(y)\le \frac{1}{\omega_n r^n}\int_{B(y,r)} u\,dV. \label{equ-H8.3} \end{equation}

Formulate similar statements for functions satisfying $\Delta u\le 0$ (in the next problem we refer to them as (1)' and (2)'). <!--\end{problem}-->

**Problem 3.**
1. Functions having property (1) (or (2) does not matter) of the
previous problem are called *subharmonic.*

- Functions having property (a)' (or (b)' does not matter) are called
*superharmonic.*

**Problem 4.**
1. Using the proof of maximum principle prove the maximum principle for
subharmonic functions and minimum principle for superharmonic
functions.

Show that minimum principle for subharmonic functions and maximum principle for superharmonic functions do not hold (

*Hint*: construct counterexamples with $f=f(r)$).Prove that if $u,v,w$ are respectively harmonic, subharmonic and superharmonic functions in the bounded domain $\Omega$, coinciding on its boundary ($u|_\Sigma=v|_\Sigma=w|_\Sigma$) then in $w\ge u \ge v$ in $\Omega$.

**Problem 5.**
Using Newton shell theorem (see Subsection 7.3.1) prove that if Earth was a homogeneous solid
ball then the gravity pull inside of it would be proportional to the distance to the center.

**Problem 6.**
Find function $u$ harmonic in $\{x^2+y^2+z^2\le 1\}$ and
coinciding with $g=z^3$ as $x^2+y^2+z^2=1$.

*Hint.* According to Section 8.1 solution
must be a harmonic polynomial of degree $3$ and it should depend only
on $x^2+y^2+z^2$ and $z$ (Explain why). The only way to achive it
(and still coincide with $g$ on $\{x^2+y^2+z^2=1\}$) is to find
\begin{equation*}
u= z^3 + az(1-x^2-y^2-z^2)
\end{equation*}
with unknown coefficient $A$.

**Problem 7.**
Apply method of descent described in Subsection 9.1.4 but to Laplace equation in $\mathbb{R}^2$ and starting from Coulomb potential in $3D$
\begin{equation}
U_3(x,y,z)=-\frac{1}{4\pi} \bigl(x^2+y^2+z^2\bigr)^{-\frac{1}{2}},
\label{equ-H9.1}
\end{equation}
derive logarithmic potential in $2D$
\begin{equation}
U_2(x,y,z)=\frac{1}{2\pi}\log \bigl(x^2+y^2\bigr)^{\frac{1}{2}},
\label{equ-H9.2}
\end{equation}
*Hint.* You will need to calculate diverging integral
$\int_0^\infty U_3 (x,y,z)$. Instead consider
$\int_0^N U_3 (x,y,z)$, subtract constant (f.e. $\int_0^N U_3 (1,0,z)$) and
then tend $N\to \infty$.

**Problem 8.**
Using method of reflection (studied earlier for different equations)
construct Green function for

- Dirichlet problem
- Neumann problem

for Laplace equation in

- half-plane
- half-space

since we know that in the whole plane and space they are just potentials \begin{gather} \frac{1}{2\pi}\log \bigl((x_1-y_1)^2+(x_2-y_2)^2\bigr)^{\frac{1}{2}},\label{equ-H9.3}\\ -\frac{1}{4\pi} \bigl((x_1-y_1)^2+(x_2-y_2)^2+(x_3-y_3)^2\bigr)^{-\frac{1}{2}} \label{equ-H9.4} \end{gather} respectively.

**Problem 9.**
Apply method of descent but now looking for stationary solution of
$-\Delta u=f(x_1,x_2,x_3)$ instead of non-stationary solution of
\begin{align*}
& u_{tt}-\Delta u=f(x_1,x_2,x_3),\\
& u|_{t=0}=g(x_1,x_2,x_3),\\
& u_t|_{t=0}=h(x_1,x_2,x_3)
\end{align*}
start from Kirchhoff formula
(9.1.12) and derive for $n=3$
(7.2.10) with $G(x,y)$ equal to (\ref{equ-H9.4})
here.