Partial Differential Equations/The heat equation

This chapter is about the heat equation, which looks like this:

$\forall (t,x)\in \mathbb {R} \times \mathbb {R} ^{d}:\partial _{t}u(t,x)-\Delta _{x}u(t,x)=f(t,x)$

for some $f:\mathbb {R} \times \mathbb {R} ^{d}\to \mathbb {R}$. Using distribution theory, we will prove an explicit solution formula (if $f$ is often enough differentiable), and we even prove a solution formula for the initial value problem.

Green's kernel and solution

Lemma 6.1:

$\int _{\mathbb {R} }e^{-x^{2}}dx={\sqrt {\pi }}$

Proof:

${\begin{aligned}\left(\int _{\mathbb {R} }e^{-x^{2}}\right)^{2}&=\left(\int _{\mathbb {R} }e^{-x^{2}}\right)\cdot \left(\int _{\mathbb {R} }e^{-y^{2}}\right)&\\&=\int _{\mathbb {R} }\int _{\mathbb {R} }e^{-(x^{2}+y^{2})}dxdy&\\&=\int _{\mathbb {R} ^{2}}e^{-\|(x,y)\|^{2}}d(x,y)&{\text{Fubini}}\\&=\int _{0}^{\infty }\int _{0}^{2\pi }re^{-r^{2}}d\varphi dr&{\text{ integration by substitution using spherical coordinates}}\\&=2\pi \int _{0}^{\infty }re^{-r^{2}}dr&\\&=2\pi \int _{0}^{\infty }{\frac {1}{2{\sqrt {r}}}}{\sqrt {r}}e^{-r}dr&{\text{integration by substitution using }}r\mapsto {\sqrt {r}}\\&=\pi &\end{aligned}}$

Taking the square root on both sides finishes the proof.$\Box$

Lemma 6.2:

$\int _{\mathbb {R} ^{d}}e^{-\|x\|^{2}/2}dx={\sqrt {2\pi }}^{d}$

Proof:

${\begin{aligned}\int _{\mathbb {R} ^{d}}e^{-{\frac {\|x\|^{2}}{2}}}dx&=\overbrace {\int _{-\infty }^{\infty }\cdots \int _{-\infty }^{\infty }} ^{d{\text{ times}}}e^{-{\frac {x_{1}^{2}+\cdots +x_{d}^{2}}{2}}}dx_{1}\cdots dx_{d}&{\text{Fubini's theorem}}\\&=\int _{-\infty }^{\infty }e^{-{\frac {x_{d}^{2}}{2}}}\cdots \int _{-\infty }^{\infty }e^{-{\frac {x_{1}^{2}}{2}}}dx_{1}\,\cdots dx_{d}&{\text{pulling the constants out of the integrals}}\end{aligned}}$

By lemma 6.1,

$\int _{\mathbb {R} }e^{-x^{2}}dx={\sqrt {\pi }}$.

If we apply to this integration by substitution (theorem 5.5) with the diffeomorphism $x\mapsto {\frac {x}{\sqrt {2}}}$, we obtain

${\sqrt {\pi }}=\int _{\mathbb {R} }{\frac {1}{\sqrt {2}}}e^{-{\frac {x^{2}}{2}}}dx$

and multiplying with ${\sqrt {2}}$

${\sqrt {2\pi }}=\int _{\mathbb {R} }e^{-{\frac {x^{2}}{2}}}dx$

Therefore, calculating the innermost integrals first and then pulling out the resulting constants,

$\overbrace {\int _{-\infty }^{\infty }e^{-{\frac {x_{d}^{2}}{2}}}\cdots \int _{-\infty }^{\infty }e^{-{\frac {x_{1}^{2}}{2}}}} ^{d{\text{ times}}}dx_{1}\cdots dx_{d}={\sqrt {2\pi }}^{d}$$\Box$

Theorem 6.3:

The function

$E:\mathbb {R} \times \mathbb {R} ^{d}\to \mathbb {R} ,E(t,x)={\begin{cases}{\sqrt {4\pi t}}^{-d}e^{-{\frac {\|x\|^{2}}{4t}}}&t>0\\0&t\leq 0\end{cases}}$

is a Green's kernel for the heat equation.

Proof:

1.

We show that $E$ is locally integrable.

Let $K\subset \mathbb {R} \times \mathbb {R} ^{d}$ a compact set, and let $T>0$ such that $K\subset (-T,T)\times \mathbb {R} ^{d}$. We first show that the integral

$\int _{(-T,T)\times \mathbb {R} ^{d}}E(s,y)d(s,y)$

exists:

${\begin{aligned}\int _{(-T,T)\times \mathbb {R} ^{d}}E(s,y)d(s,y)&=\int _{(0,T)\times \mathbb {R} ^{d}}E(s,y)d(s,y)&\forall s\leq 0:E(s,y)=0\\&=\int _{0}^{T}\int _{\mathbb {R} ^{d}}{\frac {1}{{\sqrt {4\pi s}}^{d}}}e^{-{\frac {\|y\|^{2}}{4s}}}dyds&{\text{Fubini's theorem}}\end{aligned}}$

By transformation of variables in the inner integral using the diffeomorphism $y\mapsto {\sqrt {2s}}y$, and lemma 6.2, we obtain:

$=\int _{0}^{T}\int _{\mathbb {R} ^{d}}{\frac {{\sqrt {2s}}^{d}}{{\sqrt {4\pi s}}^{d}}}e^{-{\frac {\|y\|^{2}}{2}}}dyds=\int _{0}^{T}1ds=T$

Therefore the integral

$\int _{(-T,T)\times \mathbb {R} ^{d}}E(s,y)d(s,y)$

exists. But since

$\forall (s,y)\in \mathbb {R} \times \mathbb {R} ^{d}:|\chi _{K}(s,y)E(s,y)|\leq |E(s,y)|$

, where $\chi _{K}$ is the characteristic function of $K$, the integral

$\int _{(-T,T)\times \mathbb {R} ^{d}}\chi _{K}(s,y)E(s,y)d(s,y)=\int _{K}E(s,y)d(s,y)$

exists. Since $K$ was an arbitrary compact set, we thus have local integrability.

2.

We calculate $\partial _{t}E$ and $\Delta _{x}E$ (see exercise 1).

$\partial _{t}E(t,x)=\left({\frac {\|x\|^{2}}{4t^{2}}}-{\frac {d}{4t}}\right)E(t,x)$

$\Delta _{x}E(t,x)=\left({\frac {\|x\|^{2}}{4t^{2}}}-{\frac {d}{4t}}\right)E(t,x)$

3.

We show that

$\forall \varphi \in {\mathcal {D}}(\mathbb {R} \times \mathbb {R} ^{d}),(t,x)\in \mathbb {R} \times \mathbb {R} ^{d}:(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )=\delta _{(t,x)}(\varphi )$

Let $\varphi \in {\mathcal {D}}(\mathbb {R} \times \mathbb {R} ^{d}),(t,x)\in \mathbb {R} \times \mathbb {R} ^{d}$ be arbitrary.

In this last step of the proof, we will only manipulate the term $(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )$.

${\begin{aligned}(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )&=T_{E(\cdot -(t,x))}((-\partial _{t}-\Delta _{x})\varphi )&{\text{by definition of distribution derivation}}\\&=\int _{\mathbb {R} \times \mathbb {R} ^{d}}(-\partial _{t}-\Delta _{x})\varphi (s,y)E(s-t,y-x)d(s,y)&\\&=\int _{(t,\infty )\times \mathbb {R} ^{d}}(-\partial _{t}-\Delta _{x})\varphi (s,y)E(s-t,y-x)d(s,y)&\forall t\leq 0:E(t,x)=0\\\end{aligned}}$

If we choose $R>0$ and $T>0$ such that

${\text{supp }}\varphi \subset (-\infty ,t+T)\times B_{R}(x)$

, we have even

$(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )=\int _{(t,t+T)\times B_{R}(x)}(-\partial _{t}-\Delta _{x})\varphi (s,y)E(s-t,y-x)d(s,y)$

Using the dominated convergence theorem (theorem 5.1), we can rewrite the term again:

${\begin{aligned}(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )&=\int _{(t,t+T)\times B_{R}(x)}(-\partial _{t}-\Delta _{x})\varphi (s,y)E(s-t,y-x)d(s,y)\\&=\lim _{\epsilon \downarrow 0}\int _{(t,t+T)\times B_{R}(x)}(-\partial _{t}-\Delta _{x})\varphi (s,y)E(s-t,y-x)(1-\chi _{[t,t+\epsilon ]}(s))d(s,y)\\&=\lim _{\epsilon \downarrow 0}\int _{(t+\epsilon ,t+T)\times B_{R}(x)}(-\partial _{t}-\Delta _{x})\varphi (s,y)E(s-t,y-x)d(s,y)\end{aligned}}$

, where $\chi _{[t,t+\epsilon ]}$ is the characteristic function of $[t,t+\epsilon ]$.

We split the limit term in half to manipulate each summand separately:

${\begin{aligned}\int _{(t+\epsilon ,t+T)\times B_{R}(x)}(-\partial _{t}-\Delta _{x})\varphi (s,y)E(s-t,y-x)d(s,y)\\=-\int _{(t+\epsilon ,t+T)\times B_{R}(x)}\Delta _{x}\varphi (s,y)E(s-t,y-x)d(s,y)\\-\int _{(t+\epsilon ,t+T)\times B_{R}(x)}\partial _{t}\varphi (s,y)E(s-t,y-x)d(s,y)\\\end{aligned}}$

The last integrals are taken over $(t+\epsilon ,t+T)\times B_{R}(x)$ for $\epsilon >0$. In this area and its boundary, $E(s-t,y-x)$ is differentiable. Therefore, we are allowed to integrate by parts.

${\begin{aligned}\int _{(t+\epsilon ,t+T)\times B_{R}(x)}\Delta _{x}\varphi (s,y)E(s-t,y-x)d(s,y)&=\int _{t+\epsilon }^{t+T}\int _{B_{R}(x)}\Delta _{x}\varphi (s,y)E(s-t,y-x)dyds&{\text{Fubini}}\\=\int _{t+\epsilon }^{t+T}\int _{\partial B_{R}(x)}E(s,y)n(y)\cdot \underbrace {\nabla _{x}\varphi (s,y)} _{=0}dyds&-\int _{t+\epsilon }^{t+T}\int _{B_{R}(x)}\nabla _{x}\varphi (s,y)\cdot \nabla _{x}E(s-t,y-x)dyds&{\text{integration by parts in }}y\\=\int _{t+\epsilon }^{t+T}\int _{B_{R}(x)}\varphi (s,y)\Delta _{x}E(s-t,y-x)dyds&-\int _{t+\epsilon }^{t+T}\int _{\partial B_{R}(x)}\underbrace {\varphi (s,y)} _{=0}n(y)\cdot \nabla _{x}E(s-t,y-x)dyds&{\text{integration by parts in }}y\end{aligned}}$

In the last two manipulations, we used integration by parts where $\varphi$ and $f$ exchanged the role of the function in theorem 5.4, and $\nabla _{x}f$ and $\nabla _{x}\varphi$ exchanged the role of the vector field. In the latter manipulation, we did not apply theorem 5.4 directly, but instead with subtracted boundary term on both sides.

Let's also integrate the other integral by parts.

${\begin{aligned}\int _{(t+\epsilon ,t+T)\times B_{R}(x)}\partial _{t}\varphi (s,y)E(s-t,y-x)d(s,y)&=\int _{B_{R}(x)}\int _{t+\epsilon }^{t+T}\partial _{t}\varphi (s,y)E(s-t,y-x)dsdy&{\text{Fubini}}\\=\int _{B_{R}(x)}\underbrace {\varphi (s,y)E(s-t,y-x){\big |}_{s=t+\epsilon }^{s=t+T}} _{=-\varphi (t+\epsilon ,y)E(\epsilon ,y-x)}dy&-\int _{B_{R}(x)}\int _{t+\epsilon }^{t+T}\varphi (s,y)\partial _{t}E(s-t,y-x)dsdy&{\text{integration by parts in }}s\end{aligned}}$

Now we add the two terms back together and see that

${\begin{aligned}(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )&=\lim _{\epsilon \downarrow 0}-\int _{B_{R}(x)}-\varphi (t+\epsilon ,y)E(\epsilon ,y-x)dy\\+\int _{B_{R}(x)}\int _{t+\epsilon }^{t+T}\varphi (s,y)\partial _{t}E(s-t,y-x)dsdy&-\int _{t+\epsilon }^{t+T}\int _{B_{R}(x)}\varphi (s,y)\Delta _{x}E(s-t,y-x)dyds\end{aligned}}$

The derivative calculations from above show that $\partial _{t}E=\Delta _{x}E$, which is why the last two integrals cancel and therefore

$(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )=\lim _{\epsilon \downarrow 0}\int _{B_{R}(x)}\varphi (t+\epsilon ,y)E(\epsilon ,y-x)dy$

Using that ${\text{supp }}\varphi (t+\epsilon ,\cdot )\subset B_{R}(x)$ and with multi-dimensional integration by substitution with the diffeomorphism $y\mapsto x+{\sqrt {2\epsilon }}y$ we obtain:

$\int _{B_{R}(x)}\varphi (t+\epsilon ,y)E(\epsilon ,y-x)dy=\int _{\mathbb {R} ^{d}}\varphi (t+\epsilon ,y)E(\epsilon ,y-x)dy$ $=\int _{\mathbb {R} ^{d}}\varphi (t+\epsilon ,y){\frac {1}{{\sqrt {4\pi \epsilon }}^{d}}}e^{-{\frac {\|y-x\|^{2}}{4\epsilon }}}dy$ $=\int _{\mathbb {R} ^{d}}\varphi (t+\epsilon ,x+{\sqrt {2\epsilon }}y){\frac {{\sqrt {2\epsilon }}^{d}}{{\sqrt {4\pi \epsilon }}^{d}}}e^{-{\frac {\|y\|^{2}}{2}}}dy={\frac {1}{{\sqrt {2\pi }}^{d}}}\int _{\mathbb {R} ^{d}}\varphi (t+\epsilon ,x+{\sqrt {2\epsilon }}y)e^{-{\frac {\|y\|^{2}}{2}}}dy$

Since $\varphi$ is continuous (even smooth), we have

$\forall x\in \mathbb {R} ^{d}:\lim _{\epsilon \to 0}\varphi (t+\epsilon ,x+{\sqrt {2\epsilon }}y)=\varphi (t,x)$

Therefore

${\begin{aligned}(\partial _{t}-\Delta _{x})T_{E(\cdot -(t,x))}(\varphi )&=\lim _{\epsilon \downarrow 0}{\frac {1}{{\sqrt {2\pi }}^{d}}}\int _{\mathbb {R} ^{d}}\varphi (t+\epsilon ,x+{\sqrt {2\epsilon }}y)e^{-{\frac {\|y\|^{2}}{2}}}dy&\\&={\frac {1}{{\sqrt {2\pi }}^{d}}}\int _{\mathbb {R} ^{d}}\varphi (t,x)e^{-{\frac {\|y\|^{2}}{2}}}dy&{\text{dominated convergence}}\\&=\varphi (t,x)&{\text{lemma 6.2}}\\&=\delta _{(t,x)}(\varphi )&\end{aligned}}$$\Box$

Theorem 6.4: If $f:\mathbb {R} \times \mathbb {R} ^{d}$ is bounded, once continuously differentiable in the $t$-variable and twice continuously differentiable in the $x$-variable, then

$u(t,x):=(E*f)(t,x)$

solves the heat equation

$\forall (t,x)\in \mathbb {R} \times \mathbb {R} ^{d}:\partial _{t}u(t,x)-\Delta _{x}u(t,x)=f(t,x)$

Proof:

1.

We show that $(E*f)(t,x)$ is sufficiently often differentiable such that the equations are satisfied.

2.

We invoke theorem 5.?, which states exactly that a convolution with a Green's kernel is a solution, provided that the convolution is sufficiently often differentiable (which we showed in part 1 of the proof).$\Box$

Initial Value Problem

Definition 6.5: Let $\mu :\mathbb {R} \times \mathbb {R} ^{d}\to \mathbb {R}$ and $f:\mathbb {R} ^{d}\to \mathbb {R}$ be two functions. The spatial convolution of $\mu$ and $f$ is given by:

$(\mu *_{x}f)(t,x):=\int _{\mathbb {R} ^{d}}\mu (t,y-x)f(y)dy$

Theorem and definition 6.6: Let $f:[0,\infty )\times \mathbb {R} ^{d}\to \mathbb {R}$ be bounded, once continuously differentiable in the $t$-variable and twice continuously differentiable in the $x$-variable, and let $E:\mathbb {R} ^{d}\to \mathbb {R}$ be continuous and bounded. If we define

${\tilde {f}}:\mathbb {R} \times \mathbb {R} ^{d}\to \mathbb {R} ,{\tilde {f}}(t,x)={\begin{cases}f(t,x)&t\geq 0\\0&t<0\end{cases}}$

, then the function

$u:[0,\infty )\times \mathbb {R} ^{d}\to \mathbb {R} ,u(t,x)={\begin{cases}(E*_{x}g)(t,x)+({\tilde {f}}*E)(t,x)&t>0\\g(x)&t=0\end{cases}}$

is a continuous solution of the initial value problem for the heat equation, that is

${\begin{cases}\forall (t,x)\in (0,\infty )\times \mathbb {R} ^{d}:&\partial _{t}u(t,x)-\Delta _{x}u(t,x)=f(t,x)\\\forall x\in \mathbb {R} ^{d}:&u(0,x)=g(x)\end{cases}}$

Note that if we do not require the solution to be continuous, we may just take any solution and just set it to $g$ at $t=0$.

Proof:

1.

We show

$\forall (t,x)\in (0,\infty )\times \mathbb {R} ^{d}:\partial _{t}u(t,x)-\Delta _{x}u(t,x)=f(t,x)~~~~~(*)$

From theorem 7.4, we already know that ${\tilde {f}}*E$ solves

$\forall (t,x)\in (0,\infty )\times \mathbb {R} ^{d}:\partial _{t}({\tilde {f}}*E)(t,x)-\Delta _{x}({\tilde {f}}*E)(t,x)={\tilde {f}}(t,x){\overset {t>0}{=}}f(t,x)$

Therefore, we have for $\forall (t,x)\in (0,\infty )\times \mathbb {R} ^{d}$,

${\begin{aligned}\partial _{t}u(t,x)-\Delta _{x}u(t,x)=&\partial _{t}(E*_{x}g)(t,x)+\partial _{t}({\tilde {f}}*E)(t,x)\\&-\Delta _{x}(E*_{x}g)(t,x)-\Delta _{x}({\tilde {f}}*E)(t,x)\\=&f(t,x)+\partial _{t}(E*_{x}g)(t,x)-\Delta _{x}(E*_{x}g)(t,x)\end{aligned}}$

which is why $(*)$ would follow if

$\forall (t,x)\in (0,\infty )\times \mathbb {R} ^{d}:\partial _{t}(E*_{x}g)(t,x)-\Delta _{x}(E*_{x}g)(t,x)=0$

This we shall now check.

By definition of the spatial convolution, we have

$\partial _{t}(E*_{x}g)(t,x)=\partial _{t}\int _{\mathbb {R} ^{d}}E(t,x-y)g(y)dy$

and

$\Delta _{x}(E*_{x}g)(t,x)=\Delta _{x}\int _{\mathbb {R} ^{d}}E(t,x-y)g(y)dy$

By applying Leibniz' integral rule (see exercise 2) we find that

${\begin{aligned}\partial _{t}(E*_{x}g)(t,x)-\Delta _{x}(E*_{x}g)(t,x)&=\partial _{t}\int _{\mathbb {R} ^{d}}E(t,x-y)g(y)dy-\Delta _{x}\int _{\mathbb {R} ^{d}}E(t,x-y)g(y)dy&\\&=\int _{\mathbb {R} ^{d}}\partial _{t}E(t,x-y)g(y)dy-\int _{\mathbb {R} ^{d}}\Delta _{x}E(t,x-y)g(y)dy&{\text{ Leibniz' integral rule}}\\&=\int _{\mathbb {R} ^{d}}\left(\partial _{t}E(t,x-y)-\Delta _{x}E(t,x-y)\right)g(y)dy&{\text{ linearity of the integral}}\\&=0&{\text{ exercise 1}}\end{aligned}}$

for all $(t,x)\in (0,\infty )\times \mathbb {R} ^{d}$.

2.

We show that $u$ is continuous.

It is clear that $u$ is continuous on $(0,\infty )\times \mathbb {R} ^{d}$, since all the first-order partial derivatives exist and are continuous (see exercise 2). It remains to be shown that $u$ is continuous on $\{0\}\times \mathbb {R} ^{d}$.

To do so, we first note that for all $(t,x)\in (0,\infty )\times \mathbb {R} ^{d}$

${\begin{aligned}\int _{\mathbb {R} ^{d}}E(t,x-y)dy&=\int _{\mathbb {R} ^{d}}E(t,y)dy&{\text{ integration by substitution using }}y\mapsto x-y\\&=\int _{\mathbb {R} ^{d}}{\sqrt {4\pi t}}^{-d}e^{\frac {\|y\|^{2}}{4t}}dy&\\&=\int _{\mathbb {R} ^{d}}{\sqrt {2\pi }}^{-d}e^{\frac {\|y\|^{2}}{2}}dy&{\text{ integration by substitution using }}y\mapsto {\sqrt {2t}}y\\&=1&{\text{ lemma 6.2}}\end{aligned}}$

Furthermore, due to the continuity of $g$, we may choose for arbitrary $\epsilon >0$ and any $x\in \mathbb {R} ^{d}$ a $\delta >0$ such that

$\forall y\in B_{\delta }(x):|g(y)-g(x)|<\epsilon$.

From these last two observations, we may conclude:

${\begin{aligned}|g(x)-(E*_{x}g)(t,x)|&=\left|1\cdot g(x)-\int _{\mathbb {R} ^{d}}E(t,x-y)g(x)dy\right|&\\&=\left|\int _{\mathbb {R} ^{d}}E(t,x-y)g(x)dy-\int _{\mathbb {R} ^{d}}E(t,x-y)g(x)dy\right|\\&=\left|\int _{B_{\delta }(x)}E(t,x-y)(g(y)-g(x))dy+\int _{\mathbb {R} ^{d}\setminus B_{\delta }(x)}E(t,x-y)(g(y)-g(x))dy\right|&\\&\leq \left|\int _{B_{\delta }(x)}E(t,x-y)(g(y)-g(x))dy\right|+\left|\int _{\mathbb {R} ^{d}\setminus B_{\delta }(x)}E(t,x-y)(g(y)-g(x))dy\right|&{\text{triangle ineq. in }}\mathbb {R} \\&\leq \int _{B_{\delta }(x)}|E(t,x-y)|\underbrace {|g(y)-g(x)|} _{<\epsilon }dy+\int _{\mathbb {R} ^{d}\setminus B_{\delta }(x)}|E(t,x-y)(g(y)-g(x))|dy&{\text{ triangle ineq. for }}\int \\&<\int _{\mathbb {R} ^{d}}|E(t,x-y)|\epsilon dy+\int _{\mathbb {R} ^{d}\setminus B_{\delta }(x)}|E(t,x-y)|\underbrace {(|g(y)|+|g(x)|)} _{\leq 2\|g\|_{\infty }}dy&{\text{ monotony of the }}\int \\&=\epsilon +2\|g\|_{\infty }\left|\int _{\mathbb {R} ^{d}\setminus B_{\delta }(x)}E(t,x-y)dy\right|\end{aligned}}$

But due to integration by substitution using the diffeomorphism $x\mapsto {\sqrt {2t}}x$, we obtain

$\int _{\mathbb {R} ^{d}\setminus B_{\delta }(x)}E(t,x-y)dy=\int _{\mathbb {R} ^{d}\setminus B_{\delta }(0)}E(t,x)dy=\int _{\mathbb {R} ^{d}\setminus B_{\frac {\delta }{\sqrt {2t}}}(0)}{\frac {1}{{\sqrt {2\pi }}^{d}}}e^{-{\frac {\|x\|^{2}}{2}}}dy\to 0,t\to 0$

which is why

$\lim _{t\to 0}|g(x)-(E*_{x}g)(t,x)|<\epsilon$

Since $\epsilon >0$ was arbitrary, continuity is proven.$\Box$