<TeXmacs|1.0.6>

<style|<tuple|article|mystyle>>

<\body>
  <doc-data|<doc-title|PDE Summary>>

  <section|General Stuff>

  <\itemize>
    <item><em|Standard mollifier>:

    <\equation*>
      \<eta\>(x)=exp<left|(><frac|1|x<rsup|2>-1><right|)>\<b-1\><rsub|[-1,1]>
    </equation*>

    is a <with|mode|math|C<rsup|*\<infty\>><rsub|c>> hump.

    <\equation*>
      \<eta\><rsub|\<varepsilon\>>(x)=<frac|1|\<varepsilon\><rsup|n>>\<eta\>(x/\<varepsilon\>).
    </equation*>

    Normalization (<with|mode|math|<big|int>=1>) is still missing.

    <item><em|Gamma function>:

    <\equation*>
      \<Gamma\>(\<alpha\>)=<big|int><rsub|0><rsup|\<infty\>>e<rsup|-t>t<rsup|\<alpha\>-1>\<mathd\>t.
    </equation*>

    <item><em|Volumes of sphere and ball>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|\|S<rsup|n-1>\|>|<cell|=>|<cell|\<omega\><rsub|n>r<rsup|n-1>=<frac|2\<pi\><rsup|n/2>|\<Gamma\>(n/2)>r<rsup|n-1>.>>|<row|<cell|\|B<rsup|n>\|>|<cell|=>|<cell|<frac|\<omega\><rsub|n>|n>r<rsup|n>.>>>>
    </eqnarray*>

    <item><em|Green's identities>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|<big|int><rsub|U>v\<Delta\>u>|<cell|=>|<cell|-<big|int><rsub|U>\<nabla\>v\<cdot\>\<nabla\>u+<big|int><rsub|\<partial\>U>v\<partial\><rsub|n>u>>|<row|<cell|<big|int><rsub|U>v\<Delta\>u-u\<Delta\>v>|<cell|=>|<cell|<big|int><rsub|\<partial\>U>v\<partial\><rsub|n>u-u\<partial\><rsub|n>v>>>>
    </eqnarray*>

    <item><em|Young's Inequality>:

    <\equation*>
      <norm|f\<ast\>g|L<rsup|r>|>\<leqslant\><norm|f|L<rsup|p>|><norm|g|L<rsup|q>|><space|1em><with|mode|text|with><space|1em><frac|1|p>+<frac|1|q>=1+<frac|1|r>.
    </equation*>

    In particular <with|mode|math|q=1>, <with|mode|math|r=p>.

    <item><em|Generalized Hölder>:

    <\equation*>
      <norm|f<rsub|1>\<cdot\>f<rsub|2>\<cdots\>f<rsub|m>|L<rsup|1>|>\<leqslant\><norm|f<rsub|1>|p<rsub|1>|><norm|f|p<rsub|2>|>\<cdots\><norm|f<rsub|m>|p<rsub|m>|>
    </equation*>

    if

    <\equation*>
      <frac|1|p<rsub|1>>+<frac|1|p<rsub|2>>+\<cdots\>+<frac|1|p<rsub|m>>=1.
    </equation*>

    <item><em|Interpolation Inequality for <with|mode|math|L<rsup|p>>>: If
    <with|mode|math|1\<leqslant\>s\<leqslant\>r\<leqslant\>t\<leqslant\>\<infty\>>

    <\equation*>
      <frac|1|r>=<frac|\<theta\>|s>+<frac|1-\<theta\>|t>,
    </equation*>

    <with|mode|math|u\<in\>L<rsup|s>\<cap\>L<rsup|t>>, then
    <with|mode|math|u\<in\>L<rsup|r>> and

    <\equation*>
      <norm|u|L<rsup|r>|>\<leqslant\><norm|u|L<rsup|s>|\<theta\>>+<norm|u|L<rsup|t>|1-\<theta\>>.
    </equation*>

    <item><em|Compact>: Every open cover has finite subcover. Metric
    space:<with|mode|math|\<Leftrightarrow\>>sequentially compact.
    Heine-Borel (finite-dim): <with|mode|math|\<Leftrightarrow\>>closed and
    bounded.

    <item><em|Arzelà-Ascoli>: <with|mode|math|(S,d)> compact metric space.
    <with|mode|math|M\<subset\>C(S)> with sup-norm is compact if
    <with|mode|math|M> is bounded, closed and equicontinuous.

    <item><em|Precompact>: has compact closure.

    <item><em|Compact operator>: <with|mode|math|T:B<rsub|1>\<rightarrow\>B<rsub|2>>
    <em|compact> if <with|mode|math|T> continuous and <with|mode|math|T(A)>
    precompact for every bounded <with|mode|math|A>.

    <item><em|Fredholm Alternative>: <with|mode|math|T:B\<rightarrow\>B>
    linear, continuous, compact:

    <\itemize>
      <item>either <with|mode|math|(I-T)x=0> has a nontrivial solution

      <item>or <with|mode|math|(I-T)<rsup|-1>> exists and is bounded.
    </itemize>

    ``Uniqueness and Compactness<with|mode|math|\<Rightarrow\>>Existence''.

    <item><em|Lax-Milgram>: <with|mode|math|B:H\<times\>H\<rightarrow\>\<bbb-F\>>,
    bounded above and coercive<with|mode|math|\<Rightarrow\>><with|mode|math|B[u,g]=F(u)>
    solvable in <with|mode|math|H> for every <with|mode|math|g\<in\>H>.

    Proof: Build operator <with|mode|math|T<rsub|g>:H\<rightarrow\>H<rsup|\<ast\>>>
    that gives <with|mode|math|T<rsub|g>(u)=B[u,g]> (Riesz rep.). Prove 1-1
    and onto.

    The point is: <em|no symmetry>.

    <item><em|Banach-Steinhaus/Uniform Boundedenss Principle>:\ 

    <with|mode|math|X> BR, <with|mode|math|Y> NR,
    <with|mode|math|T<rsub|i>\<in\>L(X,Y)> <with|mode|math|(i\<in\>I)>,
    <with|mode|math|sup<rsub|i\<in\>I>\<\|\|\>T<rsub|i>x\<\|\|\>\<less\>\<infty\>>
    (<with|mode|math|x\<in\>X>)

    <with|mode|math|\<Rightarrow\>><with|mode|math|sup<rsub|i>\<\|\|\>T<rsub|i>\<\|\|\>\<less\>\<infty\>>.

    Read as ``linear+pw bounded<with|mode|math|\<Rightarrow\>>uniformly
    bounded.''
  </itemize>

  <section|Equations>

  <\itemize>
    <item><em|Classification of second order equations>:

    <\equation*>
      A<rsub|i,j>\<partial\><rsub|i>\<partial\><rsub|j>u+B<rsub|i>\<partial\><rsub|i>u+C=0,
    </equation*>

    where <with|mode|math|A> is symmetric WLOG can be rewritten into one of

    <\eqnarray*>
      <tformat|<table|<row|<cell|u<rsub|x x>+u<rsub|y
      y>+<with|mode|text|l.o.d.>>|<cell|=>|<cell|F,>>|<row|<cell|u<rsub|x
      x>-u<rsub|y y>+<with|mode|text|l.o.d.>>|<cell|=>|<cell|F,>>|<row|<cell|u<rsub|x
      x>\<pm\>u<rsub|y>+<with|mode|text|l.o.d.>>|<cell|=>|<cell|F.>>>>
    </eqnarray*>

    <item><em|Minimal surface equation>:

    <\equation*>
      div<left|(><frac|D u|<sqrt|\|D u\|<rsup|2>+1>><right|)>=0
    </equation*>

    <item><em|Monge-Ampére equation>:

    <\equation*>
      det(D<rsup|2>u)=K(x)(1+\|D u\|<rsup|2>)<rsup|(n+2)/2>
    </equation*>
  </itemize>

  <section|Laplace's Equation>

  <with|mode|math|U> open.

  <\itemize>
    <item><with|mode|math|u\<in\>C<rsup|2>(U)>: harmonic, subharmonic
    <with|mode|math|\<Delta\>u\<geqslant\>0>, superharmonic.

    <item><em|Mean Value Inequality:> <with|mode|math|u> subharmonic

    <\eqnarray*>
      <tformat|<table|<row|<cell|u(x)>|<cell|\<leqslant\>>|<cell|<superpose|<big|int>||-><rsub|S(x,r)>u(y)\<mathd\>S<rsub|y>>>|<row|<cell|u(x)>|<cell|\<leqslant\>>|<cell|<superpose|<big|int>||-><rsub|S(x,r)>u(y)\<mathd\>B>>>>
    </eqnarray*>

    (implies <em|Mean Value Property> if harmonic)

    Proof: <with|mode|math|0\<leqslant\><big|int><rsub|B>\<Delta\>u=<big|int><rsub|S>\<partial\><rsub|n>u>,
    then exploit <with|mode|math|\<partial\><rsub|n>u=\<partial\><rsub|r>(x+\<rho\>n)>.
    <with|mode|math|<big|int><rsub|B>u=<big|int><rsub|r><big|int><rsub|\|\<omega\>\|=1>u=u<big|int><rsub|r>.>

    <item><em|Strong Maximum Principle:> <with|mode|math|U> bounded,
    connected, <with|mode|math|u> subharmonic,
    <with|mode|math|u(x)=sup<rsub|U>u><with|mode|math|\<Rightarrow\>><with|mode|math|u>
    constant

    Proof: Consider <with|mode|math|{u=sup}>. By MVI, <with|mode|math|u=sup>
    on any ball in <with|mode|math|U>. Thus <with|mode|math|{u=sup}> open.
    But so is <with|mode|math|{u\<less\>sup}>.
    <with|mode|math|U={u=sup}\<cup\>{u\<less\>sup}>, both
    open<with|mode|math|\<Rightarrow\>><with|mode|math|{u=sup}=U>.

    <item><em|Weak maximum principle:> <with|mode|math|u\<in\>C(<wide|U|\<bar\>>)>
    and subharmonic. Then <with|mode|math|u> assumes extrema on the boundary.

    Proof: SMP or: Suppose <with|mode|math|x\<in\>U> is max and
    <with|mode|math|\<Delta\>u\<gtr\>0>. Then <with|mode|math|D u=0> and
    <with|mode|math|D<rsup|2>u> negative semidef, contradicting
    <with|mode|math|\<Delta\>u=tr(D<rsup|2>u)\<geqslant\>0>. If only
    <with|mode|math|\<Delta\>u\<geqslant\>0>, consider
    <with|mode|math|u+\<varepsilon\>\|x\|<rsup|2>>, which is strictly
    subharmonic.

    <item>Strong<with|mode|math|\<Rightarrow\>>constant,
    Weak<with|mode|math|\<Rightarrow\>>extrema on boundary.

    <item>Uniqueness follows directly from the WMP.

    <item><em|Harnack's Inequality:> <with|mode|math|u\<geqslant\>0> (!)
    harmonic, <with|mode|math|U<rprime|'>\<subset\>\<subset\>U>
    connected<with|mode|math|\<Rightarrow\>><with|mode|math|\<exists\>C> such
    that <with|mode|math|sup u\<less\>C inf u.>

    Proof: Pick <with|mode|math|x<rsub|1>,x<rsub|2>\<in\>U>, apply MVP for
    large and small circle, respectively, then shrink/expand domain by using
    <with|mode|math|u\<geqslant\>0>, take sup/inf. Use cover of balls to
    repeat argument as necessary.

    <item><em|Fundamental solution:> look for radial symmetry

    <\equation*>
      \<psi\>=C+<choice|<tformat|<table|<row|<cell|<frac|1|2\<pi\>>log
      r>|<cell|n=2,>>|<row|<cell|<frac|1|(2-n)\<omega\><rsub|n>>r<rsup|2-n>>|<cell|n\<geqslant\>3.>>>>>
    </equation*>

    Constant chosen because it gives the right constant to prove
    <with|mode|math|\<Delta\>\<psi\>=\<delta\><rsub|0>> (use Green's second
    id on a ball surrounding the signularity).
    <with|mode|math|K(x,\<xi\>)=\<psi\>(\|x-\<xi\>\|)>.

    <item><em|Liouville's Theorem:> (only in 2D) Subharmonic functions
    bounded above are constant.

    <item><with|mode|math|u\<in\>C<rsup|2>(<wide|U|\<bar\>>):>

    <\equation>
      <label|eq:laplace-solrep>u(\<xi\>)=<big|int><rsub|U>K(x,\<xi\>)\<Delta\>u\<mathd\>x+<big|int><rsub|\<partial\>U>u\<partial\><rsub|n<rsub|x>>K(x,\<xi\>)-K(x,\<xi\>)\<partial\><rsub|n<rsub|x>>u\<mathd\>S<rsub|x>.
    </equation>

    Proof: Integrate on <with|mode|math|U\<setminus\>B<rsub|\<varepsilon\>>>,
    <with|mode|math|\<varepsilon\>\<rightarrow\>0>.

    Remains valid if <with|mode|math|K> replaced by <with|mode|math|K+w> with
    harmonic <with|mode|math|w>.

    <item><em|Green's function for Dirichlet problem:>
    <with|mode|math|\<Delta\><rsub|x>G=\<delta\><rsub|\<xi\>>>,
    <with|mode|math|G(x,\<xi\>)=0> for <with|mode|math|x\<in\>\<partial\>U>.
    Use <with|mode|math|G> in (<reference|eq:laplace-solrep>). To get one, we
    need to find <with|mode|math|w> with <with|mode|math|w=-K> on
    <with|mode|math|\<partial\>U>. (Use method of images.) For a ball, we get
    the <em|Poisson kernel>

    <\equation*>
      H(x,\<xi\>)=<frac|r<rsup|2>-\|\<xi\>\|<rsup|2>|\<omega\><rsub|n>r\|x-\<xi\>\|<rsup|n>>
    </equation*>

    <em|Poisson's integral formula:>

    <\equation*>
      u(\<xi\>)=<big|int><rsub|S(0,r)>H(x,\<xi\>)f(x)\<mathd\>S<rsub|\<xi\>>.
    </equation*>

    <item><em|Kelvin's transformation>: <with|mode|math|u>
    harmonic<with|mode|math|\<Rightarrow\>>

    <\equation*>
      \|x\|<rsup|2-n>u(x/\|x\|<rsup|2>) <with|mode|text|harmonic for
      <with|mode|math|x\<neq\>0>>.
    </equation*>

    <item>Properties of <with|mode|math|H>:

    <\itemize>
      <item><with|mode|math|H(x,\<xi\>)=H(\<xi\>,x)>

      <item><with|mode|math|H(x,\<xi\>)\<gtr\>0> on <with|mode|math|B(0,r)>

      <item><with|mode|math|\<Delta\><rsub|\<xi\>>H(x,\<xi\>)=0> for
      <with|mode|math|\<xi\>\<in\>B(0,r)> and <with|mode|math|x\<in\>S(0,r)>

      <item><with|mode|math|<big|int><rsub|S(0,1)>H(x,\<xi\>)\<mathd\>S<rsub|x>=1>
    </itemize>

    <item><em|Existence on a ball:> also gives
    <with|mode|math|C(<wide|B|\<bar\>>)>\ 

    Proof: Differentiate under integral (using DCT). Prove continuity onto
    the boundary by

    <\equation*>
      u(\<xi\>)-f(y)=<big|int><rsub|S(\<xi\>,r)>H(x,\<xi\>)(f(x)-f(y))\<mathd\>S<rsub|x>
    </equation*>

    Use <with|mode|math|\<varepsilon\>>-<with|mode|math|\<delta\>>-continuity
    of <with|mode|math|f> and split integral into
    <with|mode|math|\|x-y\|\<less\>\<delta\>> and
    <with|mode|math|\|x-y\|\<gtr\>\<delta\>>. (Method called <em|approximate
    identities>.)

    <item><em|Converse of MVP>: <with|mode|math|u\<in\>C(U)>
    harmonic<with|mode|math|\<Leftrightarrow\>>satisifes MVP for every
    <with|mode|math|B(x,r)\<subset\>U>.

    Proof: Construct a harmonic function <with|mode|math|v> on
    <with|mode|math|B(x,r)> with <with|mode|math|v=u> on
    <with|mode|math|S(x,r)>. <with|mode|math|v-u> satisfies MVP on any
    subcircle, thus it satisfies the strong maximum principle. Thus
    <with|mode|math|v=u>.

    <item><em|Real analytic:> completely represented by absolutely convergent
    Taylor series.

    <with|mode|math|\<exists\>M\<gtr\>0\<forall\>\<alpha\>:
    \|\<partial\><rsup|\<alpha\>>f(y)\|\<leqslant\><frac|M\|\<alpha\>\|!|r<rsup|\|\<alpha\>\|>>><with|mode|math|\<Leftrightarrow\>><with|mode|math|analytic>.

    Real analytic <with|mode|math|f> is completely determined by power series
    (use similar open-set method on <with|mode|math|{\<partial\><rsup|\<alpha\>>h(y)=0\<forall\>\<alpha\>}>
    as SMP)

    <item><em|Harmonic<with|mode|math|\<Rightarrow\>>Analytic:> Consider
    <with|mode|math|H(x,\<xi\>+i\<sigma\>)>. Find a region of
    <with|mode|math|\<sigma\>> where <with|mode|math|H> is differentiable.

    <item>Analyticity estimates can be obtained by the MVP applied to
    <with|mode|math|\<partial\><rsub|x<rsub|j>>u>, then coordinatewise Gauÿ,
    giving

    <\equation*>
      \|\<partial\><rsub|x<rsub|j>>u(x)\|\<leqslant\><frac|n|r>max<rsub|S(x,r)>\|u\|\<leqslant\><frac|n|r>sup<rsub|U>\|u\|.
    </equation*>

    Then iterate this estimate with <with|mode|math|1/\|\<alpha\>\|> radius
    to get

    <\equation*>
      \|\<partial\><rsup|\<alpha\>>u(x)\|\<leqslant\><left|(><frac|n\|\<alpha\>\||r><right|)><rsup|\|\<alpha\>\|>max<rsub|S(x,r)>\|u\|.
    </equation*>

    <item>Uniformly (on compact subsets of <with|mode|math|U>) converging
    sequences of harmonic functions converge to harmonic functions.

    Proof: Limit is continuous (because of uniform convergence). Now exchange
    limits (DCT) in MVP and prove harmonicity.

    <item><em|Harnack's convergence theorem:> <with|mode|math|u<rsub|k>>
    harmonic, increasing and bounded at a point. Then
    <with|mode|math|(u<rsub|k>)> converges uniformly on compact subsets to a
    harmonic function.

    Proof: above + Harnack inequality.

    <item><em|``Montel's Theorem''>--a compactness criterion:\ 

    <with|mode|math|(u<rsub|k>)> bounded,
    harmonic<with|mode|math|\<Rightarrow\>><with|mode|math|\<exists\>>uniformly
    (on compact subsets) converging subsequence
    <with|mode|math|\<rightarrow\>> harmonic limit.

    Proof: <with|mode|math|(u<rsub|k>)> is equicontinuous because of the
    derivative estimates and the assumed uniform bound.

    <item><em|Subharmonicity on C(U):> Satisfies MVI locally.

    <item>Perron's method:

    <\itemize>
      <item><with|mode|math|S<rsub|f>:={v\<in\>C(<wide|U|\<bar\>>),
      v\<leqslant\>BC, v <with|mode|text|subharmonic>}>.

      <item><with|mode|math|u\<assign\>sup S<rsub|f>> is harmonic.

      Proof:

      <\itemize>
        <item><with|mode|math|S<rsub|f>> is closed under finite max. (MVI)

        <item><em|Harmonic lifting:> <with|mode|math|v> subharmonic,

        <\equation*>
          V(x)=<choice|<tformat|<table|<row|<cell|<with|mode|text|harmonic
          function with matching BCs>>|<cell|B(\<xi\>,r),>>|<row|<cell|v>|<cell|<with|mode|text|elsewhere>.>>>>>
        </equation*>

        <with|mode|math|v\<in\>S<rsub|f>\<Rightarrow\>V\<in\>S<rsub|f>>,
        <with|mode|math|v\<leqslant\>V>.

        <item>Fix a closed ball, grab sequence
        <with|mode|math|v<rsub|k>\<rightarrow\>u> at a point
        <with|mode|math|\<xi\>>. <with|mode|math|<wide|v|\<bar\>><rsub|k>\<assign\>max(v<rsub|1>,\<ldots\>,v<rsub|k>,min
        BC)>.

        <item>Replace these by their harmonic lifting
        <with|mode|math|V<rsub|k>> around <with|mode|math|\<xi\>>.

        <item>HCT for a limit <with|mode|math|V>.

        <item>Prove <with|mode|math|V=u> on ball by finding SMP uniqueness of
        harmonic liftings of in-between (<with|mode|math|V\<less\>u>)
        functions.
      </itemize>
    </itemize>

    <item><em|Barrier function at <with|mode|math|y\<in\>\<partial\>U>/regular
    boundary point:>

    <with|mode|math|w\<in\>C(<wide|U|\<bar\>>)> subharmonic,
    <with|mode|math|w(y)=0>, <with|mode|math|w(\<partial\>U\<setminus\>{y})\<less\>0>.

    <with|mode|math|\<exists\>>tangent plane<with|mode|math|\<Rightarrow\>>regular

    <with|mode|math|\<exists\>>exterior sphere<with|mode|math|\<Rightarrow\>>barrier<with|mode|math|=><with|mode|math|K(<with|mode|text|boundary
    point>, <with|mode|text|outside center>)-K(x,<with|mode|text|outside
    center>)>

    <with|mode|math|><with|mode|math|\<exists\>>exterior
    cone<with|mode|math|\<Rightarrow\>>regular

    <item>At regular boundary points, <with|mode|math|u=BC>.

    Proof:\ 

    <\itemize>
      <item>Fix <with|mode|math|\<varepsilon\>\<gtr\>0>.
      <with|mode|math|\<delta\>> from <with|mode|math|\<varepsilon\>>-<with|mode|math|\<delta\>>
      with <with|mode|math|f>.

      <item><with|mode|math|v=BC+A\<cdot\>barrier-\<varepsilon\>>, where
      <with|mode|math|A*w\<leqslant\>-2 max BC> outside a ball around the
      boundary point in question. <with|mode|math|v> subharmonic by def.

      <item>Show <with|mode|math|v\<leqslant\>f(x)> on boundary and interior.

      <item>Do some funky tricks involving <with|mode|math|-f>, its Perron
      function, and the maximum principle to show opposite inequality.
    </itemize>

    <item>The Dirichlet problem is solvable for all continuous BC data iff
    the domain is regular.
  </itemize>

  <subsection|Energy Methods>

  <\itemize>
    <item><with|mode|math|0=<big|int>w\<Delta\>w=<big|int>\|\<nabla\>w\|<rsup|2>>
    proves uniqueness in <with|mode|math|C<rsup|2>(<wide|U|\<bar\>>)>.

    <item><em|Energy Functional:>

    <\equation*>
      I[w]=<big|int><rsub|U><frac|1|2>\|\<nabla\>w\|<rsup|2>+w*g\<mathd\>x
    </equation*>

    for <with|mode|math|g> the RHS.

    <item><em|Dirichlet's principle:> <with|mode|math|u\<in\>C<rsup|2>(<wide|U|\<bar\>>)>
    solves PDE+BC<with|mode|math|\<Leftrightarrow\>>it minimizes
    <with|mode|math|I[u]> over <with|mode|math|{w\<in\>C<rsup|2>(<wide|U|\<bar\>>),w=RHS<with|mode|text|
    on >\<partial\>\<Omega\>}>.

    Proof: PDE<with|mode|math|\<Rightarrow\>>min: Start from

    <\equation*>
      0=<big|int>(-\<Delta\>u+g)(u-w),
    </equation*>

    use Gauÿ, Cauchy-Schwarz, <with|mode|math|<sqrt|a><sqrt|b>\<leqslant\>1/2(a<rsup|2>+b<rsup|2>)>.

    min<with|mode|math|\<Rightarrow\>>PDE: <with|mode|math|w=u+t*v>, for
    <with|mode|math|v\<in\>C<rsup|\<infty\>><rsub|c>>. Differentiate by
    <with|mode|math|t>.
  </itemize>

  <subsection|Potentials>

  <\itemize>
    <item><em|Potential of a measure>:

    <\equation*>
      u<rsub|\<mu\>>(x)=<frac|2-n|\<omega\><rsub|n>><big|int><rsub|\<bbb-R\><rsup|n>>K(x,y)\<mu\>(\<mathd\>y)=<big|int><rsub|\<bbb-R\><rsup|n>>\|x-y\|<rsup|2-n>\<mu\>(\<mathd\>y)
    </equation*>

    <item>Computable for a sphere with uniform charge density (same as point
    charge), finite line, disk.

    <item><with|mode|math|u<rsub|\<mu\>>=0\<Rightarrow\>\<mu\>=0>.

    Proof: Show <with|mode|math|\<mu\>\<ast\>f=0> for any
    <with|mode|math|f\<in\>C<rsup|\<infty\>><rsub|c>> by

    <\equation*>
      \<mu\>\<ast\>f=\<mu\>\<ast\>(K\<ast\>\<Delta\>f)=(\<mu\>\<ast\>K)\<ast\>\<Delta\>f=0.
    </equation*>

    <item><em|Potentials of compact set>: Harmonic function with BC 1 on
    compact set <with|mode|math|F> and BC zero at infinity. Perron function
    on ever-increasing balls--independent of exact domains.

    <item><em|A (unique) generating (positive) measure on
    <with|mode|math|\<partial\>F> exists:>

    Proof (if <with|mode|math|\<partial\>F\<in\>C<rsup|2>>): by Poisson's
    boundary representation formula (with both <with|mode|math|u> and
    <with|mode|math|\<partial\><rsub|n>u>)

    <\equation*>
      p<rsub|F>(\<xi\>)=<big|int><rsub|\<partial\>F>K(x,\<xi\>)<wide*|\<partial\><rsub|n>p<rsub|F>\<mathd\>S<rsub|x>|\<wide-underbrace\>><rsub|<with|mode|text|measure!>>.
    </equation*>

    <with|mode|math|\<partial\><rsub|n>u\<leqslant\>0> by the max principle
    (1 on the boundary must be the max value)<with|mode|math|\<Rightarrow\>>positivity.

    Proof (if not):

    <\itemize>
      <item>Approximate <with|mode|math|F> through shrinking compact sets
      with <with|mode|math|C<rsup|\<infty\>>> boundary
      (<with|mode|math|1/k<rsup|2>>-mollified indicators of
      <with|mode|math|F<rsup|1/k>={dist(x,F)\<leqslant\>1/k}>.
      <with|mode|math|\<psi\>=\<varphi\><rsub|1/k<rsup|2>>\<ast\>\<b-1\><rsub|F<rsup|1/k>>>.
      Then consider <with|mode|math|F<rsup|1/2k>\<subset\>\<psi\><rsup|-1>([c,1])\<subset\>F<rsup|1/k>>
      and use Sard's Theorem to deduce boundary smoothness for a.e.
      <with|mode|math|c>. Generate <with|mode|math|\<mu\><rsub|k>> by above
      theorem.

      <item><with|mode|math|p<rsub|F<rsub|k>>\<rightarrow\>p<rsub|F>>
      uniformly on compact subsets (Harnack)

      <item>Prove <with|mode|math|\<mu\><rsub|k>(\<bbb-R\><rsup|n>)\<leqslant\>R<rsup|n-2>>
      by using a <with|mode|math|B(0,R)\<supset\>F<rsub|k>>--use Fubini and
      the generator of the disk potential. (``<em|Gauÿ' trick>') Thus
      <with|mode|math|\<exists\>>weak-* convergent subsequence supported on
      <with|mode|math|\<partial\>F>. Thus convergene of
      <with|mode|math|p<rsub|F<rsub|k>>\<rightarrow\>p<rsub|F>> away from
      <with|mode|math|\<partial\>F>. Uniqueness by uniqueness of potentials
      of measures.
    </itemize>
  </itemize>

  <subsection|Lebesgue's Thorn>

  <\itemize>
    <item>In 2D, Riemann mapping theorem guarantees that point regularity is
    topological, not geometric.

    <item>Lebesgue's Thorn: Using level sets of the potential of the measure
    <with|mode|math|x<rsup|\<beta\>>\<mathd\>x> on <with|mode|math|(0,1)>,
    one may construct exceptional points.
  </itemize>

  <subsection|Capacity>

  <\itemize>
    <item>

    <\equation*>
      cap(F)=\<mu\><rsub|F>(\<bbb-R\><rsup|n>)=<frac|2-n|\<omega\><rsub|n>><big|int><rsub|\<partial\>F<with|mode|text|
      or enclosing surface>>\<partial\><rsub|n>p<rsub|F>\<mathd\>S<rsub|x>.
    </equation*>

    <item>If <with|mode|math|\<partial\>F\<in\>C<rsup|2>>, Green's 1st id
    gives

    <\equation*>
      cap(F)=<frac|2-n|\<omega\><rsub|n>><big|int><rsub|U\<subset\>\<bbb-R\><rsup|n>\<setminus\>F>\|\<nabla\>p<rsub|F>\|<rsup|2>.
    </equation*>

    <item><em|Wiener's criterion:> <with|mode|math|y\<in\>\<partial\>U>
    regular <with|mode|math|\<Leftrightarrow\>>

    <\equation*>
      \<lambda\><rsup|2-n><big|sum><rsub|k=0><rsup|\<infty\>>\<lambda\><rsup|k(2-n)>cap(F<rsub|k>)<space|2em>F<rsub|k>\<assign\>{\<lambda\><rsup|k+1>\<leqslant\>\|x-y\|\<leqslant\>\<lambda\><rsup|k>}<space|2em>(\<lambda\>\<in\>(0,1)).
    </equation*>

    <item>Properties of capacity:

    <\itemize>
      <item><with|mode|math|F<rsub|1>\<subset\>F<rsub|2>\<Rightarrow\>cap(F<rsub|1>)\<leqslant\>cap(F<rsub|2>)>
      (<em|Gauÿ' Trick!>)

      \;

      <\equation*>
        cap(F<rsub|1>)=<big|int><rsub|\<bbb-R\><rsup|n>>\<mu\><rsub|1>(\<mathd\>x)=<big|int><rsub|\<bbb-R\><rsup|n>>p<rsub|2>\<mu\><rsub|1>(\<mathd\>x)=<big|int><big|int>\|x-y\|<rsup|2-n>\<mu\><rsub|2>(\<mathd\>y)\<mu\><rsub|1>(\<mathd\>y)=<big|int>p<rsub|1>\<mu\><rsub|2>(\<mathd\>y)\<leqslant\>cap(F<rsub|2>).
      </equation*>

      <item><with|mode|math|(F<rsub|k>)> nested sequence with
      <with|mode|math|<big|cap>F<rsub|k>=F>, then
      <with|mode|math|cap(F<rsub|k>)\<rightarrow\>cap(F)>.

      (smooth <with|mode|math|\<varphi\>=1> on <with|mode|math|F<rsub|1>>,
      <with|mode|math|cap(F)=<big|int>\<varphi\>\<mu\><rsub|F>\<leftarrow\><big|int>\<varphi\>\<mu\><rsub|F<rsub|k>>=cap(F<rsub|k>)>)

      <item><with|mode|math|cap(A\<cup\>B)\<leqslant\>cap(A)+cap(B)>.

      (<with|mode|math|p<rsub|\<cup\>>\<leqslant\>p<rsub|A>+p<rsub|B>> by
      WMP. Then use Gauÿ' trick.)

      <item><with|mode|math|cap(A\<cup\>B)+cap(A\<cap\>B)\<leqslant\>cap(A)+cap(B)>
    </itemize>

    <item><with|mode|math|cap(<wide|B(0,R)|\<bar\>>)=cap(S(0,R))=R<rsup|n-2>>.

    <item><em|Screening>: nested spheres <with|mode|math|A\<subset\>B>.
    <with|mode|math|cap(A\<cup\>B)=cap(B)> (think of the potentials)

    <item><with|mode|math|cap(F)=sup{\<mu\>(F):supp(\<mu\>)\<subset\>F,
    u<rsub|\<mu\>>(F)\<leqslant\>1}> (Smooth approx
    <with|mode|math|F<rsub|k>> to <with|mode|math|F> so that
    <with|mode|math|p<rsub|F<rsub|k>>=1> on <with|mode|math|\<partial\>F>.
    Then Gauÿ' trick.)

    <item><em|Coulomb energy>:

    <\equation*>
      E[\<mu\>]=<frac|1|2><big|int><big|int>\|x-y\|<rsup|2-n>\<mu\>(\<mathd\>x)\<mu\>(\<mathd\>y).
    </equation*>

    Mutual energy:

    <\equation*>
      E[\<mu\>,\<nu\>]=<frac|1|2><big|int><big|int>\|x-y\|<rsup|2-n>\<mu\>(\<mathd\>x)\<nu\>(\<mathd\>y).
    </equation*>

    <item>Properties:

    <\itemize>
      <item>If <with|mode|math|E[\|\<mu\>\|]\<less\>\<infty\>>, then pos.def.

      <item>CSU

      <item><with|mode|math|\<mu\>\<mapsto\>E[\<mu\>]> strictly convex
    </itemize>

    <item><em|Gauÿ' principle:> <with|mode|math|\<mu\>\<geqslant\>0> finite
    measure on <with|mode|math|F>.

    <\equation*>
      G[\<mu\>]=E[\<mu\>]-\<mu\>(F)\<geqslant\>-<frac|1|2>cap(F)
    </equation*>

    Proof:

    <\itemize>
      <item><with|mode|math|G(\<mu\>)> bounded below (<with|mode|math|F>
      compact<with|mode|math|\<Rightarrow\>><with|mode|math|\|x-y\|> bdd.)

      <item>Infimizing sequences are precompact (i.e. have bounded
      <with|mode|math|\<mu\><rsub|k>(F)>)

      <item><with|mode|math|G> is wlsc (take infimizing sequence
      <with|mode|math|(\<mu\><rsub|k>)>, use <with|mode|math|max(M,\|x-y\|)>
      to cut off, <with|mode|math|k\<rightarrow\>\<infty\>>,
      <with|mode|math|M\<rightarrow\>\<infty\>> (MCT), consider
      <with|mode|math|E[\<mu\>-\<mu\><rsub|k>]>)

      <item>Minimizer is unique (strict convexity)

      <item>Minimizer is <with|mode|math|\<mu\><rsub|F>> (Consider
      Euler-Lagrange Equation)

      <item>Evaluate minimum
    </itemize>

    <item><em|Kelvin's principle:>

    <\equation*>
      <frac|1|2cap(F)>=inf{E[\<mu\>]:\<mu\>\<geqslant\>0,
      supp(\<mu\>)\<subset\>F,\<mu\>(F)=1}.
    </equation*>

    Proof: Apply Gauÿ' principle to <with|mode|math|t\<mu\>>, choose
    <with|mode|math|t=cap(F)>.
  </itemize>

  <section|Heat Equation>

  <\itemize>
    <item><em|Conservation of mass>: <with|mode|math|\<partial\><rsub|t>u+div(\<b-v\>)=0>

    <item><em|Fick's law>: <with|mode|math|\<b-v\>=-\<alpha\><rsup|2>\<nabla\>u>.

    <item>Together: <with|mode|math|u<rsub|t>=\<Delta\>u>.

    <item>Parabolic scaling invariance: <with|mode|math|x\<mapsto\>\<lambda\>x>,
    <with|mode|math|t\<mapsto\>\<lambda\><rsup|2>t>.

    <item>Use conservation of mass (<with|mode|math|\<partial\><rsub|t><big|int>u=0>)
    to obtain the ansatz <with|mode|math|u(x,t)=t<rsup|-n/2>g(r*t<rsup|-1/2>)>.
    Plug in heat equation to get the heat kernel

    <\equation*>
      k(x,t)=<frac|1|(4\<pi\>t)<rsup|n/2>>e<rsup|-\|x\|<rsup|2>/4t>.
    </equation*>

    <item>Use

    <\equation*>
      2<big|int><rsub|y\<gtr\>a>e<rsup|-y<rsup|2>>\<mathd\>y\<less\>2<big|int><rsub|y\<gtr\>a><frac|y|a>e<rsup|-y<rsup|2>>=<frac|e<rsup|-a<rsup|2>>|a>.
    </equation*>

    and in-boxing the ball to show

    <\equation*>
      <big|int><rsub|\|x\|\<geqslant\>\<delta\>>k(x,t)\<mathd\>x\<rightarrow\>0<space|1em><with|mode|text|as><space|1em>t\<rightarrow\>0.
    </equation*>

    <item><with|mode|math|u=k\<ast\>f> solves
    <with|mode|math|u<rsub|t>=\<Delta\>u> for
    <with|mode|math|u\<rightarrow\>f> for <with|mode|math|t\<rightarrow\>0>.

    <item><em|Tychonoff counterexample> for uniqueness:

    <\equation*>
      u(x,t)=<big|sum><rsub|k>g<rsub|k>(t)x<rsup|2k>
    </equation*>

    <item><em|Widder's Theorem:> <with|mode|math|u\<geqslant\>0><with|mode|math|\<Rightarrow\>>uniqueness.

    <item><em|Heat ball:> <with|mode|math|E(x,t,r)={k(x-y,t-r)\<geqslant\>r<rsup|-n>}.>

    <item><with|mode|math|V<rsub|T>=U\<times\>[0,T]>,\ 

    <with|mode|math|\<partial\><rsub|1>V<rsub|T>=>all except top ``lid'',\ 

    <with|mode|math|\<partial\><rsub|2>V<rsub|T>>=lid.

    <item><em|Mean Value Property:> <with|mode|math|><with|mode|math|u\<in\>C<rsup|2>(V<rsub|T>)>,
    <with|mode|math|\<partial\><rsub|t>u-\<Delta\>u\<leqslant\>0>,
    <with|mode|math|E(\<ldots\>)\<subset\>V<rsub|T>>:

    <\equation*>
      u(x,t)\<leqslant\><frac|1|4r<rsup|n>><big|int><big|int><rsub|E(x,t,r)>u(y,s)<frac|\|x-y\|<rsup|2>|(t-s)<rsup|2>>\<mathd\>y\<mathd\>s
    </equation*>

    <\itemize>
      <item>Exists for heat spheres as well.

      <item><em|Converse:> Equality and <with|mode|math|C<rsup|2>(V<rsub|T>)>
      implies <with|mode|math|\<partial\><rsub|t>u=\<Delta\>u>.
    </itemize>

    Proof: Let RHS=<with|mode|math|\<varphi\>(r)>.
    <with|mode|math|\<varphi\>(0)=u(x,t)>,

    <\equation*>
      \<varphi\><rprime|'>(r)=-C<big|int>(\<partial\><rsub|s>u-\<Delta\>u)\<psi\>*\<mathd\>y\<mathd\>s\<geqslant\>0
    </equation*>

    with <with|mode|math|{\<psi\>\<geqslant\>0}=E(\<ldots\>)>.

    <item><em|Strong Maximum Principle:> <with|mode|math|U> open, bounded,
    connected, <with|mode|math|u\<in\>C(<wide|V|\<bar\>><rsub|T>)> and
    satisfies MVI. Then

    <\equation*>
      max<rsub|<wide|V|\<bar\>><rsub|T>>u\<leqslant\>max<rsub|\<partial\><rsub|1>V<rsub|T>>u.
    </equation*>

    If max attained at <with|mode|math|(x,t)\<in\>V<rsub|T>>, then
    <with|mode|math|u> is \ constant in <with|mode|math|<wide|V|\<bar\>><rsub|t>>.

    Proof: If max attained in interior, then <with|mode|math|u=M> on heat
    ball. Then a polygonal path reaches every point on
    <with|mode|math|V<rsub|T>>.

    <item><with|color|red|<em|Temperatures are analytic:>>

    <item><with|color|red|<em|Green's functions for the heat equation:>>

    <item><em|Strong Converse of MVP>.

    Proof: Construct parallel solution by Green's functions. Conclude
    uniqueness by MVP.
  </itemize>

  <subsection|Difference Schemes and Probabilistic Interpretation>

  <\itemize>
    <item>Work on a lattice.

    <item><em|Strong Maximum Principle> (subharmonic<with|mode|math|\<Rightarrow\>>assume
    max <with|mode|math|M> in interior<with|mode|math|\<Rightarrow\>><with|mode|math|M=u\<leqslant\>E[x+h\<omega\>]\<leqslant\>M>.)

    <item>Implies discrete Laplacian has trivial
    null-space<with|mode|math|\<Rightarrow\>><with|mode|math|\<exists\>!>

    <item>Allows <em|Discrete Poisson Integral Formula>. (by solving for
    <with|mode|math|\<delta\>> on the boundary)

    <item><em|Markov property>: <with|mode|math|E[X<rsub|m+1>\|X<rsub|1>,\<ldots\>,X<rsub|m>]=E[X<rsub|m+1>\|X<rsub|m>]>.

    <item><em|(Super)Martingale property>: <with|mode|math|u>
    subharmonic<with|mode|math|\<Rightarrow\>><with|mode|math|E[u(X<rsub|m+1>)\|X<rsub|m>]\<geqslant\>u(X<rsub|m>)>
    (just like discrete SMP) [with <with|mode|math|X<rsub|m>> a random walk]

    <item><em|Strong Martingale Property:> <with|mode|math|m> may be a
    stopping time.

    <item>If <with|mode|math|M<rsub|U>> is first passage time to
    <with|mode|math|\<partial\>U>, then <with|mode|math|><with|mode|math|u=E[f(x+W<rsub|M<rsub|U>>)]>.
    (<with|mode|math|f>=BC, <with|mode|math|u> harmonic)

    <\equation*>
      E[f(x+W<rsub|M<rsub|U>>)]=<big|sum><rsub|y\<in\>\<partial\>U<rsub|h>>H(x,y)f(y)=<big|sum><rsub|y\<in\>\<partial\>U<rsub|h>><wide*|P(hit
      y)|\<wide-underbrace\>><rsub|H>f(y).
    </equation*>

    \;

    <item><em|Method of relaxation>:

    <\equation*>
      u<rsup|(l+1)>(x)=avg(u<rsup|(l)> <with|mode|text|on pixels surrounding
      <with|mode|math|x>>)
    </equation*>

    <item><em|Brownian motion:> Same formula as above holds for
    continuous-time.

    (Central Limit Theorem, path space version of it,
    <with|mode|math|W<rsub|t>\<sim\>k(x,t/2)>. Cylinder sets. Convergence in
    weak-* topology. Law of iterated logarithm. Proof of CLT: Convolution of
    densities becomes multiplication after Fourier transform. Use
    independence. Done.)

    <item><em|Feynman-Kac formula:> <with|mode|math|u<rsub|t>=<frac|1|2>\<Delta\>u>
    with IC <with|mode|math|f>.

    <\equation*>
      E(f(x+W<rsub|t>))=u(x,t)
    </equation*>

    <item>Implications on boundary regularity:

    <\itemize>
      <item><with|mode|math|u> <em|defined> by F-K is the Perron function

      <item><with|mode|math|y\<in\>\<partial\>U> is regular iff
      <with|mode|math|P(T<rsub|y>=0)=1> (BM immediately exits
      <with|mode|math|U>.)

      <item>Littlewood's crocodile

      <item>Lebesgue's thorn
    </itemize>
  </itemize>

  <subsection|Hearing the shape of a drum>

  <\itemize>
    <item><em|Spectral measure>:

    <\equation*>
      A(\<lambda\>)=<big|sum><rsub|k=1><rsup|\<infty\>>\<b-1\><rsub|\<lambda\><rsub|k>\<leqslant\>\<lambda\>>(\<lambda\>).
    </equation*>

    <item><em|Weyl's result>:

    <\equation*>
      lim<rsub|\<lambda\>\<rightarrow\>\<infty\>><frac|A(\<lambda\>)|\<lambda\><rsup|n/2>>=<frac|\|U\||(2\<pi\>)<rsup|n/2>\<Gamma\>(n/2)>.
    </equation*>

    <item><em|Kac's result:>

    <\equation*>
      lim<rsub|t\<rightarrow\>0+>(2\<pi\>t)<rsup|n/2><big|sum><rsub|k=1>e<rsup|-\<lambda\><rsub|k>t>=(2\<pi\>t)<rsup|n/2><big|int>e<rsup|-t\<lambda\>>A(\<mathd\>\<lambda\>)=\|U\|.
    </equation*>

    (Weyl<with|mode|math|\<Rightarrow\>>Kac: Integrate by parts, rescale.
    Proof of Kac: represent Green's function in terms of eigenfunctions
    somehow.)
  </itemize>

  <section|Wave equation>

  <\itemize>
    <item><with|mode|math|u<rsub|t t>=c<rsup|2>u<rsub|x x>>

    <item><em|D'Alembert's formula>:

    <\equation*>
      u(x,t)=<frac|1|2><left|[>f(x+c*t)+f(x-c*t)+<frac|1|c><big|int><rsub|x-c*t><rsup|x+c*t>g(y)\<mathd\>y<right|]>.
    </equation*>

    <item><em|Characteristics>.

    <item><em|Parallelogram identity>:

    <\equation*>
      u(top)+u(bottom)=u(left)+u(right).
    </equation*>

    <item>Good/bad BCs, Inflow/outflow. Domain of dependence. Method of
    reflection. Odd/even extension.

    <item><em|D'Alembertian>: <with|mode|math|\<box\>u\<assign\>u<rsub|t
    t>-c<rsup|2>\<Delta\>u=0.> <with|mode|math|u=f>,
    <with|mode|math|u<rsub|t>=g>.

    <item>Fourier Analysis: <with|mode|math|<wide|u|^>(\<xi\>,t)=<wide|f|^>(\<xi\>)cos(c\|\<xi\>\|t)+<wide|g|^>(\<xi\>)sin(c\|\<xi\>\|t/\|\<xi\>\|t)=<wide|f|^>(\<xi\>)cos(c\|\<xi\>\|t)+<wide|g|^>(\<xi\>)\<partial\><rsub|t>cos(c\|\<xi\>\|t)>:

    <\equation*>
      u(x,t)=<big|int><rsub|\<bbb-R\><rsup|n>>k(x-y,t)g(y)\<mathd\>y+\<partial\><rsub|t><big|int><rsub|\<bbb-R\><rsup|n>>k(x-y,t)f(y)\<mathd\>y
    </equation*>

    Needs to coincide with solution formula.

    <item>For <with|mode|math|n=3>, <with|mode|math|k=t\<cdot\><with|mode|text|uniform
    measure on <with|mode|math|{\|x\|=c*t}>>>

    <item><em|Method of Spherical means:> Observe:

    <\equation*>
      M<rsub|u>(x,r)=<superpose|<big|int>|
      -><rsub|S(x,r)>u(y)\<mathd\>S<rsub|y>
    </equation*>

    satisfies <em|Darboux's Equation:>

    <\equation*>
      \<Delta\><rsub|x>M<rsub|u>=<with|mode|text|``<with|mode|math|\<Delta\><rsub|r>>''>M<rsub|u>=<left|(>\<partial\><rsub|r
      r>-<frac|n-1|r>\<partial\><rsub|r><right|)>M<rsub|u>.
    </equation*>

    Similarly, if <with|mode|math|u> solves <with|mode|math|u<rsub|t
    t>=u<rsub|x x>>, then <with|mode|math|M<rsub|u>> solves the
    <em|Euler-Poisson-Darboux equation>:

    <\equation*>
      (M<rsub|u>)<rsub|t t>-\<Delta\><rsub|r>M<rsub|u>=0.
    </equation*>

    In 3D, this reduces the wave equation to
    <with|mode|math|\<partial\><rsub|t><rsup|2>(r*M<rsub|u>)=\<partial\><rsub|r><rsup|2>(r*M<rsub|u>)>,
    which we can solve by D'Alembert's formula for all <with|mode|math|x>.
    Then

    <\equation*>
      u=lim<rsub|r\<rightarrow\>0><frac|M<rsub|u>|r>.
    </equation*>

    <item>

    <\equation*>
      <frac|1|(2\<pi\>)<rsup|n/2>><big|int><rsub|\|y\|=c*t>e<rsup|-i\<xi\>\<cdot\>y>\<mathd\>S<rsub|y>=<frac|sin(c\|\<xi\>\|t)|c\|\<xi\>\|>.
    </equation*>

    <item><em|Huygens' principle.>

    <item><em|Hadamard's method of descent:> Treat 2D equation as 3D
    equation, independent of third coordinate.

    <item><em|General solution for odd <with|mode|math|n\<geqslant\>3>>:
    Assume <with|mode|math|u<rprime|'>(0)=0>. Define

    <\equation*>
      v(x,t)\<assign\><big|int>k(s,t)u(x,s)\<mathd\>s
    </equation*>

    as a temporal heat kernel average. Oddly,
    <with|mode|math|\<partial\><rsub|t>v=\<Delta\><rsub|x>v>. Solve this.
    Rewrite using spherical means. Change variables as
    <with|mode|math|\<lambda\>=1/4t> and invert using the Laplace transform

    <\equation*>
      h<rsup|#>(\<lambda\>)=<big|int><rsub|0><rsup|\<infty\>>e<rsup|-\<lambda\>\<varphi\>>h(\<varphi\>)\<mathd\>\<varphi\>.
    </equation*>

    <item>Uniqueness by energy norm.
  </itemize>

  <section|Distributions/Fourier Transform>

  <with|mode|math|U\<subset\>\<bbb-R\><rsup|n>> open

  <\itemize>
    <item><with|mode|math|\<cal-D\>(U)\<assign\>C<rsup|\<infty\>><rsub|c>(U)>.
    <with|mode|math|\<varphi\><rsub|k>\<rightarrow\>\<varphi\>> iff

    <\itemize>
      <item><with|mode|math|\<exists\>> fixed compact set <with|mode|math|F>:
      <with|mode|math|supp(\<varphi\><rsub|k>)\<subset\>F>

      <item><with|mode|math|\<forall\>\<alpha\>:><with|mode|math|sup<rsub|F>\|\<partial\><rsup|\<alpha\>>\<varphi\><rsub|k>-\<partial\><rsup|\<alpha\>>\<varphi\>\|\<rightarrow\>0>.
    </itemize>

    <item><em|Distribution:> <with|mode|math|\<cal-D\><rprime|'>(U)>

    <\itemize>
      <item>Convergence: <with|mode|math|L<rsub|k><above|\<rightarrow\>|\<cal-D\>>L><with|mode|math|\<Leftrightarrow\>><with|mode|math|\<forall\>\<varphi\>\<in\>\<cal-D\>(U):><with|mode|math|<ip|L<rsub|k>|\<varphi\>||>\<rightarrow\><ip|L|\<varphi\>||>>.
    </itemize>

    <item>Examples: <with|mode|math|L<rsup|p><rsub|loc>\<subset\>\<cal-D\><rprime|'>(U)>.
    Aside: <with|mode|math|L<rsup|p><rsub|loc>\<subset\>L<rsup|q><rsub|loc>>
    for <with|mode|math|p\<geqslant\>q>. (not for
    <with|mode|math|L<rsup|p>>), <em|Radon measure> (A Borel measure that is
    finite on compact sets.), <with|mode|math|\<delta\>> function, Cauchy
    Principal value.

    <item><em|Derivative:> <with|mode|math|<ip|\<partial\><rsup|\<alpha\>>L|\<varphi\>||>=(-1)<rsup|\|\<alpha\>\|><ip|L|\<partial\><rsup|\<alpha\>>\<varphi\>||>>.

    <item>Differentiation is continuous.

    <item><em|Partial differential operator>:
    <with|mode|math|P=<big|sum><rsub|\|\<alpha\>\|\<leqslant\>N>c<rsub|\<alpha\>>(x)\<partial\><rsup|\<alpha\>>>,
    adjoint, fundamental solution: <with|mode|math|P*K=\<delta\>>.

    <item><em|Schwartz class:> <with|mode|math|\<cal-S\>(\<bbb-R\><rsup|n>)\<subset\>C<rsup|\<infty\>>(\<bbb-R\><rsup|n>)>

    <\equation*>
      <norm|\<varphi\>|\<alpha\>,\<beta\>|>\<assign\>sup<rsub|x>\|x<rsup|\<alpha\>>\<partial\><rsup|\<beta\>>\<varphi\>(x)\|\<less\>\<infty\><space|1em>\<forall\>\<alpha\>,\<beta\>.
    </equation*>

    A polynormed, metrizable space (Use <with|mode|math|<big|sum>2<rsup|-k><big|sum><rsub|\|\<alpha\>\|+\|\<beta\>\|=k><frac|<norm|\<cdot\>|<rsub|\<alpha\>,\<beta\>>|>|1+<norm|\<cdot\>|\<alpha\>,\<beta\>|>>>).
    Complete, too. (Arzelà-Ascoli).

    <item>Examples:

    <\itemize>
      <item><with|mode|math|\<cal-D\>\<subset\>\<cal-S\>> (convergence
      carries over, too.)\ 

      <item><with|mode|math|exp(-\|x\|<rsup|2>)\<in\>\<cal-S\>>, but not
      <with|mode|math|\<in\>\<cal-D\>>.
    </itemize>

    <item><em|Fourier Transform:>

    <\equation*>
      <wide|\<varphi\>|^>(\<xi\>)=\<cal-F\>\<varphi\>(\<xi\>)=<frac|1|(2\<pi\>)<rsup|n/2>><big|int><rsub|\<bbb-R\><rsup|n>>e<rsup|-i*\<xi\>\<cdot\>x>\<varphi\>(x)\<mathd\>x
    </equation*>

    <item>Basic estimates:

    <\eqnarray*>
      <tformat|<table|<row|<cell|<norm|<wide|\<varphi\>|^>(\<xi\>)|L<rsup|\<infty\>>|>>|<cell|\<leqslant\>>|<cell|C\<\|\|\>(1+\|x\|)<rsup|n+1>\<varphi\>(x)\<\|\|\><rsub|L<rsup|\<infty\>>>\<leqslant\>C<norm|\<varphi\>|L<rsup|1>|>\<less\>\<infty\>,>>|<row|<cell|\<\|\|\>\<partial\><rsup|\<beta\>><rsub|\<xi\>><wide|\<varphi\>|^>(\<xi\>)\<\|\|\><rsub|L<rsup|\<infty\>>>>|<cell|\<leqslant\>>|<cell|C\<\|\|\>(1+\|x\|)<rsup|n+1>x<rsup|\<beta\>>\<varphi\>\<\|\|\><rsub|L<rsup|\<infty\>>>>>|<row|<cell|\<\|\|\>\<xi\><rsup|\<alpha\>><wide|\<varphi\>|^>(\<xi\>)\<\|\|\><rsub|L<rsup|\<infty\>>>>|<cell|\<leqslant\>>|<cell|C\<\|\|\>(1+\|x\|)<rsup|n+1>\<partial\><rsup|\<alpha\>><rsub|x>\<varphi\>\<\|\|\><rsub|L<rsup|\<infty\>>>>>|<row|<cell|\<\|\|\><wide|\<varphi\>|^>\<\|\|\><rsub|\<alpha\>,\<beta\>>>|<cell|\<leqslant\>>|<cell|C\<\|\|\>(1+\|x\|)<rsup|n+1>x<rsup|\<beta\>>\<partial\><rsub|x><rsup|\<alpha\>>\<varphi\>\<\|\|\><rsub|L<rsup|\<infty\>>><space|1em>\<Rightarrow\><space|1em><wide|\<varphi\>|^>\<in\>C<rsup|\<infty\>>.>>>>
    </eqnarray*>

    <item><em|Dilation:> <with|mode|math|\<sigma\><rsub|\<lambda\>>\<varphi\>(x)=\<varphi\>(x/\<lambda\>)>.
    <with|mode|math|(\<cal-F\>\<sigma\><rsub|\<lambda\>>\<varphi\>)=\<lambda\><rsup|n>\<sigma\><rsub|1/\<lambda\>>\<cal-F\>\<varphi\>>.

    <item><em|Translation:> <with|mode|math|\<tau\><rsub|h>\<varphi\>(x)=\<varphi\>(x-h)>.
    <with|mode|math|(\<cal-F\>\<tau\><rsub|h>\<varphi\>)=e<rsup|-i*h\<cdot\>\<xi\>>\<cal-F\>\<varphi\>>.

    <item><em|Inversion formula:>

    <\equation*>
      \<varphi\>(x)=<frac|1|(2\<pi\>)<rsup|n/2>><big|int><rsub|\<bbb-R\><rsup|n>>e<rsup|i*x\<cdot\>\<xi\>><wide|\<varphi\>|^>(\<xi\>)\<mathd\>\<xi\>=\<cal-F\><rsup|\<ast\>><wide|\<varphi\>|^>=\<cal-F\>\<cal-R\><wide|\<varphi\>|^>,
    </equation*>

    where <with|mode|math|\<cal-R\>\<varphi\>(x)=\<varphi\>(-x)>.

    <with|mode|math|\<cal-F\>> is an isomorphism of
    <with|mode|math|\<cal-S\>>, with <with|mode|math|\<cal-F\>\<cal-F\><rsup|\<ast\>>=Id>.

    Proof: Prove \ <with|mode|math|(\<cal-F\>\<cal-F\><rsup|\<ast\>>-Id)e<rsup|-\|x\|<rsup|2>>=0>,
    then for dilations and translations, linear comb. of which are dense in
    <with|mode|math|\<cal-S\>>. <with|mode|math|\<cal-F\>> is 1-1,
    <with|mode|math|\<cal-F\><rsup|\<ast\>>> is onto, but
    <with|mode|math|\<cal-F\><rsup|\<ast\>>=\<cal-R\>\<cal-F\>>.

    <item><with|mode|math|\<cal-F\>> isometry of <with|mode|math|L<rsup|2>>,
    <with|mode|math|\<cal-F\>> continuous from <with|mode|math|L<rsup|p>> to
    <with|mode|math|L<rsup|q>>, where\ 

    <\equation*>
      <frac|1|p>+<frac|1|q>=1,<space|1em>p\<in\>[1,2].
    </equation*>

    In particular <with|mode|math|p=1>, <with|mode|math|q=\<infty\>>.

    Proof: Show <with|mode|math|\<cal-S\>> dense in
    <with|mode|math|L<rsup|p>> (see below), extend
    <with|mode|math|\<cal-F\>>, use Plancherel for
    <with|mode|math|L<rsup|2>>.

    <item><em|Mollifier:> <with|mode|math|\<eta\>\<in\>C<rsup|\<infty\>><rsub|c>>.<with|mode|math|<big|int>\<eta\>=1>.
    <with|mode|math|\<eta\><rsub|N>(x)\<assign\>N<rsup|n>\<eta\>(N*x)>.

    <item><with|mode|math|C<rsup|\<infty\>><rsub|c>(\<bbb-R\><rsup|n>)> is
    dense in <with|mode|math|L<rsup|p>(\<bbb-R\><rsup|n>)>.
    (<with|mode|math|1\<leqslant\>p\<less\>\<infty\>>)

    Proof: <with|mode|math|<norm|\<eta\><rsub|N>\<ast\>f-f|L<rsup|p>|>\<rightarrow\>0>
    holds for step functions. Step functions are dense in
    <with|mode|math|L<rsup|p>(\<bbb-R\><rsup|n>)>.\ 

    <with|mode|math|<norm|f\<ast\>\<eta\><rsub|N>|L<rsup|p>|>\<leqslant\>C<norm|f|L<rsup|p>|>>
    (Young's)

    Pick <with|mode|math|g> a step function such that
    <with|mode|math|<norm|f-g|L<rsup|p>|>\<less\>\<varepsilon\>>. Now measure

    <\equation*>
      <norm|f\<ast\>\<eta\><rsub|N>-f|L<rsup|p>|>=<norm|f\<ast\>\<eta\><rsub|N>-g\<ast\>\<eta\><rsub|N>+g\<ast\>\<eta\><rsub|N>-g+g-f|L<rsup|p>|>.
    </equation*>

    <item><with|mode|math|C<rsup|\<infty\>><rsub|c>(\<bbb-R\><rsup|n>)> is
    dense in <with|mode|math|\<cal-S\>>.

    Proof: Smooth cutoff.

    <item><em|Plancherel's Theorem>: <with|mode|math|<ip|\<cal-F\>f|\<cal-F\>g|L<rsup|2>|>=<ip|f|g|L<rsup|2>|>.>

    Proof: by Fubini.

    <item><with|mode|math|\<cal-F\>:L<rsup|1>(\<bbb-R\><rsup|n>)\<rightarrow\><wide|C|\<dot\>>(\<bbb-R\><rsup|n>)>,
    with <with|mode|math|<wide|C|\<dot\>>\<assign\>{h:\<bbb-R\><rsup|n>\<rightarrow\>\<bbb-R\>:
    h(x)\<rightarrow\>0 (x\<rightarrow\>\<infty\>)}>.

    Proof: <with|mode|math|\<cal-S\>> is dense in <with|mode|math|L<rsup|1>>.
    Well-defined: Take <with|mode|math|\<varphi\><rsub|k>,\<psi\><rsub|k>\<rightarrow\>f\<in\>L<rsup|1>>,
    show <with|mode|math|\<cal-F\>\<varphi\><rsub|k>-\<cal-F\>\<psi\><rsub|k>\<rightarrow\>0>
    in <with|mode|math|L<rsup|\<infty\>>>.\ 

    Goes to <with|mode|math|<wide|C|\<dot\>>>: <with|color|red|unproven>.

    <item><em|Linear operator of type <with|mode|math|(r,s)>>:

    <\equation*>
      <norm|K\<varphi\>|L<rsup|s>|>\<leqslant\>C(r,s)<norm|\<varphi\>|L<rsup|r>|>.
    </equation*>

    <with|mode|math|\<cal-F\>> is of type <with|mode|math|(1,\<infty\>)> and
    <with|mode|math|(2,2)>.

    <item><em|Riesz-Thorin Convexity Theorem:> <with|mode|math|\<cal-F\>> of
    type <with|mode|math|(r<rsub|0>,s<rsub|0>)> and
    <with|mode|math|(r<rsub|1>,s<rsub|1>)>

    <\eqnarray*>
      <tformat|<table|<row|<cell|<frac|1|r>>|<cell|=>|<cell|<frac|\<theta\>|r<rsub|0>>+<frac|1-\<theta\>|r<rsub|1>>>>|<row|<cell|<frac|1|s>>|<cell|=>|<cell|<frac|\<theta\>|s<rsub|0>>+<frac|1-\<theta\>|s<rsub|1>>>>>>
    </eqnarray*>

    Then <with|mode|math|\<cal-F\>> of type <with|mode|math|(r,s)> for
    <with|mode|math|\<theta\>\<in\>[0,1]>.
  </itemize>

  <subsection|Tempered Distributions>

  <\itemize>
    <item><em|Tempered Distributions:> <with|mode|math|\<cal-S\><rprime|'>>,
    convergence as in <with|mode|math|\<cal-D\><rprime|'>>.
    <with|mode|math|\<cal-D\>\<subset\>\<cal-S\>\<subset\>\<cal-S\><rprime|'>\<subset\>\<cal-D\><rprime|'>>.

    Examples: <with|mode|math|L<rsup|1>> functions,
    <with|mode|math|e<rsup|\|x\|<rsup|2>>> not,
    <with|mode|math|e<rsup|-\|x\|<rsup|2>>>,
    <with|mode|math|<norm|(1+\|x\|<rsup|2>)<rsup|-M>f|L<rsup|1>|>\<less\>\<infty\>>.

    <item>A tempered distribution is no worse than a certain derivative
    coupled with a monomial multiplication.

    <with|mode|math|L\<in\>\<cal-S\><rprime|'>><with|mode|math|\<Rightarrow\>><with|mode|math|\<exists\>C,N\<forall\>\<varphi\>\<in\>\<cal-S\>:><with|mode|math|\|<ip|L|\<varphi\>||>\|\<leqslant\><big|sum><rsub|\|\<alpha\>\|,\|\<beta\>\|\<leqslant\>N><norm|x<rsup|\<alpha\>>\<partial\><rsup|\<beta\>>\<varphi\>|L<rsup|\<infty\>>|>>
    (continuity).

    <item><with|mode|math|<ip|\<eta\>\<ast\>L|\<varphi\>||>=<ip|L|(\<cal-R\>\<eta\>)\<ast\>\<varphi\>||>>
    for <with|mode|math|L\<in\>D<rprime|'>>, <with|mode|math|\<cal-R\>> is
    reflection and <with|mode|math|\<eta\>> a mollifier

    <item><with|mode|math|\<eta\>\<ast\>L> is a
    <with|mode|math|C<rsup|\<infty\>>> function, namely
    <with|mode|math|\<gamma\>(x)=<ip|L|\<tau\><rsub|x>\<cal-R\>\<eta\>||>>,
    where <with|mode|math|\<tau\><rsub|x>f(y)=(y-x)>.

    Proof: 1. <with|mode|math|\<gamma\>> maps to <with|mode|math|\<bbb-R\>>.
    2. <with|mode|math|\<gamma\>> sequentially continuous. 3.
    <with|mode|math|\<gamma\>\<in\>C<rsup|1>> (FD). 4.
    <with|mode|math|\<gamma\>\<in\>C<rsup|\<infty\>>> (induction). 5.
    <with|mode|math|<ip|\<eta\>\<ast\>L|\<varphi\>||>=<ip|\<gamma\>|\<varphi\>||>>
    (Riemann sums).

    <item><with|mode|math|\<cal-D\>> is dense in
    <with|mode|math|\<cal-D\><rprime|'>>.

    Proof: <with|mode|math|\<chi\><rsub|m>\<assign\>\<b-1\><rsub|[-m,m]>>.
    Fix <with|mode|math|L\<in\>\<cal-D\><rprime|'>>,
    <with|mode|math|L<rsub|m>\<assign\>\<chi\><rsub|m>(\<eta\><rsub|m>\<ast\>L)\<in\>D\<rightarrow\>L>
    in <with|mode|math|\<cal-D\><rprime|'>>.\ 

    <item><with|mode|math|\<cal-S\>> is dense in
    <with|mode|math|\<cal-S\><rprime|'>>.

    (because <with|mode|math|\<cal-D\>> is already dense in
    <with|mode|math|\<cal-D\><rprime|'>>.)

    <item><em|Transpose> <with|mode|math|K<rsup|t>:\<cal-S\>\<rightarrow\>\<cal-S\>>
    for <with|mode|math|K:\<cal-S\>\<rightarrow\>\<cal-S\>> as by
    <with|mode|math|<ip|K<rsup|t>L|\<varphi\>||>\<assign\><ip|L|K\<varphi\>||>>.

    <item><with|mode|math|K:\<cal-S\>\<rightarrow\>\<cal-S\>> linear and
    continuous. <with|mode|math|K<rsup|t>\|<rsub|\<cal-S\>>> continuous.
    <with|mode|math|\<exists\>!>unique, continuous extension of
    <with|mode|math|K<rsup|t>> onto <with|mode|math|\<cal-S\><rprime|'>>.

    <item><with|mode|math|\<cal-F\>:\<cal-S\><rprime|'>\<rightarrow\>\<cal-S\><rprime|'>>
    continuous.

    <item><with|mode|math|\<cal-F\>\<delta\>=1/(2\<pi\>)<rsup|n/2>>.

    <item><with|mode|math|0\<less\>\<beta\>\<less\>n>,
    <with|mode|math|C<rsub|\<beta\>>=\<Gamma\>((n-\<beta\>)/2)>

    <\equation*>
      \<cal-F\>(C<rsub|\<beta\>>\|x\|<rsup|-\<beta\>>)=C<rsub|n-\<beta\>>\|x\|<rsup|-(n-\<beta\>)>.
    </equation*>

    Use this to solve Laplace's equation.
  </itemize>

  <\section>
    Hyperbolic Equations
  </section>

  <\itemize>
    <item>General constant coefficient problem.
    <with|mode|math|P(D,\<tau\>)=\<tau\><rsup|m>+\<tau\><rsup|m-1>P<rsub|1>(D)+\<cdots\>+P<rsub|m>(D)>

    <item><em|Duhamel's principle:> Solve <with|mode|math|P(D,\<tau\>)u=f> by
    solving the <em|standard problem> <with|mode|math|P(D,\<tau\>)u<rsub|s>=0>,
    <with|mode|math|u<rsub|s>(0)=0>, <with|mode|math|\<partial\><rsub|t><rsup|m-1>u<rsub|s>(0)=g>
    and finding

    <\equation*>
      u(x,t)=<big|int><rsub|0><rsup|t>u<rsub|s>\<mathd\>s.
    </equation*>

    <item>Treat remaining ICs by solving standard problems for
    <with|mode|math|\<tau\><rsup|m-1>P<rsub|1>,\<ldots\>,\<tau\><rsup|0>P<rsub|m>>,
    each time adding to the right hand side, which can finally be killed with
    the above approach.

    <item>Fourier-transforms to <with|mode|math|P(i\<xi\>,\<tau\>)<wide|u|^>=0>,
    with <with|mode|math|\<tau\>=\<partial\><rsub|t>>.

    Initial conditions <with|mode|math|\<tau\><rsup|0\<ldots\>m-2><wide|u|^>(\<xi\>,0)>,
    <with|mode|math|\<tau\><rsup|m-1><wide|u|^>(\<xi\>,0)>.

    <item>Representation of the solution:

    <\eqnarray*>
      <tformat|<table|<row|<cell|Z(\<xi\>,t)>|<cell|=>|<cell|<frac|1|2\<pi\>><big|int><rsub|\<Gamma\>><frac|e<rsup|i\<lambda\>t>|P(i\<xi\>,i\<lambda\>)>\<mathd\>\<lambda\>>>|<row|<cell|P(i\<xi\>,\<tau\>)Z>|<cell|=>|<cell|<frac|1|2\<pi\>><big|int><rsub|\<Gamma\>>P(i\<xi\>,i\<lambda\>)<frac|e<rsup|i\<lambda\>t>|P(i\<xi\>,i\<lambda\>)>\<mathd\>\<lambda\>=<frac|1|2\<pi\>><big|int><rsub|\<Gamma\>>e<rsup|i\<lambda\>t>\<mathd\>\<lambda\>=0,>>>>
    </eqnarray*>

    where <with|mode|math|\<Gamma\>> is a path around the roots.

    <item>Classical solution requires <with|mode|math|u\<in\>C<rsup|m>><with|mode|math|>.
    Requires <with|mode|math|\<forall\>T\<exists\>C<rsub|T>,N:>

    <\equation*>
      \|\<tau\><rsup|k>Z(\<xi\>,t)\|\<leqslant\>C<rsub|T>(1+\|\<xi\>\|)<rsup|N>.
    </equation*>

    <item><em|Hyperbolicity:> A standard problem is
    <em|hyperbolic><with|mode|math|:\<Leftrightarrow\>><with|mode|math|\<exists\>>a
    <with|mode|math|C<rsup|m>> solution for all
    <with|mode|math|g\<in\>\<cal-S\>(\<bbb-R\><rsup|n>)>.

    <item><em|Gårding's Criterion:> It's hyperbolic iff
    <with|mode|math|\<exists\>c\<in\>\<bbb-R\>:><with|mode|math|P(i\<xi\>,i\<lambda\>)\<neq\>0>
    for all <with|mode|math|\<xi\>> and <with|mode|math|Im
    \<lambda\>\<leqslant\>-c>.

    Proof: Estimate around in the above representation for
    <with|mode|math|Z>.

    <item><em|Paley-Wiener Theorem:> <with|mode|math|g\<in\>L<rsup|1>><with|mode|math|\<Rightarrow\>><with|mode|math|<wide|g|^>>
    entire.
  </itemize>

  <section|Conservation Laws>

  <\itemize>
    <item><with|mode|math|u<rsub|t>+f(u)<rsub|x>=0>.

    Why are they called called conservation laws?

    <\equation*>
      <frac|\<mathd\>|\<mathd\>t><big|int>u=<big|int><rsub|>u<rsub|t>=<big|int>f(u)<rsub|x>=f(b)-f(a)\<rightarrow\>0.
    </equation*>

    <item><em|Inviscid Burgers' Equation:>
    <with|mode|math|u<rsub|t>+(u<rsup|2>)<rsub|x>=0>.

    <item><em|Characteristics:> Assume <with|mode|math|u=u(x(t),t)>,

    <\equation*>
      <frac|\<mathd\>u|\<mathd\>t>=<frac|\<partial\>u|\<partial\>x>\<cdot\><frac|\<mathd\>x|\<mathd\>t>+<frac|\<partial\>u|\<partial\>t>
    </equation*>

    Compare shape with

    <\equation*>
      0=u<rsub|x>f<rprime|'>(u)+u<rsub|t>,
    </equation*>

    obtain <with|mode|math|\<mathd\>x/\<mathd\>t=f<rprime|'>(u)>.

    <item><em|Weak solution>: slap test function onto equation, integrate by
    parts.

    <item><em|Rankine-Hugoniot>:

    <\equation*>
      <with|mode|text|shock speed>=<frac|<jump|f(u)>|<jump|u>>
    </equation*>

    Apply weak solution formula across a jump, consider normal geometrically
    to obtain speed.

    <item><em|Riemann problem>: Jump IC. <with|mode|math|\<rightarrow\>>
    non-uniqueness of the weak solution for jump up: rarefaction wave or
    shock with correct speed?

    <item><em|Hopf's treatment of Burger's Equation:>

    <\itemize>
      <item>Add viscosity to get <with|mode|math|u<rsub|t>+(u<rsup|2>/2)<rsub|x>=\<varepsilon\>u<rsub|x
      x>>.

      <item>Put <with|mode|math|U> as an antiderivative of
      <with|mode|math|u>.

      <item>Gives Hamilton-Jacobi PDE <with|mode|math|U<rsub|t>+U<rsub|x><rsup|2>/2=\<varepsilon\>U<rsub|x
      x>>.

      <item>Now try to rewrite that into a linear equation, by assuming
      <with|mode|math|\<psi\>=\<psi\>(u)>. Yields ODE
      <with|mode|math|C\<varphi\><rprime|''>+C\<varphi\><rprime|'>=0>,
      solution <with|mode|math|\<psi\>=exp(-U/2\<varepsilon\>)>.

      <item>This gives the heat equation <with|mode|math|\<psi\><rsub|t>=\<varepsilon\>\<psi\><rsub|x
      x>>.

      <item>

      <\equation*>
        u=2\<varepsilon\><frac|\<psi\><rsub|x>|\<psi\>>=<frac|<big|int><frac|x-y|t>exp(-G/2\<varepsilon\>)\<mathd\>y|<big|int>exp(-G/2\<varepsilon\>)\<mathd\>y>=<frac|x|t>-<frac|\<langle\>y\<rangle\>|t>\<rightarrow\><frac|x|t>-<frac|argmin
        G|t>
      </equation*>

      with <with|mode|math|G=(x-y)<rsup|2>/2t+U<rsub|0>>.
    </itemize>

    <item><with|mode|math|a<rsub|->=inf argmin G>,
    <with|mode|math|a<rsub|+>=sup argmin G>.

    <item><em|Properties:> well-defined, increasing,
    <with|mode|math|a<rsub|+>(\<leftarrow\>)\<leqslant\>a<rsub|->(\<rightarrow\>)>,
    <with|mode|math|a<rsub|->> left-continuous, <with|mode|math|a<rsub|+>>
    right-continuous, go to <with|mode|math|\<pm\>\<infty\>>. Equal except
    for a countable set of shocks.

    <item><em|Hopf's theorem:>

    <\equation*>
      <frac|x-a<rsub|+>|t>\<leqslant\>liminf<rsub|\<varepsilon\>\<rightarrow\>0>u<rsup|\<varepsilon\>>\<leqslant\>limsup<rsub|\<varepsilon\>\<rightarrow\>0>u<rsup|\<varepsilon\>>\<leqslant\><frac|x-a<rsub|->|t>
    </equation*>

    <item><with|mode|math|u<rsub|0>\<in\>BC> (bounded,
    continuous)<with|mode|math|\<Rightarrow\>><with|mode|math|u(\<cdot\>,t)\<in\>BV<rsub|loc>>.
    <with|color|red|Globally BV?>

    Proof: <with|mode|math|x,a<rsub|+>,a<rsub|->> are
    increasing<with|mode|math|\<Rightarrow\>>differences in
    <with|mode|math|BV<rsub|loc>>.

    <item>Vanishing viscosity solutions are weak solutions.

    Proof: Pass to vanishing viscosity under integral using DCT and
    boundedness.

    <item>Cole-Hopf solutions produce rarefaction <with|mode|math|x/t> for
    jump up, shock for jump down.

    <item>More properties:

    <\itemize>
      <item><with|mode|math|lim<rsub|\<varepsilon\>\<rightarrow\>0>u<rsup|\<varepsilon\>>>
      exists except for a countable set. <with|mode|math|u=lim
      u<rsup|\<varepsilon\>>\<in\>BV<rsub|loc>> with left and right limits.

      Proof: <with|mode|math|u> is a difference of increasing functions.

      <item><em|Lax-Oleinik entropy condition:>
      <with|mode|math|u(x<rsub|->,t)\<gtr\>u(x<rsub|+>,t)> at jumps.\ 

      ``Characteristics never leave a shock.''

      Proof: Travelling waves for Burgers with viscosity only exist for
      <with|mode|math|u<rsub|->\<gtr\>u<rsub|+>>.

      <item><with|mode|math|x> a shock location:

      <\eqnarray*>
        <tformat|<table|<row|<cell|<with|mode|text|shock
        speed>>|<cell|=>|<cell|<frac|<jump|f(u)>|u>=<frac|1|2>(u<rsub|+>+u<rsub|->),>>|<row|<cell|<with|mode|text|shock
        speed>=<frac|1|2><left|(>u(x<rsub|->,t)+u(x<rsub|+>,t)<right|)>>|<cell|=>|<cell|<superpose|<big|int>|
        -><rsub|a<rsub|->><rsup|a<rsub|+>>u<rsub|0>(y)\<mathd\>y,>>|<row|<cell|(a<rsub|+>-a<rsub|->)<with|mode|text|shock
        speed>>|<cell|=>|<cell|<big|int><rsub|a<rsub|->><rsup|a<rsub|+>>u<rsub|0>(y)\<mathd\>y>>>>
      </eqnarray*>

      The last equation here is a momentum conservation equality.

      Proof: <with|mode|math|G(a<rsup|+>)=G(a<rsup|->)>.
    </itemize>

    <item><em|Entropy/entropy-flux pair>:
    <with|mode|math|\<Phi\>,\<Psi\>:\<bbb-R\><rsup|m>\<rightarrow\>\<bbb-R\>>
    smooth are an e/ef pair for <with|mode|math|u<rsub|t>+f(u)<rsub|x>=0><with|mode|math|:\<Leftrightarrow\>><with|mode|math|\<Phi\>>
    convex, <with|mode|math|\<Phi\><rprime|'>f<rprime|'>=\<Psi\><rprime|'>>.
    Then <with|mode|math|\<Phi\>(u)<rsub|t>+\<Psi\>(u)<rsub|x>=0> for
    perfectly smooth solutions, otherwise
    <with|mode|math|\<Phi\>(u<rsub|t>)+\<Psi\>(u)<rsub|x>\<leqslant\>0> in
    the distributional sense, which means

    <\equation*>
      <big|int><rsub|0><rsup|\<infty\>><big|int><rsub|-\<infty\>><rsup|\<infty\>>\<Phi\>(u)v<rsub|t>+\<Psi\>(u)v<rsub|x>\<mathd\>x\<mathd\>t\<geqslant\>0.
    </equation*>

    for smooth non-negative <with|mode|math|v>.

    <item>By the vanishing viscosity method, we get an entropy solution.

    Proof: Multiply the viscosity-added c.law by
    <with|mode|math|\<Phi\><rprime|'>>. Use chain rule on
    <with|mode|math|\<Phi\>(u<rsup|\<varepsilon\>>)<rsub|x x>>. Use convexity
    of <with|mode|math|\<Phi\>> to show one term involving
    <with|mode|math|\<varphi\><rprime|''>> non-negative. Multiply by a
    non-negative smooth function, let <with|mode|math|\<varepsilon\>\<rightarrow\>0>
    to obtain entropy inequality.

    <item><em|Entropy solution>: <with|mode|math|u> is an entropy solution of
    a c.law if <with|mode|math|u> is a weak solution that satisfies the
    entropy condition for every e/ef pair.

    <item><em|Dissipation measure:>

    <\eqnarray*>
      <tformat|<table|<row|<cell|<frac|\<mathd\>|\<mathd\>t><big|int>(u<rsup|\<varepsilon\>>)<rsup|2>>|<cell|=>|<cell|-2\<varepsilon\><big|int>(u<rsup|\<varepsilon\>><rsub|x>)<rsup|2>.>>>>
    </eqnarray*>

    Assuming a traveling wave solution of the form

    <\equation*>
      u<rsup|\<varepsilon\>>=v<left|(><frac|x-c*t|\<varepsilon\>><right|)>,
    </equation*>

    we find

    <\equation*>
      <frac|\<mathd\>|\<mathd\>t><big|int>(u<rsup|\<varepsilon\>>)<rsup|2>=<frac|(u<rsub|->-u<rsub|+>)<rsup|3>|6>.
    </equation*>

    <item><em|Kruºkov's Uniqueness Theorem:>
    <with|mode|math|L<rsup|\<infty\>>> Entropy solutions <with|mode|math|u>,
    <with|mode|math|v>, <with|mode|math|S<rsub|t>> cuts of the event cone
    (given by max. speed <with|mode|math|c<rsup|\<ast\>>=max<rsub|range u>
    \|f<rprime|'>\|>. Then for <with|mode|math|t<rsub|1>\<less\>t<rsub|2>>

    <\equation*>
      <big|int><rsub|S<rsub|t<rsub|2>>>\|u-v\|\<leqslant\><big|int><rsub|S<rsub|t<rsub|1>>>\|u-v\|.
    </equation*>

    Proof: Doubling trick, clever choice of test functions.

    Implies uniqueness.
  </itemize>

  <section|Hamilton-Jacobi Equations>

  <\itemize>
    <item><with|mode|math|u<rsub|t>+H(D u,x)=0>.

    <item>Example: Curve evolving with normal velocity:
    <with|mode|math|u<rsub|t>+<sqrt|1+\|D<rsub|x>u\|<rsup|2>>=0>.

    <item>Non-Example: Motion by mean curvature
    <with|mode|math|u<rsub|t>=u<rsub|x x>/(1+u<rsub|x>)<rsup|2>> (parabolic).

    <item>Example: Substitute <with|mode|math|U=<big|int>u> in conservation
    laws.

    <item>PDE is infinitely-many-particle limit of Hamilton ODE

    <\eqnarray*>
      <tformat|<table|<row|<cell|<wide|x|\<dot\>>>|<cell|=>|<cell|\<partial\><rsub|p>H(p,x)>>|<row|<cell|<wide|p|\<dot\>>>|<cell|=>|<cell|-\<partial\><rsub|x>H(p,x),>>>>
    </eqnarray*>

    which coincides with characteristic equation of PDE.

    <item>Mechanics motivation:

    <\itemize>
      <item><with|mode|math|L(q,x)=T-V>

      <item>Lagrange's Equation

      <\equation*>
        <frac|\<mathd\>|\<mathd\>t><left|(><frac|\<partial\>L|\<partial\>q><right|)>=<frac|\<partial\>L|\<partial\>x>.
      </equation*>

      Way to see this: If <with|mode|math|RHS=0>, then <with|mode|math|L>
      symmetric in <with|mode|math|x>, so LHS becomes conserved. (Noether's
      theorem.)

      Equivalent to Hamilton's ODE (Proof:
      <with|mode|math|H=max<rsub|q>(q*p-L(x,q,t))>, where
      <with|mode|math|q=q(x,p,t)> is the solution of
      <with|mode|math|p=\<partial\><rsub|q>L(v)>.

      <item>Action, given path <with|mode|math|x(t)>:

      <\equation*>
        S(x)=<big|int><rsub|0><rsup|t>L(<wide|x|\<dot\>>,x,t)\<mathd\>t
      </equation*>

      <item>Principle of least action: <with|mode|math|min
      S><with|mode|math|\<Leftrightarrow\>>Lagrange's Equation.

      Proof: <with|mode|math|u+\<varepsilon\>v>, derivative by
      <with|mode|math|\<varepsilon\>>, the usual.

      <item>Generalized momentum: <with|mode|math|p=\<partial\><rsub|q>L>.
      Assumed solvable for <with|mode|math|q>.

      <item>Hamiltonian: <with|mode|math|H=T+V=p\<cdot\>q-L=2T-(T-V)=T+V>.

      <item><em|Legendre transform:> More general way of obtaining
      <with|mode|math|H>. Assume <with|mode|math|L(q)> (dropping
      dependencies!) convex, <with|mode|math|lim<rsub|\|q\|\<rightarrow\>\<infty\>>L(q)/\|q\|=\<infty\>>.
      Then

      <\equation*>
        H(p)=L<rsup|\<ast\>>(p)=sup<rsub|q>{p\<cdot\>q-L(q)}.
      </equation*>

      Solved when <with|mode|math|p=\<partial\><rsub|q>L>, but in a more
      general sense.\ 

      Duality: Edge<with|mode|math|\<leftrightarrow\>>Corner.
      Subdifferentials.

      <item><with|mode|math|L> convex<with|mode|math|\<Rightarrow\>><with|mode|math|L<rsup|\<ast\>\<ast\>>=L>.

      Proof: Prove convexity and superlinearity of
      <with|mode|math|L<rsup|\<ast\>>>. Use symmetry

      <\equation*>
        H(p)+L(q)\<geqslant\>p\<cdot\>q
      </equation*>

      to prove two sides of the equality <with|mode|math|H<rsup|\<ast\>>=L>.
    </itemize>

    <item><em|Hopf-Lax formula>: <with|mode|math|g> is IC

    <\equation*>
      u(x,t)=inf<left|{><big|int>L(<wide|x|\<dot\>>)\<mathd\>x+g(y),x(0)=y,x(t)=x<right|}>=min<left|{>t*L<left|(><frac|x-y|t><right|)>+g(y)<right|}>.
    </equation*>

    Proof: Inf bounded above by straight-line characteristic. Lower bound
    works by Jensen's inequality.

    <item><em|Semigroup Property>.

    Proof: Always pick particular solutions, prove both sides of the
    inequality.

    <item><with|mode|math|u> defined by Hopf-Lax is Lipschitz if
    <with|mode|math|g> is Lipschitz.

    Proof: Lipschitzicity for given <with|mode|math|t> is immediate (pick
    good <with|mode|math|z>). Transform problem to comparison with
    <with|mode|math|t=0> by semigroup property. Temporal estimate is screwy,
    involves special choices in inf.

    <item><with|mode|math|u> by Hopf-Lax is differentiable a.e. and satisfies
    the H-J PDE where it is.

    Proof: Rademacher's Theorem. Prove <with|mode|math|u<rsub|t>+H(D
    u)\<leqslant\>0> for forward in time by taking increments
    <with|mode|math|\<rightarrow\>0>, using inequality with Legendre
    transform.

    <item>Lipschitz+Differentiable solution a.e. is not sufficient for
    uniqueness. (45-degree angle trough vs. 90-degree trough)

    <item><with|mode|math|f:\<bbb-R\><rsup|n>\<rightarrow\>\<bbb-R\>>
    <em|semiconcave> if

    <\equation*>
      f(x+z)-2f(x)+f(x-z)\<leqslant\>C\|z\|<rsup|2>
    </equation*>

    for some <with|mode|math|z>.

    <with|mode|math|\<Leftrightarrow\>><with|mode|math|f(z)-C/2\|z\|<rsup|2>>
    is concave.

    <with|mode|math|\<Leftrightarrow\>>``can be forced into concavity by
    subtracting a parabola.''

    <with|mode|math|\<Leftarrow\>><with|mode|math|C<rsup|2>> and bounded
    second derivatives implies semiconcavity.

    <item><with|mode|math|g> semiconcave<with|mode|math|\<Rightarrow\>><with|mode|math|u>
    semiconcave.

    Clever choice of test locations in Hopf-Lax.

    <item><with|mode|math|H:\<bbb-R\><rsup|n>\<rightarrow\>\<bbb-R\>>
    <em|uniformly convex><with|mode|math|:\<Leftrightarrow\>>

    <\equation*>
      <big|sum><rsub|i,j>H<rsub|p<rsub|i>p<rsub|j>>\<xi\><rsub|i>\<xi\><rsub|j>\<geqslant\>j\|\<xi\>\|<rsup|2>.
    </equation*>

    <item>If <with|mode|math|H> uniformly convex. Then <with|mode|math|u> is
    semiconcave (indep. of initial data)

    Proof: Taylor, mess about with Hopf-Lax.

    <item>Now <with|mode|math|H(p)\<rightarrow\>H(p,x)> nonconvex.

    <item><em|Vanishing Viscosity Method:> Use <with|mode|math|u<rsub|t>+H(D
    u,x)=\<varepsilon\>\<Delta\>u>. Locally uniform convergence follows from
    Arzelà-Ascoli.\ 

    <item><with|mode|math|u> is a <em|viscosity
    solution><with|mode|math|:\<Leftrightarrow\>><with|mode|math|u=g> on
    <with|mode|math|\<bbb-R\><rsup|n>\<times\>{t=0}>, for each
    <with|mode|math|v\<in\>C<rsup|\<infty\>>(\<bbb-R\><rsup|n>\<times\>(0,\<infty\>))>

    <\quote-env>
      <with|mode|math|u-v> has a local maximum at
      <with|mode|math|(x<rsub|0>,t<rsub|0>)><with|mode|math|\<Rightarrow\>><with|mode|math|v<rsub|t>(x<rsub|0>,t<rsub|0>)+H(D
      v(x<rsub|0>,t<rsub|0>))\<leqslant\>0> (and
      min<with|mode|math|\<rightarrow\>><with|mode|math|\<geqslant\>>).
    </quote-env>

    <item>If <with|mode|math|u> is a <em|vanishing> viscosity solution, then
    it is a viscosity solution.

    Proof: Convergence is locally uniform as
    <with|mode|math|\<varepsilon\><rsub|j>\<rightarrow\>0>. Thus for each
    fixed ball around a local strict maximum in <with|mode|math|u-v>, a local
    maximum in <with|mode|math|u<rsup|\<varepsilon\>>-v> exists if
    <with|mode|math|\<varepsilon\>> is small enough. There,
    <with|mode|math|v<rsub|x>=u<rsub|x><rsup|\<varepsilon\>>> and
    <with|mode|math|v<rsub|t>=u<rsub|t><rsup|\<varepsilon\>>> and
    <with|mode|math|-\<Delta\>u<rsup|\<varepsilon\>>\<geqslant\>-\<Delta\>v>.
    <with|mode|math|v<rsub|t>+H(D v)\<leqslant\>0> follows. Generalize to
    non-strict maxima by adding parabolas.

    <item>A classical solution of a H-J PDE is a viscosity solution.

    Proof: Maximum of <with|mode|math|u-v><with|mode|math|\<Rightarrow\>>derivatives
    are equal<with|mode|math|\<Rightarrow\>>PDE.

    <item><em|Touching by <with|mode|math|C<rsup|1>> function>:
    <with|mode|math|u> continuous. <with|mode|math|u> differentiable at
    <with|mode|math|x<rsub|0>>. Then <with|mode|math|\<exists\>v\<in\>C<rsup|1>:v(x<rsub|0>)=u(x<rsub|0>)>,
    <with|mode|math|u-v> has a strict local max.

    <item><with|mode|math|u> viscosity solution<with|mode|math|\<Rightarrow\>><with|mode|math|u>
    satisfies H-J wherever it is differentiable

    Proof: Mollify touching function, <with|mode|math|u-v<rsup|\<varepsilon\>>>
    maintains strict max., verify definition of Viscosity solution.
    (Mollification necessary because test functions are required to be
    <with|mode|math|C<rsup|\<infty\>>>.)

    <item><em|Uniqueness:> <with|mode|math|H\<in\>Lip<rsub|p>(C)\<cap\>Lip<rsub|x>(C1+\|p\|)><with|mode|math|\<Rightarrow\>>uniqueness.

    Proof: doubling trick again.
  </itemize>

  <section|Sobolev Spaces>

  <with|mode|math|1\<leqslant\>p\<less\>\<infty\>>.

  <\itemize>
    <item><with|mode|math|<norm|u|k,p;\<Omega\>|>=<big|sum><rsub|\|\<alpha\>\|\<leqslant\>k><norm|D<rsup|\<alpha\>>u|p|>>.

    <item><with|mode|math|W<rsup|k,p>(\<Omega\>)\<assign\>{u\<in\>\<cal-D\><rprime|'>(\<Omega\>):D<rsup|\<alpha\>>u\<in\>L<rsup|p>(\<Omega\>),\|\<alpha\>\|\<leqslant\>k}>
    Banach space.

    <item><with|mode|math|W<rsup|k,p><rsub|0>(\<Omega\>)\<assign\>cl(\<cal-D\>(\<Omega\>),<norm|\<cdot\>|k,p;\<Omega\>|>)>.

    <item><with|mode|math|u\<in\>W<rsup|l,p>(\<Omega\>)>.
    <with|mode|math|\<Omega\><rprime|'>\<subset\>\<subset\>\<Omega\>>
    open<with|mode|math|\<Rightarrow\>><with|mode|math|\<exists\>u<rsub|k>\<in\>C<rsup|\<infty\>><rsub|c>(\<Omega\><rprime|'>):<norm|u<rsub|k>-u|l,p;\<Omega\><rprime|'>|>\<rightarrow\>0.>

    Proof: Mollification, throw derivatives onto <with|mode|math|u> by
    integration by parts.

    <item><with|mode|math|u\<in\>W<rsup|l,p>(\<Omega\>)>,
    <with|mode|math|\<Omega\>> bounded<with|mode|math|\<Rightarrow\>><with|mode|math|\<exists\>u<rsub|k>\<in\>C<rsup|\<infty\>>(\<Omega\>)\<cap\>W<rsup|l,p>(\<Omega\>):<norm|u<rsub|k>-u|l,p;\<Omega\>|>\<rightarrow\>0.>

    Proof: Exhaust <with|mode|math|\<Omega\>> by
    <with|mode|math|U<rsub|k>\<assign\>{dist(x,\<partial\>U)\<gtr\>1/k}>.
    Consider smooth partition of unity <with|mode|math|\<zeta\><rsub|i>>
    subordinate to <with|mode|math|V<rsub|i>\<assign\>\<Omega\><rsub|i+3>\<setminus\><wide|\<Omega\>|\<bar\>><rsub|i+1>>.
    <with|mode|math|u<rsub|i>\<assign\>\<eta\><rsub|\<varepsilon\><rsub|i>>\<ast\>(\<zeta\><rsub|i>u)>
    s.t. <with|mode|math|<norm|u<rsub|i>-\<zeta\><rsub|i>u|l,p|>\<less\>\<delta\>2<rsup|-i-1>>.
    Give one more set of wiggle room on each side for mollification.
    <with|mode|math|v\<assign\><big|sum>\<zeta\><rsub|i>u<rsub|i>\<in\>C<rsup|\<infty\>>>
    because there's only a finite number of terms for fixed point/set. Then
    estimate <with|mode|math|<norm|u-v|l,p|>>.

    <item>Typical idea: Consider

    <\equation*>
      f<rsup|\<ast\>>(x)=lim<rsub|r\<rightarrow\>0><superpose|<big|int>|
      -><rsub|B(x,r)>f(y)\<mathd\>y.
    </equation*>

    <item><with|mode|math|u\<in\>W<rsup|1,p>(\<Omega\>)>,
    <with|mode|math|\<Omega\><rprime|'>\<subset\>\<subset\>\<Omega\>>. Then

    <\itemize>
      <item>There exists a representative on
      <with|mode|math|\<Omega\><rprime|'>> that is absolutely continuous on a
      line and whose classical derivative agrees a.e. with the weak one.

      <item>If the above is true of a function, then
      <with|mode|math|u\<in\>W<rsup|1,p>(\<Omega\>)>.
    </itemize>

    Proof: WLOG <with|mode|math|p=1> (Jensen). <with|color|red|WTF?>

    Consequences: <with|mode|math|W<rsup|1,p>> closed wrt. max, min, abs.
    value, <with|mode|math|\<cdot\><rsup|+>>. <with|mode|math|\<Omega\>>
    connected, <with|mode|math|D u=0><with|mode|math|\<Rightarrow\>><with|mode|math|u>
    constant.
  </itemize>

  <subsection|Campanato>

  <\itemize>
    <item><em|Oscillation:>

    <\equation*>
      osc<rsub|U>=sup<rsub|x,y\<in\>U>\|u(x)-u(y)\|.
    </equation*>

    <item><with|mode|math|C<rsup|0,\<alpha\>>\<assign\>{\|u(x)-u(y)\|\<leqslant\>C\|x-y\|<rsup|\<alpha\>>}>.
    <with|mode|math|<norm|u|C<rsup|0,\<alpha\>>|>\<assign\><norm|u|C(<wide|U|\<bar\>>)|>+sup<rsub|x\<neq\>y>\|u(x)-u(y)\|/\|x-y\|<rsup|\<alpha\>>>.

    <item><with|mode|math|C<rsup|k,\<alpha\>>\<assign\>D<rsup|\<alpha\>>\<in\>C<rsup|0,\<alpha\>>>.
    Norm: sum over multi-indices.

    <item><em|Campanato's Inequality:> <with|mode|math|u\<in\>L<rsup|1><rsub|loc>(\<Omega\>)>,
    <with|mode|math|0\<less\>\<alpha\>\<leqslant\>1>,
    <with|mode|math|\<exists\>M\<gtr\>0>:

    <\equation*>
      <superpose|<big|int>| -><rsub|B>\|u(x)-<wide|u|\<bar\>><rsub|B>(x)\|\<mathd\>x\<leqslant\>M*r<rsup|\<alpha\>>.
    </equation*>

    Then <with|mode|math|u\<in\>C<rsup|0,\<alpha\>>(\<Omega\>)> and
    <with|mode|math|osc<rsub|B(x,r/2)>u\<leqslant\>C*M*r<rsup|\<alpha\>>.>
    <with|mode|math|<wide|u|\<bar\>><rsub|B>> is the mean over
    <with|mode|math|B>.

    Proof: <with|mode|math|x> a Lebesgue point of <with|mode|math|u>,
    <with|mode|math|B(x,r/2)\<subset\>B(z,r)>. Then
    <with|mode|math|\|<wide|u|\<bar\>><rsub|B(x,r/2)>-<wide|u|\<bar\>><rsub|B(z,r)>\|\<leqslant\>2<rsup|n>M*r<rsup|\<alpha\>>.>
    Iteration via geometric series and Lebesgue-pointy-ness yields

    <\equation*>
      \|u(x)-<wide|u|\<bar\>><rsub|B(z,r)>\|\<leqslant\>C(n,\<alpha\>)M*r<rsup|\<alpha\>>.
    </equation*>

    For two Lebesgue points,

    <\eqnarray*>
      <tformat|<table|<row|<cell|\|u(x)-u(y)\|\<leqslant\>\|u(x)-<wide|u|\<bar\>><rsub|B(z,r)>\|+\|<wide|u|\<bar\>><rsub|B(z,r)>-u(y)\|<rsub|>>|<cell|\<leqslant\>>|<cell|C(n,\<alpha\>)M*r<rsup|\<alpha\>>.>>>>
    </eqnarray*>
  </itemize>

  <subsection|Sobolev>

  <\itemize>
    <item><em|Gagliardo-Nirenberg-Sobolev>:
    <with|mode|math|u\<in\>C<rsup|1><rsub|c>(\<bbb-R\><rsup|n>)>,
    <with|mode|math|1\<leqslant\>p\<less\>n><with|mode|math|><with|mode|math|\<Rightarrow\>>

    <\equation*>
      <norm|u|p<rsup|\<ast\>>|>\<leqslant\>C<norm|D u|p|>,
    </equation*>

    where

    <\equation*>
      <frac|1|p<rsup|*\<ast\>>>+<frac|1|n>=<frac|1|p><space|1em>\<Rightarrow\><space|1em>p<rsup|\<ast\>>\<gtr\>p.
    </equation*>

    <item>Considering what happens when you scale functions
    <with|mode|math|u\<rightarrow\>u<rsub|\<lambda\>>(x)\<assign\>u(\<lambda\>x)>,
    these exponents are the only ones possible.

    <item>If we choose <with|mode|math|p=1>, then the best constant comes to
    light by choosing <with|mode|math|u=\<b-1\><rsub|B(0,1)>>, giving the
    isoperimetric inequality.

    <item>Proof: Suppose <with|mode|math|p=1> at first. Compact
    support<with|mode|math|\<Rightarrow\>>

    <\equation*>
      u(x)\<leqslant\><big|int><rsub|-\<infty\>><rsup|\<infty\>>\|D
      u(x\<ldots\>x,y<rsub|i>,x,\<ldots\>,x)\|\<mathd\>y<rsub|i><space|1em>(i=1,\<ldots\>,n).
    </equation*>

    Then

    <\equation*>
      \|u(x)\|<rsup|n/(n-1)>\<leqslant\><left|(><big|prod><rsub|i><big|int>\<ldots\>\<mathd\>y<rsub|i><right|)><rsup|1/(n-1)>.
    </equation*>

    Integrating this gives

    <\equation*>
      <big|int>\|u\|<rsup|n/(n-1)>\<mathd\>x<rsub|1>\<leqslant\><left|(><big|int>\|D
      u\|\<mathd\>x<rsub|1><right|)><rsup|1/(n-1)><left|(><big|prod><rsub|i=2><big|int><big|int>\|D
      u\|\<mathd\>x<rsub|1>\<mathd\>y<rsub|i><right|)><rsup|1/(n-1)>
    </equation*>

    by pulling out an independent part and using generalized Hölder. Then
    iterate the same trick. To obtain for general <with|mode|math|p>, use on
    <with|mode|math|v=\|u\|<rsup|\<gamma\>>> with suitable
    <with|mode|math|\<gamma\>>.
  </itemize>

  <subsection|Poincaré and Morrey>

  <\itemize>
    <item><em|Riesz potential>: <with|mode|math|0\<less\>\<alpha\>\<less\>n>

    <\equation*>
      I<rsub|\<alpha\>>(x)=\|x\|<rsup|\<alpha\>-n>\<in\>L<rsup|1><rsub|loc>(\<bbb-R\><rsup|n>).
    </equation*>

    <item><with|mode|math|<norm|I<rsub|1>\<ast\>f|L<rsup|p>|>\<leqslant\>C<norm|f|L<rsup|p>|>.>

    <item><em|Poincaré's Inequality>: <with|mode|math|\<Omega\>> convex,
    <with|mode|math|\|\<Omega\>\|\<less\>\<infty\>>,
    <with|mode|math|d=diam(\<Omega\>)>, <with|mode|math|u\<in\>W<rsup|1,p>(\<Omega\>)>.
    Then

    <\equation*>
      <left|(><superpose|<big|int>| -|><rsub|\<Omega\>>\|u(x)-<wide|u|\<bar\>><rsub|\<Omega\>>\|<rsup|p><right|)><rsup|1/p>\<leqslant\>C*d<left|(><superpose|<big|int>|
      -|><rsub|\<Omega\>>\|D u\|<rsup|p><right|)><rsup|1/p>.
    </equation*>

    Proof: Use calculus to derive

    <\equation*>
      \|u(x)-<wide|u|\<bar\>>\|\<leqslant\><frac|d<rsup|n>|n><superpose|<big|int>|
      -|><rsub|\<Omega\>><frac|\|D u(y)\||\|x-y\|<rsup|n-1>>\<mathd\>y.
    </equation*>

    Then use potential estimate.

    <item><em|Morrey's Inequality:> <with|mode|math|u\<in\>W<rsup|1,1><rsub|loc>(\<Omega\>)>,
    <with|mode|math|0\<less\>\<alpha\>\<leqslant\>1>. If
    <with|mode|math|\<exists\>M\<gtr\>0> with

    <\equation*>
      <big|int><rsub|B(x,r)>\|D u\|\<leqslant\>M*r<rsup|n-1+\<alpha\>>,
    </equation*>

    for all <with|mode|math|B(x,r)\<subset\>\<Omega\>>. Then
    <with|mode|math|u\<in\>C<rsup|0,\<alpha\>>(\<Omega\>)> and
    <with|mode|math|osc<rsub|B(x,r)>u\<leqslant\>C*M*r<rsup|\<alpha\>>>.

    <item>Morrey=Poincaré+Campanato in <with|mode|math|W<rsup|1,1>>.

    <item><em|More general Morrey:> <with|mode|math|u\<in\>W<rsup|1,p>(\<bbb-R\><rsup|n>)>,
    <with|mode|math|n\<less\>p\<leqslant\>\<infty\>>. Then
    <with|mode|math|u\<in\>C<rsup|0,1-n/p><rsub|loc>(\<bbb-R\><rsup|n>)> and

    <\equation*>
      osc<rsub|B(x,r)>u\<leqslant\>r<rsup|1-n/p><norm|D u|L<rsup|p>|>.
    </equation*>

    If <with|mode|math|p=\<infty\>>, <with|mode|math|u> is locally Lipschitz.

    Proof: Use Jensen <with|mode|math|(\<cdot\>)<rsup|p\<cdot\><frac|1|p>>>
    on Poincaré's RHS. Then apply Campanato.
  </itemize>

  <subsection|BMO>

  <\itemize>
    <item><em|BMO seminorm>:

    <\equation*>
      [u]<rsub|BMO>\<assign\>sup<rsub|B><superpose|<big|int>|
      -><rsub|B>\|u-<wide|u|\<bar\>><rsub|B>\|\<mathd\>x
    </equation*>

    <item><with|mode|math|BMO\<assign\>{[u]<rsub|BMO>\<less\>\<infty\>}>.

    <item><em|John-Nirenberg>: <with|mode|math|W<rsup|1,n>(\<bbb-R\><rsup|n>)(*\<cap\>L<rsup|1>(\<bbb-R\><rsup|n>))\<subset\>BMO(\<bbb-R\><rsup|n>)>.

    Proof: Poincaré-then-Jensen.

    <item>For a compact domain, <with|mode|math|L<rsup|p>\<subset\>L<rsup|\<infty\>>\<subset\>BMO>.
  </itemize>

  <subsection|Imbeddings>

  <\itemize>
    <item><em|Imbedding> <with|mode|math|B<rsub|1>\<rightarrow\>B<rsub|2>>:
    <with|mode|math|\<exists\>>continuous, linear, injective map.

    <item><with|mode|math|W<rsup|1,p>(\<bbb-R\><rsup|n>)\<rightarrow\>L<rsup|p<rsup|\<ast\>>>>
    for <with|mode|math|1\<leqslant\>p\<less\>n> (Sobolev inequality)

    <item><with|mode|math|W<rsup|1,p>(\<bbb-R\><rsup|n>)\<rightarrow\>BMO>
    for <with|mode|math|p=n>

    <item><with|mode|math|W<rsup|1,p>(\<bbb-R\><rsup|n>)\<rightarrow\>C<rsup|0,1-n/p><rsub|loc>>
    (Morrey)
  </itemize>

  <with|mode|math|\<Omega\>> bounded now.

  <\itemize>
    <item><with|mode|math|W<rsup|1,p>(\<Omega\>)\<rightarrow\>L<rsup|q>(\<Omega\>)>
    for <with|mode|math|1\<less\>p\<less\>n> and
    <with|mode|math|1\<leqslant\>q\<less\>p<rsup|\<ast\>>>. Proof:
    Hölder-then-Sobolev:

    <\equation*>
      <norm|u|L<rsup|q>|>\<leqslant\><norm|u|L<rsup|p<rsup|\<ast\>>>|>\|\<Omega\>\|<rsup|1-q/p<rsup|\<ast\>>>\<leqslant\><norm|D
      u|W<rsup|1,p>|>.
    </equation*>

    <item><with|mode|math|W<rsup|1,p><rsub|0>(\<Omega\>)\<rightarrow\>C<rsup|0,1-n/p>(<wide|\<Omega\>|\<bar\>>)>
    for <with|mode|math|n\<less\>p\<leqslant\>\<infty\>>.

    <item><em|Compact imbedding> <with|mode|math|B<rsub|1>\<hookrightarrow\>B<rsub|2>>:
    The image of every bounded set in <with|mode|math|B<rsub|1>> is
    precompact in <with|mode|math|B<rsub|2>>.

    (precompact: has compact closure)

    <item>Rellich-Kondrachev:

    <\itemize>
      <item><with|mode|math|W<rsup|1,p>(\<Omega\>)\<hookrightarrow\>L<rsup|q>(\<Omega\>)>
      for <with|mode|math|1\<less\>p\<less\>n> and
      <with|mode|math|1\<leqslant\>q\<less\>p<rsup|\<ast\>>>.

      <with|color|red|In Evans, we need <with|mode|math|\<partial\>U\<in\>C<rsup|1>>.
      Our notes do not.>

      Proof:

      <\itemize>
        <item>Grab a <with|mode|math|W<rsup|1,p>>-bounded sequence
        <with|mode|math|u<rsub|m>>.

        <item>Mollify it to <with|mode|math|u<rsub|m><rsup|\<varepsilon\>>>

        <item>Use an <with|mode|math|\<varepsilon\>>-derivative trick to show
        <with|mode|math|<norm|u<rsup|\<varepsilon\>><rsub|m>-u<rsub|m>|L<rsup|1>|>\<leqslant\>\<varepsilon\><norm|D
        u<rsub|m>|L<rsup|p>|>\<rightarrow\>0>

        <item>Interpolation inequality for <with|mode|math|L<rsup|p>>:
        <with|mode|math|<norm|u<rsup|\<varepsilon\>><rsub|m>-u<rsub|m>|L<rsup|q>|>\<leqslant\><norm|u<rsup|\<varepsilon\>><rsub|m>-u<rsub|m>|L<rsup|1>|\<theta\>><norm|u<rsub|m><rsup|\<varepsilon\>>-u<rsub|m>|L<rsup|p<rsup|\<ast\>>>|>\<rightarrow\>0>,
        also using GNS.

        <item>For fixed <with|mode|math|\<varepsilon\>>,
        <with|mode|math|u<rsup|\<varepsilon\>><rsub|m>> is bounded and
        equicontinuous (directly mess with convolution).

        <item>Use Arzelà-Ascoli and a diagonal argument to finish off.
      </itemize>

      <item><with|mode|math|W<rsup|1,p><rsub|0>(\<Omega\>)\<hookrightarrow\>C<rsup|0>(<wide|\<Omega\>|\<bar\>>)\<subset\>L<rsup|p>(\<Omega\>)>
      for <with|mode|math|n\<less\>p\<leqslant\>\<infty\>>.

      Proof: Morrey's Inequality, then Arzelà-Ascoli.
    </itemize>
  </itemize>

  <section|Scalar Elliptic Equations>

  <\itemize>
    <item><with|mode|math|L u=div(A*D u+b u)+c\<cdot\>D u+d*u>.

    <item>Motivation: Calculus of Variations.

    <item><em|Weak Formulation>: <with|mode|math|u\<in\>W<rsup|1,2>(\<Omega\>)>,
    <with|mode|math|v\<in\>C<rsup|1><rsub|c>(\<Omega\>)>

    <\equation*>
      B[u,v]\<assign\><big|int><rsub|\<Omega\>>(D v<rsup|T>A*D u+b\<cdot\>D
      v*u)-(c\<cdot\>D u+d*u)v*\<mathd\>x.
    </equation*>

    <item><em|Generalized Dirichlet Problem>: <with|mode|math|L u=g+div f> on
    <with|mode|math|\<Omega\>>, <with|mode|math|u=\<varphi\>> on
    <with|mode|math|\<partial\>\<Omega\>>, i.e. <with|mode|math|B[u,v]=F(v)>
    with

    <\equation*>
      F(v)\<assign\><big|int><rsub|\<Omega\>>D v\<cdot\>f-g*v d\<mathd\>x.
    </equation*>

    <item>Assumptions:

    <\description>
      <item*|(<with|mode|math|E<rsub|1>>)><em|Strict ellipticity:>
      <with|mode|math|\<exists\>\<lambda\>\<gtr\>0:
      \<xi\><rsup|T>A\<xi\>\<geqslant\>\<lambda\>\|\<xi\>\|<rsup|2>>

      <item*|(<with|mode|math|E<rsub|2>>)><em|Boundedness:>
      <with|mode|math|A,b,c,d\<in\>L<rsup|\<infty\>>>, i.e.
      <with|mode|math|<norm|Tr(A<rsup|T>A)|L<rsup|\<infty\>>|>\<leqslant\>\<Lambda\><rsup|2>>,
      <with|mode|math|<frac|1|\<lambda\><rsup|2>>(<norm|b|\<infty\>|>+<norm|c|\<infty\>|>)+<frac|1|\<lambda\>>(<norm|d|\<infty\>|>)\<leqslant\>\<nu\>>.

      <item*|(<with|mode|math|E<rsub|3>>)><with|mode|math|div
      b+d\<leqslant\>0> weakly, i.e.

      <\equation*>
        <big|int><rsub|\<Omega\>>d*v-b\<cdot\>D v\<mathd\>x\<leqslant\>0
      </equation*>

      for <with|mode|math|v\<in\>C<rsup|1><rsub|c>(\<Omega\>)>,
      <with|mode|math|v\<geqslant\>0>.
    </description>

    <item><em|``<with|mode|math|\<leqslant\>>'' on the boundary>:
    <with|mode|math|u\<leqslant\>v\<Leftrightarrow\>(u-v)\<leqslant\>0><with|mode|math|:\<Leftrightarrow\>><with|mode|math|(u-v)<rsup|+>\<in\>W<rsup|1,2><rsub|0>(\<Omega\>)>.

    <item><em|``sup'' on the boundary>: <with|mode|math|sup<rsub|\<partial\>\<Omega\>>u=inf{k\<in\>\<bbb-R\>:u\<leqslant\>k
    <with|mode|text|on> \<partial\>\<Omega\>}>.

    <item><em|<with|mode|math|u> is a subsolution><with|mode|math|:\<Leftrightarrow\>><with|mode|math|B[u,v]\<leqslant\>F(v)><with|mode|math|\<Leftrightarrow\>><with|mode|math|L
    u\<geqslant\>g+div f>.

    <item><em|Non-divergence form:>

    <\equation*>
      0=A*D<rsup|2>U+b*\<cdot\>D u+d*u
    </equation*>

    (Not equivalent!)

    <item><em|Classical Maximum Principle:> Holds if
    <with|mode|math|d\<leqslant\>0>.

    <item><em|Weak Maximum Principle:> <with|mode|math|L
    u\<geqslant\>0\<Leftrightarrow\>B[u,v]\<leqslant\>0> for
    <with|mode|math|v\<geqslant\>0> and <with|mode|math|(E<rsub|1>)>,
    <with|mode|math|(E<rsub|2>)>, <with|mode|math|(E<rsub|3>)>. Then
    <with|mode|math|sup<rsub|\<Omega\>>u\<leqslant\>sup<rsub|\<partial\>\<Omega\>>u<rsup|+>>.

    Proof:\ 

    <\itemize>
      <item>Use <with|mode|math|B[u,v]\<leqslant\>0> for
      <with|mode|math|v\<geqslant\>0> and <with|mode|math|(E<rsub|3>)> to
      establish

      <\equation*>
        <big|int>D v<rsup|T>A*D v-(b+c)D u\<cdot\>v\<leqslant\><big|int>d(u*v)-b\<cdot\>D(u*v)\<leqslant\>0.
      </equation*>

      Note that <with|mode|math|u*v> is the new test function in
      <with|mode|math|(E<rsub|3>)>. Consequently

      <\equation*>
        <big|int>D v<rsup|T>A*D v\<leqslant\><big|int>(b+c)D u\<cdot\>v.
      </equation*>

      <item>Suppose <with|mode|math|l=sup<rsub|\<partial\>\<Omega\>>u\<leqslant\>k\<less\>sup<rsub|\<Omega\>>u>.
      Set <with|mode|math|\<Gamma\>\<assign\>{u\<gtr\>k}> and achieve a
      <with|mode|math|<norm|D v|L<rsup|2>|>\<leqslant\>C<norm|v|L<rsup|2>|>>
      estimate by using ellipticity, the above and boundedness. Use the
      Sobolev inequality to get <with|mode|math|<norm|v|L<rsup|2<rsup|\<ast\>>>|>\<leqslant\>\<cdots\>\<leqslant\>\|\<Gamma\>\|<rsup|1/n><norm|v|L<rsup|2<rsup|\<ast\>>>|>>,
      and so <with|mode|math|\|\<Gamma\>\|\<gtr\>0> independently of
      <with|mode|math|k>. Let <with|mode|math|k\<rightarrow\>sup<rsub|\<Omega\>>>
      to obtain a contradiction. (Note <with|mode|math|sup<rsub|\<Omega\>>\<less\>\<infty\>>
      because <with|mode|math|u\<in\>W<rsup|1,2>(\<Omega\>)>.)
    </itemize>

    Remarks:

    <\itemize>
      <item>Implies uniqueness.

      <item>No assumptions on boundedness, smoothness or connectedness of
      <with|mode|math|\<Omega\>>.
    </itemize>

    <item>Implies uniqueness.
  </itemize>

  <subsection|Existence Theory>

  <\itemize>
    <item><em|Existence:> <with|mode|math|\<Omega\>> bounded,
    <with|mode|math|(E<rsub|1>)>, <with|mode|math|(E<rsub|2>)>,
    <with|mode|math|(E<rsub|3>)>. Then <with|mode|math|\<exists\>!> solution
    of the generalized Dirichlet problem.

    <\itemize>
      <item>Reduce BC to <with|mode|math|H<rsup|1,2><rsub|0>> by subtracting
      arbitrary function and handling RHS.

      <item>Prove coercivity estimate

      <\equation*>
        B[u,u]\<geqslant\><frac|\<lambda\>|2><big|int><rsub|\<Omega\>>\|D
        u\|<rsup|2>\<mathd\>x-\<lambda\>\<nu\><rsup|2><big|int><rsub|\<Omega\>>\|u\|<rsup|2>\<mathd\>x.
      </equation*>

      (Uses: <with|mode|math|(E<rsub|1>)>, <with|mode|math|(E<rsub|2>)>,
      <with|mode|math|2*a*b\<leqslant\>\<lambda\>a<rsup|2>+b<rsup|2>/\<lambda\>>.

      <with|color|red|(In Evans, Poincaré enters here. How?)>

      (For <with|mode|math|\<Delta\>>, Poincaré suffices to show coercivity.)

      <item><with|mode|math|Id:W<rsup|1,2><rsub|0>\<rightarrow\>(W<rsup|1,2><rsub|0>)<rsup|\<ast\>>>
      is compact.\ 

      <\equation*>
        Id=*<wide*|(L<rsup|2>\<rightarrow\>\<cal-H\><rsup|\<ast\>>)|\<wide-underbrace\>><rsub|continous>\<circ\><wide*|(\<cal-H\>\<rightarrow\>L<rsup|2>)|\<wide-underbrace\>><rsub|<with|mode|text|compact>>.
      </equation*>

      <item><with|mode|math|L<rsub|\<sigma\>>\<assign\>L-\<sigma\>Id>.
      (<with|mode|math|L\<approx\>\<Delta\>> has negative eigenvalues
      already. But they might be pushed upward by the first- and zeroth-order
      junk. So we might have to make them even more negative to succeed.)

      <item><with|mode|math|\<rightarrow\>><with|mode|math|B<rsub|\<sigma\>>[u,v]=B[u,v]+\<sigma\><ip|u|v|L<rsup|2>|>>,
      coercivity is maintained.

      <item>Lax-Milgram shows existence of inverse
      <with|mode|math|L<rsub|\<sigma\>><rsup|-1>> for the not-so-bad operator
      <with|mode|math|L<rsub|\<sigma\>>>.

      <item>Start with <with|mode|math|L u=g+div f>, introduce
      <with|mode|math|L<rsub|\<sigma\>>>, multiply by
      <with|mode|math|L<rsub|\<sigma\>><rsup|-1>> and see what happens.

      <item>Weak maximum principle provides uniqueness for
      <with|mode|math|L>, so that the Fredholm alternative provides existence
      when combined with Rellich.
    </itemize>
  </itemize>

  <subsection|Regularity>

  <\itemize>
    <item>Assumptions:

    <\itemize>
      <item><with|mode|math|(R<rsub|1>)>: <with|mode|math|(E<rsub|1>)>,
      <with|mode|math|(E<rsub|2>)>.

      <item><with|mode|math|(R<rsub|2>)>:
      <with|mode|math|f\<in\>L<rsup|q>(\<Omega\>)>,
      <with|mode|math|g\<in\>L<rsup|q/2>>, <with|mode|math|q\<gtr\>n>.
    </itemize>

    <item><with|mode|math|(R<rsub|1>)>, <with|mode|math|L*u=g>.
    \ <with|mode|math|A>, <with|mode|math|b> Lipschitz. Then for
    <with|mode|math|\<Omega\><rprime|'>\<subset\>\<subset\>\<Omega\>> we have

    <\equation*>
      <norm|u|W<rsup|2,2>(\<Omega\><rprime|'>)|>\<leqslant\>C<left|(><norm|u|W<rsup|1,2>(\<Omega\>)|>+<norm|g|L<rsup|2>(\<Omega\>)|><right|)>.
    </equation*>

    Proof:

    <\itemize>
      <item>Finite Differences.
    </itemize>
  </itemize>

  <subsection|Harnack Inequality Stuff>

  <\itemize>
    <item><em|(Ladyzhenskaya/Uraltseva)>: <with|mode|math|(R<rsub|1>)>,
    <with|mode|math|(R<rsub|2>)>. <with|mode|math|u\<in\>W<rsup|1,2>> a
    subsolution, <with|mode|math|u\<leqslant\>0> on
    <with|mode|math|\<partial\>\<Omega\>>. Then:

    <\equation*>
      sup<rsub|\<Omega\>>u\<leqslant\>C<left|(><norm|u<rsup|+>|L<rsup|2>(\<Omega\>)|>+k<right|)>,
    </equation*>

    where

    <\equation*>
      k=<frac|1|\<lambda\>><left|(><norm|f|L<rsup|q>|>+<norm|g|L<rsup|q/2>|><right|)>.
    </equation*>

    Proof:

    <\itemize>
      <item>
    </itemize>

    <item><em|Local Boundedness>: <with|mode|math|(R<rsub|1>)>,
    <with|mode|math|(R<rsub|2>)>. <with|mode|math|u\<in\>W<rsup|1,2>> a
    subsolution. Then:

    <\equation*>
      sup<rsub|B(y,R)>u\<leqslant\>C<left|(>R<rsup|-n/p><norm|u<rsup|+>|L<rsup|2>(\<Omega\>)|>+k(R)<right|)>,
    </equation*>

    where

    <\equation*>
      k(R)=<frac|R<rsup|1-n/q>|\<lambda\>><left|(><norm|f|L<rsup|q>|>+R<rsup|1-n/q><norm|g|L<rsup|q/2>|><right|)>.
    </equation*>

    <item><em|Weak Harnack Inequality>: <with|mode|math|(R<rsub|1>)>,
    <with|mode|math|(R<rsub|2>)>, <with|mode|math|u\<in\>W<rsup|1,2>(\<Omega\>)>
    a supersolution and <with|mode|math|u\<geqslant\>0> in
    <with|mode|math|B(y,4R)\<subset\>\<Omega\>>. Then

    <\equation*>
      R<rsup|-n/p><norm|u|L<rsup|p>(B(2R))|>\<leqslant\>C<left|(>inf<rsub|y\<in\>B(y,R)>u+k(R)<right|)>.
    </equation*>

    <item><em|Strong Harnack Inequality:> <with|mode|math|(R<rsub|1>)>,
    <with|mode|math|(R<rsub|2>)>, <with|mode|math|u\<in\>W<rsup|1,2>(\<Omega\>)>
    a solution with <with|mode|math|u\<geqslant\>0>. Then

    <\equation*>
      sup<rsub|B(y,R)>u\<leqslant\>C<left|(>inf<rsub|B(y,R)>u+k(R)<right|)>.
    </equation*>

    <item><em|Strong Maximum Principle:> <with|mode|math|(R<rsub|1>)>,
    <with|mode|math|(R<rsub|2>)>, <with|mode|math|(E<rsub|3>)>,
    <with|mode|math|\<Omega\>> connected, <with|mode|math|u\<in\>W<rsup|1,2>>
    a subsolution <with|mode|math|L u\<geqslant\>0>. If

    <\equation*>
      sup<rsub|B>u=sup<rsub|\<Omega\>>u,
    </equation*>

    then <with|mode|math|u=const>.

    Proof: Weak Harnack shows <with|mode|math|{u=sup<rsub|\<Omega\>>u}> is
    open. <with|mode|math|{u=sup<rsub|\<Omega\>>u}> is relatively closed in
    <with|mode|math|\<Omega\>>. Therefore
    <with|mode|math|{u=sup<rsub|\<Omega\>>u}=\<Omega\>>.

    <with|color|red|Why is <with|mode|math|L const=0>?>

    <with|color|red|How dow we know the ``relatively closed'' part?>

    <item><em|DeGiorgi/Nash>: <with|mode|math|(R<rsub|1>)>,
    <with|mode|math|(R<rsub|2>)>, <with|mode|math|u\<in\>W<rsup|1,2>>
    solution of <with|mode|math|L u=g+div f>. Then <with|mode|math|f> is
    locally Hölder and

    <\equation*>
      osc<rsub|B(y,R)>u\<leqslant\>C*R<rsup|\<alpha\>><left|(>R<rsub|0><rsup|-\<alpha\>>sup<rsub|B(y,R<rsub|0>)>\|u\|+k<right|)>
    </equation*>

    if <with|mode|math|0\<less\>R\<leqslant\>R<rsub|0>>.

    Proof:

    <\itemize>
      <item>
    </itemize>
  </itemize>

  <section|Calculus of Variations>

  <with|mode|math|\<Omega\>> open, bounded.

  <\itemize>
    <item>Idea: solution <with|mode|math|u>, smooth variation
    <with|mode|math|\<varphi\>>, functional <with|mode|math|I>.
    <with|mode|math|\<partial\><rsub|\<varepsilon\>>I(*u+\<varepsilon\>\<varphi\>)\|<rsub|\<varepsilon\>=0>=0>.
    Integrate by parts, <with|mode|math|\<varphi\>> was
    arbitrary<with|mode|math|\<rightarrow\>>PDE.

    <item><with|mode|math|u:\<Omega\>\<rightarrow\>\<bbb-R\><rsup|m>>
    deformation, <with|mode|math|D u:\<Omega\>\<rightarrow\>\<bbb-R\><rsup|m\<times\>n>>,
    <with|mode|math|F:\<bbb-R\><rsup|m\<times\>n>[\<times\>\<bbb-R\><rsup|n>]\<rightarrow\>\<bbb-R\>>.

    <\equation*>
      I[u]=<big|int><rsub|\<Omega\>>F(D u(x),u)\<mathd\>x.
    </equation*>

    Looking for <with|mode|math|inf<rsub|u\<in\>\<cal-A\>>I[u]>, where
    <with|mode|math|\<cal-A\>=W<rsup|1,2><rsub|0>(\<Omega\>)>.

    <item><em|Example: Dirichlet's Principle:> <with|mode|math|\<Omega\>>
    open, bounded

    <\equation*>
      I[u]=<big|int><rsub|\<Omega\>><left|(><frac|1|2>\|D u\|<rsup|2>-g
      u<right|)>\<mathd\>x.
    </equation*>

    <\itemize>
      <item>Bounded below: <with|mode|math|\<varepsilon\>a<rsup|2>+b<rsup|2>/\<varepsilon\>>,
      Sobolev (<with|mode|math|2<rsup|\<ast\>>\<gtr\>2>), Hölder as
      <with|mode|math|<norm|u|L<rsup|2>|>\<leqslant\><norm|u|L<rsup|2<rsup|\<ast\>>>|>\|\<Omega\>\|<rsup|1/n>>,
      gives

      <\equation*>
        I[u]\<geqslant\>c<norm|u|W<rsup|1,2><rsub|0>|2>-<frac|1|2\<varepsilon\>><norm|g|L<rsup|2>|2>.
      </equation*>

      <item>Bounded above by <with|mode|math|<norm|u|W<rsup|1,2><rsub|0>|2>+<norm|g|L<rsup|2>|>>.

      <item><with|mode|math|I> wlsc because <with|mode|math|F> convex.

      <item>strictly convex (unproven)<with|mode|math|\<Rightarrow\>>uniqueness.
    </itemize>

    <item><em|Weak lower semicontinuity>:
    <with|mode|math|u<rsub|k>\<rightharpoonup\>u\<Rightarrow\>I[u]\<leqslant\>liminf<rsub|k\<rightarrow\>\<infty\>>I[u<rsub|k>]>.

    <item><with|mode|math|F> convex<with|mode|math|\<Rightarrow\>><with|mode|math|I>
    wlsc in <with|mode|math|W<rsup|1,p><rsub|0>(\<Omega\>)>.

    Proof: Use representation of convex <with|mode|math|F> as limit of
    increasing sequence <with|mode|math|{F<rsub|N>}> of piecewise affine
    functions. Implies <with|mode|math|<big|int>F<rsub|N>(D
    u<rsub|k>)<above|\<rightarrow\>|k><big|int>F<rsub|N>(D u)> (weak
    convergence <with|mode|math|\<heartsuit\>> linear/affine functions). Then

    <\equation*>
      <big|int>F<rsub|N>(D u)<above|\<leqslant\>|<with|mode|text|<with|mode|math|F<rsub|N>
      <with|mode|text|incr.>>>>liminf<rsub|k\<rightarrow\>\<infty\>><big|int>F(D
      u<rsub|k>)=liminf<rsub|k\<rightarrow\>\<infty\>>I[u<rsub|k>]
    </equation*>

    and MCT.

    <item><em|Jensen>: <with|mode|math|F(<with|mode|text|w-<with|mode|math|\<ast\>>>lim
    g<rsub|k>)\<leqslant\><with|mode|text|w-<with|mode|math|\<ast\>>>lim
    F(g<rsub|k>)>.

    <item><em|Euler-Lagrange Equation>: Weak form obtained from
    <with|mode|math|i(\<tau\>)=I[u+\<tau\>v]>, where <with|mode|math|u=argmin
    I[u]> and looking at <with|mode|math|i<rprime|'>(0)=0.>

    <\equation*>
      -div(F<rsub|p>(D u))+F<rsub|u>(D u,u)=0.
    </equation*>

    Also <with|mode|math|i<rprime|''>(0)\<geqslant\>0>.

    <item><em|Motivation for Convexity>: <with|mode|math|\<rho\>(s)=<with|mode|text|0-1
    sawtooth>>. <with|mode|math|\<rho\><rprime|'>=1> a.e..
    <with|mode|math|v<rsub|\<varepsilon\>>(x)=\<varepsilon\>\<zeta\>(x)\<rho\>(x\<cdot\>\<xi\>/\<varepsilon\>)>.

    <\equation*>
      <frac|\<partial\>v<rsub|\<varepsilon\>>|\<partial\>x<rsub|i>>(x)\<approx\>\<zeta\>(x)\<rho\><rprime|'>(x\<cdot\>\<xi\>/\<varepsilon\>)\<approx\>\<zeta\>(x)\<xi\>.
    </equation*>

    Consider <with|mode|math|i<rprime|''>(0)\<geqslant\>0><with|mode|math|\<Rightarrow\>><with|mode|math|\<xi\><rsup|T>D<rsup|2>F\<xi\>\<geqslant\>0>
    pops out.

    <item><with|mode|math|m=1><with|mode|math|\<Rightarrow\>>(wlsc<with|mode|math|\<Leftrightarrow\>>convexity).

    Proof: ``<with|mode|math|\<Leftarrow\>>'': shown above.
    ``<with|mode|math|\<Rightarrow\>>'': <with|mode|math|2<rsup|n*k>> cube
    grid on <with|mode|math|[0,1]<rsup|n>>,
    <with|mode|math|v\<in\>C<rsup|\<infty\>><rsub|c>>.

    <\eqnarray*>
      <tformat|<table|<row|<cell|u<rsub|k>(x)>|<cell|=>|<cell|<frac|1|2<rsup|k>>v(2<rsup|k>(x-cell
      center))+z\<cdot\>x.>>|<row|<cell|D u<rsub|k>(x)>|<cell|=>|<cell|D
      v(2<rsup|k>(x-cell center))+z.>>>>
    </eqnarray*>

    <with|mode|math|u<rsub|k>\<rightarrow\>z\<cdot\>x>, <with|mode|math|D
    u<rsub|k>\<rightharpoonup\>D u>. Then

    <\equation*>
      F(z)<above|\<leqslant\>|wlsc>liminf<rsub|k\<rightarrow\>\<infty\>><big|sum><rsub|l><big|int><rsub|Q<rsub|l>>F(D
      u<rsub|k>)=<big|int><rsub|[0,1]<rsup|n>>F(z+D v)
    </equation*>

    Thus <with|mode|math|I[u]> has a minimum at the straight line, and for
    <with|mode|math|i(\<tau\>)=I[u+\<tau\>v]>,
    <with|mode|math|i<rprime|'>(0)=0>, <with|mode|math|i<rprime|''>(0)\<geqslant\>0>,
    convexity follows as above.
  </itemize>

  <subsection|Quasiconvexity>

  <with|mode|math|m\<geqslant\>2>, <with|mode|math|\<cal-A\>=W<rsup|1,p>\<cap\>{u=g}<rsub|\<partial\>\<Omega\>>>.
  <with|mode|math|1\<less\>p\<less\>\<infty\>>.
  <with|mode|math|F\<in\>C<rsup|\<infty\>>>,
  <with|mode|math|F(A)\<geqslant\>c<rsub|1>\|A\|<rsup|p>-c<rsub|2>>.

  <\equation*>
    I[u]=<big|int><rsub|\<Omega\>>F(D u(x))\<mathd\>x
  </equation*>

  <\itemize>
    <item>Sawtooth calculation yields <em|rank-one convexity>

    <\equation*>
      (\<eta\>\<otimes\>\<xi\>)<rsup|T>D<rsup|2>F(P)(\<eta\>\<otimes\>\<xi\>)
    </equation*>

    <with|mode|math|\<Leftrightarrow\>><with|mode|math|F(P+t(\<eta\>\<otimes\>\<xi\>))>
    convex in <with|mode|math|t>.

    <item><em|Quasiconvexity>: <with|mode|math|F>
    <em|quasiconvex><with|mode|math|:\<Leftrightarrow\>><with|mode|math|\<forall\>A\<in\>\<bbb-R\><rsup|m\<times\>n>>,
    <with|mode|math|v\<in\>C<rsup|\<infty\>><rsub|c>([0,1]<rsup|n>,\<bbb-R\><rsup|m>)>:

    <\equation*>
      F(A)\<leqslant\><big|int><rsub|[0,1]<rsup|n>>F(A+D v)
    </equation*>

    <item>If <with|mode|math|\|F(A)\|\<leqslant\>C(1+\|A\|<rsup|p>)>, then
    <with|mode|math|F> QC<with|mode|math|\<Leftrightarrow\>><with|mode|math|I>
    wlsc.

    <\itemize>
      <item>``<with|mode|math|\<Rightarrow\>>'': Subdivide domain into cubes,

      <\equation*>
        <big|int><rsub|\<Omega\>>F(D u)\<approx\><big|int><rsub|\<Omega\>>F(<with|mode|text|affine
        approx to <with|mode|math|D u>>)<above|\<leqslant\>|QC><big|int><rsub|\<Omega\>>F
        (D u<rsub|k>)+errors.
      </equation*>

      Use measure theory to keep concentrations (Dirac bumps?) of
      <with|mode|math|D u> or <with|mode|math|D u<rsub|k>> away from cube
      boundaries. Mop up the error terms.

      <item>``<with|mode|math|\<Leftarrow\>>'': cubes calculation above.
    </itemize>

    <item><em|Polyconvex>: <with|mode|math|F> is a convex function of minors
    of <with|mode|math|A>.

    <item>Convex<with|mode|math|\<Rightarrow\>>PC<with|mode|math|\<Rightarrow\>>QC<with|mode|math|\<Rightarrow\>>R1C
    (converse false).

    Proof of PC<with|mode|math|\<Rightarrow\>>QC:
    PC<with|mode|math|\<Rightarrow\>>wlsc (use
    convex<with|mode|math|\<Rightarrow\>>wlsc argument for each minor).
    wlsc<with|mode|math|\<Rightarrow\>>QC.

    <item><with|mode|math|\|D F(A)\|\<leqslant\>C(1+\|A\|<rsup|p-1>)>.

    Proof: Exploit growth estimate above, and
    QC<with|mode|math|\<Rightarrow\>>R1C. Use
    <with|mode|math|f(t)=F(A+t(\<eta\>\<otimes\>\<xi\>))>, which is
    convex<with|mode|math|\<Rightarrow\>>locally
    Lip<with|mode|math|\<Rightarrow\>>locally
    <with|mode|math|\|f<rprime|'>(0)\|\<leqslant\>max\|f\|>.
  </itemize>

  <subsection|Null Lagrangians, Determinants>

  <\itemize>
    <item><with|mode|math|F(D u)> is a <em|null Lagrangian> if E-L

    <\equation*>
      div(D F(D u))=\<partial\>x<rsub|j><left|(>\<partial\><rsub|A<rsub|i,j>>F(D
      u))=0
    </equation*>

    holds for every <with|mode|math|u\<in\>C<rsup|2>>.

    <item><with|mode|math|F> null Langrangian. Then

    <\equation*>
      u=<wide|u|~> <with|mode|text|on> \<partial\>\<Omega\><space|1em>\<Rightarrow\><space|1em>I[u]=I[<wide|u|~>].
    </equation*>

    Proof: <with|mode|math|i(\<tau\>)\<assign\>I[\<tau\>u+(1-\<tau\>)u]>.
    <with|mode|math|i<rprime|'>(\<tau\>)=0> by E-L.

    <item><em|Cofactor matrix:>

    <\itemize>
      <item><with|mode|math|cof(A<rsub|i,j>)=det(A<rsub|\<setminus\>i,\<setminus\>j>)>.

      <item><with|mode|math|A<rsup|-1>=<frac|1|det A>*(cof A)<rsup|T>>.

      <item><with|mode|math|\<Rightarrow\>A<rsup|T>cof A=det A\<cdot\>Id>

      <item><with|mode|math|\<Rightarrow\>\<partial\><rsub|A<rsub|i,j>>det(A)=(cof
      A)<rsub|i,j>>
    </itemize>

    <item><with|mode|math|det(D u)> is a null Lagrangian, i.e.

    <\itemize>
      <item><with|mode|math|div(D det(D u))>=<with|mode|math|div(cof(D u))=0>

      <item>Plug and chug if <with|mode|math|det(D u)\<neq\>0>, otherwise add
      <with|mode|math|\<varepsilon\>Id>.
    </itemize>

    <item><with|mode|math|u<rsub|k>\<rightharpoonup\>u> in
    <with|mode|math|W<rsup|1,p>>, <with|mode|math|n\<less\>p\<less\>\<infty\>>.
    <with|mode|math|\<Rightarrow\>><with|mode|math|det(D
    u<rsub|k>)\<rightharpoonup\>det(D u)> in <with|mode|math|L<rsup|p/n>>.
    (Morrey/Reshetnyak)

    <\itemize>
      <item>Reduce dimension of problem by one by reducing to ``does the
      cofactor matrix converge''?

      <item>Use <with|mode|math|det(D u)=div<left|(><frac|1|n>cof(D
      u)<rsup|T>u<right|)>>.

      <item>Morrey (<with|mode|math|n\<less\>p>!) implies uniform boundedness
      in <with|mode|math|C<rsup|0,1-n/p>>, then use A-A to extract uniformly
      converging subsequence, settling the deal for the leftover
      <with|mode|math|u> besides the cofactor matrix.

      <item>(also holds for <with|mode|math|p=n> if <with|mode|math|det(D
      u<rsub|k>)\<geqslant\>0>--no proof.)
    </itemize>

    <item><em|No Retract Theorem>: <with|mode|math|B=B(0,1)>. There is no
    continuous map <with|mode|math|u:<wide|B|\<bar\>>\<rightarrow\>\<partial\><wide|B|\<bar\>>>
    with <with|mode|math|u(x)=x> on <with|mode|math|\<partial\>B>.

    Proof: Suppose there is a retract <with|mode|math|w>. By comparison with
    <with|mode|math|Id> and identity on the boundary,

    <\equation*>
      <big|int>det(D w)=\|B\|.
    </equation*>

    OTOH, <with|mode|math|\|w\|<rsup|2>=1\<Rightarrow\>(D
    w)<rsup|T>w=0\<Rightarrow\>det(D w)=0>. Lose smoothness requirement by
    continuously extending by Id, mollifying and using
    <with|mode|math|B(0,2)> then.

    <item><em|Brouwer's Fixed Point Theorem>:
    <with|mode|math|u:<wide|B|\<bar\>>\<rightarrow\><wide|B|\<bar\>>>
    continuous. <with|mode|math|\<exists\>x\<in\><wide|B|\<bar\>>:u(x)=x>.

    Proof: Assume no fixed point. <with|mode|math|w:B\<rightarrow\>\<partial\>B>
    is the point on <with|mode|math|\<partial\>B> hit by the ray from
    <with|mode|math|u(x)> to <with|mode|math|x>.

    <with|mode|math|w> is a retract because <with|mode|math|w> hits
    <with|mode|math|\<partial\>B> in <with|mode|math|x> if
    <with|mode|math|x\<in\>\<partial\>B>. <with|mode|math|w> is continuous.

    <item><em|Degree of a map:> <with|mode|math|u\<in\>W<rsup|1,1>>

    <\equation*>
      deg(u)=<superpose|<big|int>| -><rsub|B>det(D u).
    </equation*>

    Definable for continuous functions by approximation. Is an integer.
  </itemize>

  <section|Navier-Stokes Equations>

  <with|mode|math|G> open, <with|mode|math|<wide|G|^>\<assign\>G\<times\>(0,\<infty\>)>
  space-time.

  <\itemize>
    <item><em|Navier-Stokes Equations>:

    <\equation*>
      D u/D t=(\<nu\>\<Delta\>u-\<nabla\>p)+f,
    </equation*>

    <with|mode|math|(\<ast\>)> is the material derivative <with|mode|math|D
    \<ast\>/D t=\<partial\><rsub|t>\<ast\>+u\<cdot\>\<nabla\>\<ast\>>.

    <item><with|mode|math|\<nu\>=0><with|mode|math|\<Rightarrow\>>Euler
    equation. <with|mode|math|\<nu\>\<neq\>0><with|mode|math|\<Rightarrow\>>may
    as well assume <with|mode|math|\<nu\>=1>.

    <item><em|Conservation of mass>: <with|mode|math|\<partial\><rsub|t>\<rho\>+div(\<rho\>u)=0>.
    Assume <with|mode|math|D\<rho\>/D t=0><with|mode|math|\<Rightarrow\>><with|mode|math|\<nabla\>\<cdot\>u=0>.

    <item><em|Pressures in a smooth incompressible flow are superharmonic>:
    Take div of NSE.

    <item><em|Steady flows>: <with|mode|math|u\<cdot\>\<nabla\>u+\<nabla\>p=\<nu\>\<Delta\>u>.

    <item><em|Bernoulli's Theorem>: ideal (<with|mode|math|\<nu\>=0>), steady
    flow <with|mode|math|u\<cdot\>\<nabla\>u+\<nabla\>p=0><with|mode|math|\<Rightarrow\>><with|mode|math|\<nabla\>(u<rsup|2>/2+p)=0>

    <with|mode|math|\<Rightarrow\>><with|mode|math|u<rsup|2>/2+p=const>
    (still need conservation of mass <with|mode|math|\<nabla\>\<cdot\>u=0>)

    <item><em|Vorticity>: <with|mode|math|\<omega\>=curl u>

    <\eqnarray*>
      <tformat|<table|<row|<cell|\<partial\><rsub|t>\<omega\>+\<nabla\>\<times\>(u\<cdot\>\<nabla\>u)>|<cell|=>|<cell|\<Delta\>\<omega\>,>>|<row|<cell|\<nabla\>\<cdot\>u>|<cell|=>|<cell|0,>>|<row|<cell|\<nabla\>\<times\>u>|<cell|=>|<cell|\<omega\>.>>>>
    </eqnarray*>

    In 2D, <with|mode|math|\<nabla\>\<times\>(u*\<cdot\>\<nabla\>u)> becomes
    <with|mode|math|u\<cdot\>\<nabla\>u>.

    <item><em|Helmholtz Projection>: <with|mode|math|P=L<rsup|2><with|mode|text|->closure
    {\<nabla\>\<varphi\>:\<varphi\>\<in\>C<rsup|\<infty\>><rsub|c>(\<bbb-R\><rsup|n>)}>.
    <with|mode|math|P<rsup|\<perp\>>> (note <with|mode|math|P> closed!) is
    <em|divergence-free>.

    <with|mode|math|L<rsup|2>=P\<oplus\>P<rsup|\<perp\>>>

    Example: Divergence-free field from sem. 1 final: (continuous)
    boundary-normal field matters, (discontinuous) tangential field does not.

    <item><em|Weak formulation>:

    <\itemize>
      <item>Take <with|mode|math|a\<in\>C<rsup|\<infty\>><rsub|c>(<wide|G|^>,\<bbb-R\><rsup|n>)>
      div-free, dot NSE with it,

      <item>i. by parts second term, popping the derivative onto
      <with|mode|math|a*u>-product, pull apart, one term is zero,

      <item><with|mode|math|<big|int>a\<cdot\>\<nabla\>p=-<big|int>(div
      a)p=0>
    </itemize>

    gives

    <\eqnarray*>
      <tformat|<table|<row|<cell|(W1)<space|1em>-<big|int><rsub|<wide|G|^>>\<partial\><rsub|t>a*\<cdot\>u+\<nabla\>a\<cdot\>(u\<otimes\>u)+\<Delta\>a\<cdot\>u\<mathd\>x\<mathd\>t>|<cell|=>|<cell|0>>|<row|<cell|(W2)<space|1em><big|int><rsub|<wide|G|^>>\<nabla\>\<varphi\>\<cdot\>u>|<cell|=>|<cell|0<space|1em>(\<varphi\>\<in\>C<rsup|\<infty\>><rsub|c>(\<bbb-R\><rsup|n>))>>>>
    </eqnarray*>

    (where <with|mode|math|A*\<cdot\>B=tr(A<rsup|T>B)>)

    <item><with|mode|math|V\<assign\><norm|\<cdot\>|V|><with|mode|text|->closure{a\<in\>C<rsup|\<infty\>><rsub|c>(<wide|G|^>,\<bbb-R\><rsup|n>),\<nabla\>\<cdot\>a=0}>

    <\equation*>
      <norm|a|V|>\<assign\><big|int><rsub|<wide|G|^>>\|a\|<rsup|2>+\|\<nabla\>a\|<rsup|2>\<mathd\>x\<mathd\>t.
    </equation*>

    <item>Space for ICs: <with|mode|math|P<rsub|0>\<assign\>P\<cap\>L<rsup|2><with|mode|text|->closure
    {C<rsup|\<infty\>><rsub|c>}> to replicate <with|mode|math|u=0> on
    <with|mode|math|\<partial\>G>.

    <item><em|Existence, Energy Inequality:>
    <with|mode|math|u<rsub|0>\<in\>P<rsub|0><rsup|\<perp\>>>.
    <with|mode|math|\<exists\>u\<in\>V>:\ 

    <\itemize>
      <item>(W1), (W2)

      <item><em|continuous docking to IC>:
      <with|mode|math|<norm|u(t,\<cdot\>)-u<rsub|0>|L<rsup|2>(G)|>\<rightarrow\>0>
      as <with|mode|math|t\<rightarrow\>0>,

      <item><em|energy equality>:\ 

      <\equation*>
        <frac|\<mathd\>|\<mathd\>t><norm|u|L<rsup|2>|>=2<norm|\<nabla\>u|L<rsup|2>|>.
      </equation*>

      Equivalently for <with|mode|math|t\<gtr\>0>,

      <\equation*>
        <big|int><rsub|G>\|u(x,t)\|<rsup|2>+<big|int><rsub|0><rsup|t><big|int><rsub|G>\|\<nabla\>u(x,s)\|<rsup|2>\<mathd\>x\<mathd\>s\<leqslant\><frac|1|2><big|int><rsub|G>\|u<rsub|0>(x)\|<rsup|2>\<mathd\>x.
      </equation*>
    </itemize>
  </itemize>
</body>

<\initial>
  <\collection>
    <associate|language|american>
    <associate|page-type|letter>
    <associate|par-first|0>
    <associate|preamble|false>
  </collection>
</initial>

<\references>
  <\collection>
    <associate|auto-1|<tuple|1|1>>
    <associate|auto-10|<tuple|4.2|8>>
    <associate|auto-11|<tuple|5|8>>
    <associate|auto-12|<tuple|6|9>>
    <associate|auto-13|<tuple|6.1|10>>
    <associate|auto-14|<tuple|7|11>>
    <associate|auto-15|<tuple|8|12>>
    <associate|auto-16|<tuple|9|14>>
    <associate|auto-17|<tuple|10|15>>
    <associate|auto-18|<tuple|10.1|16>>
    <associate|auto-19|<tuple|10.2|16>>
    <associate|auto-2|<tuple|2|2>>
    <associate|auto-20|<tuple|10.3|17>>
    <associate|auto-21|<tuple|10.4|17>>
    <associate|auto-22|<tuple|10.5|18>>
    <associate|auto-23|<tuple|11|18>>
    <associate|auto-24|<tuple|11.1|19>>
    <associate|auto-25|<tuple|11.2|20>>
    <associate|auto-26|<tuple|11.3|20>>
    <associate|auto-27|<tuple|12|21>>
    <associate|auto-28|<tuple|12.1|22>>
    <associate|auto-29|<tuple|12.2|22>>
    <associate|auto-3|<tuple|3|2>>
    <associate|auto-30|<tuple|13|23>>
    <associate|auto-31|<tuple|14|?>>
    <associate|auto-4|<tuple|3.1|4>>
    <associate|auto-5|<tuple|3.2|4>>
    <associate|auto-6|<tuple|3.3|5>>
    <associate|auto-7|<tuple|3.4|5>>
    <associate|auto-8|<tuple|4|6>>
    <associate|auto-9|<tuple|4.1|7>>
    <associate|eq:elliptic-equation-hstar|<tuple|2|?>>
    <associate|eq:laplace-solrep|<tuple|1|3>>
    <associate|eq:maxprinciple-proof|<tuple|2|?>>
    <associate|eq:moser-main-inequality|<tuple|3|22>>
    <associate|eq:moser-test-deriv|<tuple|2|22>>
    <associate|eq:reg-rewritten-fd|<tuple|3|?>>
    <associate|eq:reg-rewritten-weak|<tuple|2|?>>
    <associate|lem:coercivity|<tuple|3|?>>
  </collection>
</references>

<\auxiliary>
  <\collection>
    <\associate|toc>
      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|1<space|2spc>General
      Stuff> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-1><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|2<space|2spc>Equations>
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-2><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|3<space|2spc>Laplace's
      Equation> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-3><vspace|0.5fn>

      <with|par-left|<quote|1.5fn>|3.1<space|2spc>Energy Methods
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-4>>

      <with|par-left|<quote|1.5fn>|3.2<space|2spc>Potentials
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-5>>

      <with|par-left|<quote|1.5fn>|3.3<space|2spc>Lebesgue's Thorn
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-6>>

      <with|par-left|<quote|1.5fn>|3.4<space|2spc>Capacity
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-7>>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|4<space|2spc>Heat
      Equation> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-8><vspace|0.5fn>

      <with|par-left|<quote|1.5fn>|4.1<space|2spc>Difference Schemes and
      Probabilistic Interpretation <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-9>>

      <with|par-left|<quote|1.5fn>|4.2<space|2spc>Hearing the shape of a drum
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-10>>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|5<space|2spc>Wave
      equation> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-11><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|6<space|2spc>Distributions/Fourier
      Transform> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-12><vspace|0.5fn>

      <with|par-left|<quote|1.5fn>|6.1<space|2spc>Tempered Distributions
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-13>>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|7<space|2spc>Hyperbolic
      Equations> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-14><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|8<space|2spc>Conservation
      Laws> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-15><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|9<space|2spc>Hamilton-Jacobi
      Equations> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-16><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|10<space|2spc>Sobolev
      Spaces> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-17><vspace|0.5fn>

      <with|par-left|<quote|1.5fn>|10.1<space|2spc>Campanato
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-18>>

      <with|par-left|<quote|1.5fn>|10.2<space|2spc>Sobolev
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-19>>

      <with|par-left|<quote|1.5fn>|10.3<space|2spc>Poincaré and Morrey
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-20>>

      <with|par-left|<quote|1.5fn>|10.4<space|2spc>BMO
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-21>>

      <with|par-left|<quote|1.5fn>|10.5<space|2spc>Imbeddings
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-22>>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|11<space|2spc>Scalar
      Elliptic Equations> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-23><vspace|0.5fn>

      <with|par-left|<quote|1.5fn>|11.1<space|2spc>Existence Theory
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-24>>

      <with|par-left|<quote|1.5fn>|11.2<space|2spc>Regularity
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-25>>

      <with|par-left|<quote|1.5fn>|11.3<space|2spc>Harnack Inequality Stuff
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-26>>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|12<space|2spc>Calculus
      of Variations> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-27><vspace|0.5fn>

      <with|par-left|<quote|1.5fn>|12.1<space|2spc>Quasiconvexity
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-28>>

      <with|par-left|<quote|1.5fn>|12.2<space|2spc>Null Lagrangians,
      Determinants <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-29>>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|13<space|2spc>Navier-Stokes
      Equations> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-30><vspace|0.5fn>
    </associate>
  </collection>
</auxiliary>