% Jean-Pierre Demailly
% Universit\'e de Grenoble I, Institut Fourier
% \'Equations Diff\'erentielles et Analyse Num\'erique
% Collection Grenoble-Sciences, Presses Universitaires de Grenoble

\newif\ifcolorized \colorizedfalse
\newif\iftocfirst \tocfirsttrue

\input hyperbasics.tex

\ifcolorized
\magnification=1200 \hsize=13cm\hoffset=0cm \vsize=18.8cm\voffset=0.3cm
\else
%% \magnification=1000 \hsize=13cm\hoffset=1.5cm \vsize=18.8cm\voffset=2cm
\magnification=1200 \hsize=13cm\hoffset=0cm \vsize=18.8cm\voffset=0.3cm
\fi

\pretolerance=500 \tolerance=1000 \brokenpenalty=5000
\parindent=0mm
\parskip=5pt plus 1pt minus 1pt
\baselineskip=12.5pt

\catcode`\@=11
\def\makeheadline{\vbox to\z@{\vskip-25.4\p@
  \line{\vbox to8.5\p@{}\the\headline}\vss}\nointerlineskip}
\def\makefootline{\baselineskip28\p@\lineskiplimit\z@\line{\the\footline}}
\catcode`\@=12

% new fonts definitions

\font\twentyfourss=cmss10 at 24.8832pt
\font\twentyss=cmss10 at 20.736pt
\font\seventeenss=cmss10 at 17.28pt
\font\fourteenss=cmss10 at 14.4pt
\font\twelvess=cmss10 at 12pt
\font\tenss=cmss10
\font\sevenss=cmss10 at 7pt

\font\twentyfourbss=cmssbx10 at 24.8832pt
\font\twentybss=cmssbx10 at 20.736pt
\font\seventeenbss=cmssbx10 at 17.28pt
\font\fourteenbss=cmssbx10 at 14.4pt
\font\twelvebss=cmssbx10 at 12pt
\font\tenbss=cmssbx10
\font\sevenbss=cmssbx10 at 7pt

\font\twelveline=line10 at 12pt
\font\twentyrm=cmr10 at 20.736pt
\font\seventeensy=cmsy10 at 17.28pt
\font\seventeenrm=cmr10 at 17.28pt
\font\seventeenbf=cmbx10 at 17.28pt
\font\seventeeni=cmmi10 at 17.28pt
\font\fourteensy=cmsy10 at 14.4pt
\font\fourteenrm=cmr10 at 14.4pt
\font\fourteenbf=cmbx10 at 14.4pt
\font\fourteeni=cmmi10 at 14.4pt
\font\twelvesy=cmsy10 at 12pt
\font\twelvebf=cmbx10 at 12pt
\font\twelvei=cmmi10 at 12pt
\font\eightbf=cmbx10 at 8pt
\font\sixbf=cmbx10 at 6pt

\font\twelvei=cmmi10 at 12pt
\font\eighti=cmmi10 at 8pt
\font\sixi=cmmi10 at 6pt

\font\twelverm=cmr10 at 12pt
\font\eightrm=cmr10 at 8pt
\font\sixrm=cmr10 at 6pt


\font\eightsy=cmsy8
\font\sixsy=cmsy6

\font\eightit=cmti8
\font\eighttt=cmtt8
\font\eightsl=cmsl8

\font\seventeenbsy=cmbsy10 at 17.28pt
\font\fourteenbsy=cmbsy10 at 14.4pt
\font\twelvebsy=cmbsy10 at 12pt
\font\tenbsy=cmbsy10
\font\eightbsy=cmbsy8
\font\sevenbsy=cmbsy7
\font\sixbsy=cmbsy6
\font\fivebsy=cmbsy5

\font\tenmsa=msam10
\font\eightmsa=msam8
\font\sevenmsa=msam7
\font\fivemsa=msam5
\newfam\msafam
  \textfont\msafam=\tenmsa
  \scriptfont\msafam=\sevenmsa
  \scriptscriptfont\msafam=\fivemsa
\def\msa{\fam\msafam\tenmsa}

\font\fourteenmsb=msbm10 at 14.4pt
\font\tenmsb=msbm10
\font\eightmsb=msbm8
\font\sevenmsb=msbm7
\font\fivemsb=msbm5
\newfam\msbfam
  \textfont\msbfam=\tenmsb
  \scriptfont\msbfam=\sevenmsb
  \scriptscriptfont\msbfam=\fivemsb
\def\Bbb{\fam\msbfam\tenmsb}

\font\fourteenCal=eusm10 at 14.4pt
\font\twelveCal=eusm10 at 12pt
\font\tenCal=eusm10
\font\sevenCal=eusm7
\font\fiveCal=eusm5
\newfam\Calfam
  \textfont\Calfam=\tenCal
  \scriptfont\Calfam=\sevenCal
  \scriptscriptfont\Calfam=\fiveCal
\def\Cal{\fam\Calfam\tenCal}

\font\teneuf=eusm10
\font\teneuf=eufm10
\font\seveneuf=eufm7
\font\fiveeuf=eufm5
\newfam\euffam
  \textfont\euffam=\teneuf
  \scriptfont\euffam=\seveneuf
  \scriptscriptfont\euffam=\fiveeuf
\def\euf{\fam\euffam\teneuf}

\font\seventeenbfit=cmmib10 at 17.28pt
\font\fourteenbfit=cmmib10 at 14.4pt
\font\twelvebfit=cmmib10 at 12pt
\font\tenbfit=cmmib10
\font\eightbfit=cmmib8
\font\sevenbfit=cmmib7
\font\sixbfit=cmmib6
\font\fivebfit=cmmib5
\newfam\bfitfam
  \textfont\bfitfam=\tenbfit
  \scriptfont\bfitfam=\sevenbfit
  \scriptscriptfont\bfitfam=\fivebfit
\def\bfit{\fam\bfitfam\tenbfit}

\def\gS{\hbox{\teneuf S}}

% changing font sizes
\catcode`\@=11
\def\eightpoint{%
  \textfont0=\eightrm \scriptfont0=\sixrm \scriptscriptfont0=\fiverm
  \def\rm{\fam\z@\eightrm}%
  \textfont1=\eighti \scriptfont1=\sixi \scriptscriptfont1=\fivei
  \def\mit{\fam\@ne\eighti}%
  \def\oldstyle{\fam\@ne\eighti}%
  \textfont2=\eightsy \scriptfont2=\sixsy \scriptscriptfont2=\fivesy
  \textfont\itfam=\eightit
  \def\it{\fam\itfam\eightit}%
  \textfont\slfam=\eightsl
  \def\sl{\fam\slfam\eightsl}%
  \textfont\bffam=\eightbf \scriptfont\bffam=\sixbf
  \scriptscriptfont\bffam=\fivebf
  \def\bf{\fam\bffam\eightbf}%
  \textfont\ttfam=\eighttt
  \def\tt{\fam\ttfam\eighttt}%
  \textfont\msbfam=\eightmsb
  \def\Bbb{\fam\msbfam\eightmsb}%
  \abovedisplayskip=9pt plus 2pt minus 6pt
  \abovedisplayshortskip=0pt plus 2pt
  \belowdisplayskip=9pt plus 2pt minus 6pt
  \belowdisplayshortskip=5pt plus 2pt minus 3pt
  \smallskipamount=2pt plus 1pt minus 1pt
  \medskipamount=4pt plus 2pt minus 1pt
  \bigskipamount=9pt plus 3pt minus 3pt
  \normalbaselineskip=9pt
  \setbox\strutbox=\hbox{\vrule height7pt depth2pt width0pt}%
  \let\bigf@ntpc=\eightrm \let\smallf@ntpc=\sixrm
  \normalbaselines\rm}
\catcode`\@=12

\def\eightpointbf{%
 \textfont0=\eightbf   \scriptfont0=\sixbf   \scriptscriptfont0=\fivebf
 \textfont1=\eightbfit \scriptfont1=\sixbfit \scriptscriptfont1=\fivebfit
 \textfont2=\eightbsy  \scriptfont2=\sixbsy  \scriptscriptfont2=\fivebsy
 \eightbf
 \baselineskip=10pt}

\def\tenpointbf{%
 \textfont0=\tenbf   \scriptfont0=\sevenbf   \scriptscriptfont0=\fivebf
 \textfont1=\tenbfit \scriptfont1=\sevenbfit \scriptscriptfont1=\fivebfit
 \textfont2=\tenbsy  \scriptfont2=\sevenbsy  \scriptscriptfont2=\fivebsy
 \bf}
        
\def\twelvepointbf{%
 \textfont0=\twelvebf   \scriptfont0=\eightbf   \scriptscriptfont0=\sixbf
 \textfont1=\twelvebfit \scriptfont1=\eightbfit \scriptscriptfont1=\sixbfit
 \textfont2=\twelvebsy  \scriptfont2=\eightbsy  \scriptscriptfont2=\sixbsy
 \twelvebss
 \baselineskip=14.4pt}

\def\fourteenpointbf{%
 \textfont0=\fourteenbf   \scriptfont0=\tenbf   \scriptscriptfont0=\sevenbf
 \textfont1=\fourteenbfit \scriptfont1=\tenbfit \scriptscriptfont1=\sevenbfit
 \textfont2=\fourteenbsy  \scriptfont2=\tenbsy  \scriptscriptfont2=\sevenbsy
 \fourteenbss
 \baselineskip=17.28pt}

\def\seventeenpointbf{%
 \textfont0=\seventeenbf  \scriptfont0=\twelvebf  \scriptscriptfont0=\eightbf
 \textfont1=\seventeenbfit\scriptfont1=\twelvebfit\scriptscriptfont1=\eightbfit
 \textfont2=\seventeenbsy \scriptfont2=\twelvebsy \scriptscriptfont2=\eightbsy
 \seventeenbss
 \baselineskip=20.736pt}

\def\twelvepoint{%
 \textfont0=\twelverm   \scriptfont0=\eightrm  \scriptscriptfont0=\sixrm
 \textfont1=\twelvei    \scriptfont1=\eighti   \scriptscriptfont1=\sixi
 \textfont2=\twelvesy   \scriptfont2=\eightsy  \scriptscriptfont2=\sixsy
 \twelvess
 \baselineskip=14.4pt}

\def\fourteenpoint{%
 \textfont0=\fourteenrm  \scriptfont0=\tenrm  \scriptscriptfont0=\sevenrm
 \textfont1=\fourteeni  \scriptfont1=\teni    \scriptscriptfont1=\seveni
 \textfont2=\fourteensy  \scriptfont2=\tensy  \scriptscriptfont2=\sevensy
 \fourteenss
 \baselineskip=17.28pt}

\def\seventeenpoint{%
 \textfont0=\seventeenrm \scriptfont0=\twelverm  \scriptscriptfont0=\eightrm
 \textfont1=\seventeeni  \scriptfont1=\twelvei   \scriptscriptfont1=\eighti
 \textfont2=\seventeensy \scriptfont2=\twelvesy  \scriptscriptfont2=\eightsy
 \seventeenss
 \baselineskip=20.736pt}

\def\twentypoint{%
 \textfont0=\twentyrm  \scriptfont0=\fourteenrm \scriptscriptfont0=\eightrm
 \textfont1=\twentyi   \scriptfont1=\fourteeni  \scriptscriptfont1=\eighti
 \textfont2=\twentysy \scriptfont2=\fourteensy  \scriptscriptfont2=\eightsy
 \twentyss
 \baselineskip=24.8832pt}
 
% main item macros

\newdimen\srdim \srdim=\hsize
\newdimen\irdim \irdim=\hsize
\def\NOSECTREF#1{\noindent\hbox to \srdim{\null\dotfill ???(#1)}}
\def\SECTREF#1{\noindent\hbox to \srdim{\csname REF\romannumeral#1\endcsname}}
\def\INDREF#1{\noindent\hbox to \irdim{\csname IND\romannumeral#1\endcsname}}
\newlinechar=`\^^J
\def\openauxfile{
  \immediate\openin1\jobname.aux
  \ifeof1
  \message{^^JCAUTION\string: you MUST run TeX a second time^^J}
  \let\sectref=\NOSECTREF \let\indref=\NOSECTREF
  \else
  \input \jobname.aux
  \message{^^JCAUTION\string: if the file has just been modified you may 
    have to run TeX twice^^J}
  \let\sectref=\SECTREF \let\indref=\INDREF
  \fi
  \message{to get correct page numbers displayed in Contents or Index 
    Tables^^J}
  \immediate\openout1=\jobname.aux
  \let\END=\end \def\end{\immediate\closeout1\END}}
        
\newcount\notthispage \notthispage=1

\newbox\titlebox   \setbox\titlebox\hbox{\hfil}
\newbox\sectionbox \setbox\sectionbox\hbox{\hfil}
\newbox\chapterbox \setbox\chapterbox\hbox{\hfil}
\newbox\bookbox \setbox\bookbox\hbox{\hfil}
\newbox\ebookbox \setbox\ebookbox\hbox{\hfil}

\def\fakeskip{$\rlap{\phantom{\vrule height18pt depth2pt width0pt}}$}
\def\blackline{\rlap{\vrule width \hsize height -4.7pt depth 5pt}}
\def\upperblackline{\rlap{\vrule width \hsize height 11pt depth -10.7pt}}

\def\foliotoc{\ifnum\pageno=\notthispage \hfil
           \else \ifodd\pageno \blackline%
           {\tenss \copy\chapterbox%
           \hfill\romannumeral\pageno}\else\blackline%
           {\tenss\romannumeral\pageno\hfill\copy\titlebox}\fi\fi}
\def\foliochap{\ifnum\pageno=\notthispage \hfil
           \else \ifodd\pageno \blackline%
           {\tenss \copy\sectionbox%
           \hfill\number\pageno}\else\blackline%
           {\tenss\number\pageno\hfill\copy\chapterbox}\fi\fi}
\def\folioref{\ifnum\pageno=\notthispage \hfil
           \else \ifodd\pageno \blackline%
           {\tenss \copy\chapterbox%
           \hfill\number\pageno}\else\blackline%
           {\tenss\number\pageno\hfill\copy\titlebox}\fi\fi}
\def\foliopapebook{\blackline{\copy\bookbox\hfill\copy\ebookbox}}

\footline={\hfil}

\let\forceheader\eject
\def\blankline{\phantom{}\hfil\vskip0pt}
\def\chapterspacing{\phantom{$\ $}\vskip2.5cm}
\def\titlerunning#1{\setbox\titlebox\hbox{\tenss #1}}
\def\chapterrunning#1{\notthispage=\pageno
    \setbox\chapterbox\hbox{\tenss #1}}
\def\title#1{\noindent\hfil$\smash{\hbox{\seventeenpointbf #1}}$\hfil
             \titlerunning{#1}\medskip}
\def\titleleft#1{\noindent$\smash{\hbox{\seventeenpointbf #1}}$\hfil
                 \titlerunning{#1}\medskip}

\def\supersection#1{%
  \par\vskip1.5cm\penalty -100 
  \vbox{\baselineskip=17.28pt\noindent{{\fourteenpointbf #1}}}
  \vskip3pt
  \penalty 500
  \titlerunning{#1}}

\newcount\numbersection \numbersection=-1
\def\sectionrunning#1{\setbox\sectionbox\hbox{\tenss #1}
  \write1{\string\def\string\REF
      \romannumeral\numbersection\string{%
      \noexpand#1 \string\dotfill\string\space\number\pageno\string}}}

\def\RGBColor#1#2{\special{color push rgb #1}#2\special{color pop}}

\ifcolorized
\def\maincolor#1{\RGBColor{0 0 0.5}{#1}}
\else
\def\maincolor#1{#1}
\fi
\def\colorstate#1{\noindent\maincolor{{\bf #1.}}}

\def\section#1{%
  \removelastskip
  \vskip1cm\penalty -100
  \vbox{\baselineskip=17.28pt\noindent{\maincolor{\seventeenpoint #1\vskip0pt}}}
  \vskip1pt
  \penalty 500
  \advance\numbersection by 1
  \sectionrunning{#1}}

\def\subsection#1{%
  \removelastskip
  \vskip0.5cm\penalty -100
  \vbox{\noindent{\maincolor{\fourteenpoint #1\vskip0pt}}}
  \penalty 500}

\newcount\numberindex \numberindex=0  
\def\index#1#2{%
  \advance\numberindex by 1
  \write1{\string\def \string\IND #1%
     \romannumeral\numberindex \string{%
     \noexpand#2 \string\dotfill \space \number\numbersection, 
     p.\string\ \space\number\pageno \string}}}

\newcount\numberchap 
\iftocfirst \numberchap=-2 \else \numberchap=-1 \fi

\def\CHpage{
   \immediate\write1{\string\def\string\CH \romannumeral\numberchap\string{%
   \number\pageno\string}}}

\def\chapterjump{
  \vfill\eject
  \ifodd\pageno \else {\headline={\hfil}\null\vskip0pt\vfill\eject} \fi
  \advance\numberchap by 1
  \CHpage
}

\def\numpage#1{~\dotfill~#1}

\newdimen\dpp
\newbox\claimbox \setbox\claimbox\hbox{\hfil}

\long\def\claim#1#2\endclaim{\par\vskip 5pt\noindent 
{{\tenpointbf\maincolor{#1.}}\ {\it #2}\vskip-18pt}
\strut\kern\hsize
\ifcolorized
\special{" gsave 0.8 0.8 0.1 setrgbcolor -1 -1 scale
0 6 moveto 7 0 rlineto 0 -3 rlineto -4 0 rlineto 0 -4 rlineto
-3 0 rlineto closepath fill grestore}%
\else
\special{" gsave 0.667 setgray -1 -1 scale
0 6 moveto 7 0 rlineto 0 -3 rlineto -4 0 rlineto 0 -4 rlineto
-3 0 rlineto closepath fill grestore}%
\fi
\par\vskip11pt}

\long\def\exo#1#2\endexo{\par\vskip 5pt\noindent 
{\tenpointbf\maincolor{#1.}}\ {\rm #2}\par}

\def\ceqno(#1){\eqno{{\maincolor{(#1)}}}}
\def\cleqno(#1){\leqno{{\maincolor{(#1)}}}}

\def\eqleqno#1\hfill#2{\leqno
\hbox to 0.0001pt{\rlap{\rlap{\maincolor{#1}}%
\kern\hsize\llap{\maincolor{#2}}}}}

\def\joinrel{\mathrel{\mkern-3.5mu}}
\def\llra{\relbar\joinrel\relbar\joinrel\longrightarrow}
\def\llmapsto{\mapstochar\relbar\joinrel\relbar\joinrel\longrightarrow}
\def\vlra#1|{\mathrel{\hbox to#1mm{\rightarrowfill}}}


% Usual sets of numbers  
\def\bC{{\Bbb C}}
\def\bK{{\Bbb K}}
\def\bN{{\Bbb N}}
\def\bP{{\Bbb P}}
\def\bQ{{\Bbb Q}}
\def\bR{{\Bbb R}}
\def\bZ{{\Bbb Z}}

% Calligraphic capital letters
\def\cA{{\Cal A}}
\def\cC{{\Cal C}}
\def\cD{{\Cal D}}
\def\cE{{\Cal E}}
\def\cF{{\Cal F}}
\def\cG{{\Cal G}}
\def\cI{{\Cal I}}
\def\cL{{\Cal L}}
\def\cN{{\Cal N}}
\def\cP{{\Cal P}}
\def\cR{{\Cal R}}
\def\cS{{\Cal S}}
\def\cT{{\Cal T}}
\def\cV{{\Cal V}}

\def\tC{\smash{\tilde C}}

\def\tvi{{\vrule height 10pt depth 5pt width 0pt}}
\def\tv{\tvi\vrule}
\def\hexnbr#1{\ifnum#1<10 \number#1\else
 \ifnum#1=10 A\else\ifnum#1=11 B\else\ifnum#1=12 C\else
 \ifnum#1=13 D\else\ifnum#1=14 E\else\ifnum#1=15 F\fi\fi\fi\fi\fi\fi\fi}
\def\msatype{\hexnbr\msafam}
\def\msbtype{\hexnbr\msbfam}

\mathchardef\smallsetminus="2\msbtype72   \let\ssm\smallsetminus
\mathchardef\restr="3\msatype16
\mathchardef\supsetneqq="3\msbtype25
\mathchardef\subsetneq="3\msbtype28
\mathchardef\supsetneq="3\msbtype29
\mathchardef\leqslant="3\msatype36
\mathchardef\geqslant="3\msatype3E
\mathchardef\complement="0\msatype7B
\let\ge=\geqslant
\let\le=\leqslant
\let\dsp=\displaystyle

\def\VVert{{|}\kern-1.2pt{|}\kern-1.2pt{|}}

\def\square{\maincolor{\hfill \hbox{
\vrule height 1.453ex  width 0.093ex  depth 0ex
\vrule height 1.5ex  width 1.3ex  depth -1.407ex\kern-0.1ex
\vrule height 1.453ex  width 0.093ex  depth 0ex\kern-1.35ex
\vrule height 0.093ex  width 1.3ex  depth 0ex}}}
\def\qed{\phantom{$\quad$}$\square$}

\def\frac#1#2{{#1\over #2}}
\def\dfrac#1#2{{#1\over #2}}
\def\fracs#1#2{{\scriptstyle{#1\over#2}}}
\def\fracss#1#2{{\scriptscriptstyle{#1\over#2}}}

\catcode`@=11
\def\ccmalign#1#2{\null\,\vcenter{\normalbaselines\m@th
    \ialign{$\displaystyle ##$\hfil&&\kern#1$\displaystyle ##$\hfil\crcr
      \mathstrut\crcr\noalign{\kern-\baselineskip}
      #2\crcr\mathstrut\crcr\noalign{\kern-\baselineskip}}}\,}
\def\cmalign#1{\null\,\vcenter{\normalbaselines\m@th
    \ialign{\hfil$##$&&$##$\hfil\crcr
      \mathstrut\crcr\noalign{\kern-\baselineskip}
      #1\crcr\mathstrut\crcr\noalign{\kern-\baselineskip}}}\,}
\catcode`@=12

\def\coloritem#1{\item{\maincolor{#1}}}

\def\itemv{%
\item{\llap{$\square\kern7pt$}\llap{$\raise0.6pt\hbox{$\times$}\kern6pt$}}}
\def\itemf{%
\item{\llap{$\square\kern7pt$}}}

\def\bul{$\scriptstyle\bullet$}

\def\bigzero{\hbox{~\seventeenrm 0~~~}}
\def\branch{\raise2pt\hbox{
   \vbox{\hbox{\vrule height 0.4pt depth 0pt width 12.2pt}
         \hbox{\raise-11.6pt\hbox{\twelveline\char"00}}}}}
\def\?{\hbox{$\,$}}
\def\lguil{\hbox{%
\raise1pt\hbox{$\scriptscriptstyle\langle\!\langle\kern1pt$}}}
\def\rguil{\hbox{%
\raise1pt\hbox{$\kern1pt\scriptscriptstyle\rangle\!\rangle$}}}

\def\sqind#1{\kern1.5pt\rlap{\raise5pt\hbox{$\scriptstyle#1$}}\kern-1.5pt}

\def\KH{{\rm KH}}
\def\Riemann{{\rm Riemann}}
\def\eq{\mathop{\rm =}}
\def\Bij{\mathop{\rm Bij}}
\def\tr{\mathop{\rm tr}}
\def\Log{\mathop{\rm Log}}
\def\Mat{\mathop{\rm Mat}}
\def\Ker{\mathop{\rm Ker}}
\def\Vect{\mathop{\rm Vect}}
\def\rang{\mathop{\rm rang}\nolimits}
\def\Inv{\mathop{\rm Inv}}
\def\Sp{\mathop{\rm Sp}\nolimits}
\def\Arccos{\mathop{\rm Arccos}}
\def\Arcsin{\mathop{\rm Arcsin}}
\def\Arctan{\mathop{\rm Arctan}}
\def\Argcosh{\mathop{\rm Argcosh}}
\def\Argsinh{\mathop{\rm Argsinh}}
\def\cotan{\mathop{\rm cotan}}
\def\cotanh{\mathop{\rm cotanh}}
\def\supess{\mathop{\rm sup\,ess}}
\def\infess{\mathop{\rm inf\,ess}}
\def\card{\mathop{\rm card}}
\def\Id{\mathop{\rm Id}\nolimits}
\def\Re{\mathop{\rm Re}}
\def\Im{\mathop{\rm Im}}
\def\Jac{\mathop{\rm Jac}}
\def\Supp{\mathop{\rm Supp}}
\def\oscil{\mathop{\rm oscil}}
\def\vol{\mathop{\rm vol}}
\def\aspect{\mathop{\rm aspect}}
\def\aire{\mathop{\rm aire}}
\def\longueur{\mathop{\rm longueur}}
\def\diam{\mathop{\rm diam}}
\def\sign{\mathop{\rm sign}}
\def\grad{\mathop{\rm grad}}
\def\rot{\mathop{\rm rot}}
\def\div{\mathop{\rm div}}
\def\Var{\mathop{\rm Var}}
\def\pf{p\kern-1pt f}
\def\dplus{\mathrel{\dot +}}

\def\note#1#2#3{\footnote{}%
{\baselineskip=8pt\leftskip=5.4mm\rlap{\strut\kern-5.4mm${}^{\rm(#2)}$}%
\hyperdef#1{}{}{}{\eightpoint #3}\vskip-15pt}%
${}^{\rm(\hyperref#1{#2})}$}
\def\\{\hfil\break}
\def\demi{\textstyle{1\over 2}}
\def\ovl{\overline}
\def\ovr{\overrightarrow}
\def\ul#1{$\underline{\smash{\hbox{#1}}}$}

\def\build#1^#2_#3{\mathop{#1}\limits^{#2}_{#3}}
\def\demo{\noindent{\it D\'emonstration.}\ }

\newbox\formulabox  \setbox\formulabox\hbox{\hfil}
\newdimen\wdd \newdimen\htt
\def\boxed#1#2{\setbox\formulabox\hbox{$\displaystyle #2$}
\wdd=8pt \advance \wdd by \wd\formulabox
\htt=4pt \advance \htt by \ht\formulabox
\dpp=4pt \advance \dpp by \dp\formulabox
\ifcolorized
\smash{\rlap{\hbox{\RGBColor{1 0.92 0.94}{\kern-4pt\vrule width \wdd height \htt depth \dpp}}}}\else%
\smash{\rlap{\hbox{\RGBColor{0.85 0.85 0.85}{\kern-4pt\vrule width \wdd height \htt depth \dpp}}}}\fi%
\copy\formulabox}
\def\boxit#1#2{\hbox{\vrule
 \vbox{\hrule\kern#1
  \vtop{\hbox{\kern#1 #2\kern#1}%
   \kern#1\hrule}}%
 \vrule}}
\def\boxmat#1#2{\boxit{#1}{$#2$}} 

% inclusion of PostScript files
\special{header=mdrlib.ps}

\long\def\InsertFig#1 #2 #3 #4\EndFig{\par
\hbox{\hskip #1mm$\vbox to#2mm{\vfil\special{" 
#3}}#4$}}

\long\def\InsertPSFile#1 #2 #3 #4 #5 #6\EndFig{\par
\hbox{\hskip #1mm$\vbox to#2mm{%
\vfil\special{psfile=#5 hscale=#3 vscale=#4}}#6$}}

\long\def\InsertImage#1 #2 #3 #4 #5 #6 #7 #8\EndFig{\par
\hbox{\hskip #1mm$\vbox to#2mm{%
\vfil\special{psfile="`img2eps file #7 height #4 mm width #3 mm gamma #5
angle #6}}#8$}}

\long\def\LabelTeX#1 #2 #3\ELTX{\rlap{\kern#1mm\raise#2mm\hbox{#3}}}

\def\psrotate#1#2{\special{ps: /a{pop pop}def gsave #1 neg rotate}#2
\special{ps: grestore /a{moveto}bind def}}

% bibibliography
\def\bibitem#1&#2&#3&#4&%
{\hangindent=1.2cm\hangafter=1
\noindent\rlap{\hbox{\eightpointbf #1}}\kern1.2cm{\rm #2}{\it #3}{\rm #4.}} 

%%% Texte principal

\setbox\bookbox\hbox{\tenss 
\maincolor{%
Analyse num\'erique et \'equations diff\'erentielles}}
\setbox\ebookbox\hbox{\tenss 
\maincolor{%
\href{https://grenoble-sciences.ujf-grenoble.fr/pap-ebook/demailly/6-equations-differentielles-resolution}{\RGBColor{0 0 1}{$\strut\,\Rightarrow$ site web compagnon$\,$}}}}

\footline={\upperblackline
\maincolor{%
\tenss Solution des probl\`emes du chap.~6\hfill
\copyright\ Grenoble Sciences, J.-P.~Demailly $\scriptstyle \bullet$ \folio}}

\openauxfile
\headline={\foliopapebook}

\notthispage=\pageno
\strut\vskip4mm
\maincolor{%
\title{Solution des probl\`emes du chapitre 6}}
\titlerunning{\maincolor{%
Solutions des probl\`emes du chapitre 6}}
\vskip10mm

\maincolor{{\bf 5.11.}} On \'etudie ici une m\'ethode alternative pour 
d\'eterminer l'\'equation de la cha\^{\i}nette.\\
Dans un plan vertical $Oxy$ assimil\'e \`a $\bR^2$ (o\`u $Oy$ est la
direction verticale pointant vers le haut), on consid\`ere un fil souple de 
longueur $\ell > 0$ dont les extr\'emit\'es sont fix\'ees \`a des points 
de m\^eme hauteur, not\'es $(-a,0)$ et $(a,0)$, avec $a>0$. 
La position d'\'equilibre 
du fil est le graphe d'une fonction $f$ qui minimise l'\'energie 
potentielle $P(f)=\smash{\int gy\,dm}$ sous la contrainte de longueur
$L(f) = \smash{\int ds}=\ell$. On a plus pr\'ecis\'ement
$ds=\sqrt{dx^2+dy^2}$, $dm=\mu\,ds$ o\`u $\mu$ est la masse lin\'eique,
d'o\`u
$$
  L(f) = \int_{-a}^a \sqrt{1+f'(x)^2}\,dx~, \qquad
  P(f) = \mu g\int_{-a}^a f(x)\sqrt{1+f'(x)^2}\,dx~, \qquad
$$
On consid\`ere l'espace de Banach
$$
  E \,=\, \Bigl\{f \in \cC^1([-a,a],\bR)\,\Big|\, f(-a) = f(a) = 0
  \Bigr\}~,
$$
muni de la norme
$$
  \|f\| \,=\, \sup_{x\in[-a,a]}|f(x)| + \sup_{x\in[-a,a]}|f'(x)|~.
$$
Les fonctionnelles $L$, $P$ d\'efinissent des applications 
$L,\,P : E \to \bR$.

\maincolor{{\bf(a)}} Calculons les diff\'e\-rentielles respectives de
$L$ et $P$ en tout point $f \in E$. Pour cela, on \'evalue
$L(f+h)-L(f)$ et $P(f+h)-P(f)$ \`a l'aide de la formule de Taylor \`a
l'ordre~$2$. On est amen\'e \`a estimer les diff\'erences
$$
\eqalign{
\delta_L(x)=&\sqrt{1+(f'(x)+h'(x))^2}-\sqrt{1+f'(x)^2},\cr
\delta_P(x)=&(f(x)+h(x))\sqrt{1+(f'(x)+h'(x))^2}-f(x)\sqrt{1+f'(x)^2}.\cr}
$$
Pour la diff\'erence $\delta_L(x)$, la fonction $g(t)=\sqrt{1+t^2}$ v\'erifie
$g'(t)={t\over\sqrt{1+t^2}}$ et $g''(t)={1\over(1+t^2)^{3/2}}$, de sorte
que $g''(t)\le 1$ pour tout $t$. La premi\`ere diff\'erence s'\'ecrit alors
$$
\delta_L(x)=h'(x)g'(f'(x))+{1\over 2}h'(x)^2g''\big(f'(x)+\theta(x)h'(x)\big)
$$
avec $\theta(x)\in{}]0,1[$. On a d'autre part
$$
\delta_P(x)=h(x)\sqrt{1+f'(x)^2}+f(x)\delta_L(x)+h(x)\delta_L(x).
$$
En s\'eparant les termes lin\'eaires et les termes d'erreur, on en d\'eduit
$$
\eqalign{
dL(f)(h) &= \int_{-a}^a {h'(x)f'(x)\over\sqrt{1+f'(x)^2}}\,dx~, \cr
dP(f)(h) &= \mu g\int_{-a}^a h(x)\sqrt{1+f'(x)^2}\,dx+
f(x){h'(x)f'(x)\over\sqrt{1+f'(x)^2}}\,dx,\cr}
$$
avec des termes d'erreur qui valent
$$
\eqalign{
&\hbox{(pour $dL$)}\quad\int_{-a}^a{1\over 2}
h'(x)^2g''\big(f'(x)+\theta(x)h'(x)\big)\,dx\cr
&\hbox{(pour $dP$)}\quad\int_{-a}^a\Big(f(x)
{1\over 2}h'(x)^2g''\big(f'(x)+\theta(x)h'(x)\big)+
h(x)\delta_L(x)\Big)dx\cr}
$$
Dans le cas de $dL$, on voit gr\^ace \`a l'estimation $0\le g''\le 1$ que
le terme d'erreur est major\'e par
$$
2a\times {1\over 2}\Vert h'\Vert_\infty^2\le a\Vert h\Vert^2=o(\Vert h\Vert).
$$
Dans le cas de $dP$, on a un premier terme d'erreur qui est major\'e par
$a\Vert f\Vert\,\Vert h\Vert^2$ d'apr\`es ce qui pr\'ec\`ede, et un
terme suppl\'ementaire $\int_{-a}^ah(x)\delta_L(x)\,dx$. Or
$h(x)\delta_L(x)$ est la somme d'un terme
$h(x){h'(x)f'(x)\over\sqrt{1+f'(x)^2}}$ major\'e par $\Vert h\Vert^2\Vert f\Vert$,
plus un nouveau terme d'erreur encore \lguil\?plus petit\?\rguil
$$
\int_{-a}^ah(x){1\over 2}h'(x)^2g''\big(f'(x)+\theta(x)h'(x)\big)\,dx
$$
major\'e par $a\Vert h\Vert^3$. On en d\'eduit que $L$ et $P$ sont bien
diff\'erentiables, avec les expressions donn\'ees pr\'ec\'edemment pour les
diff\'erentielles. Sous l'hypoth\`ese que $f$ est de classe $\cC^2$, et
compte tenu des conditions $h(a)=h(-a)=0$ pour $h\in E$, on
peut effectuer une int\'egration par parties pour obtenir
$$
\eqalign{
dL(f)(h) &= -\int_{-a}^a h(x){d\over dx}\left(
{f'(x)\over\sqrt{1+f'(x)^2}}\right)\,dx~, \cr
dP(f)(h) &= \mu g\int_{-a}^a h(x)\sqrt{1+f'(x)^2}\,dx-
h(x){d\over dx}\left({f(x)f'(x)\over\sqrt{1+f'(x)^2}}\right)\,dx.\cr}
$$

\maincolor{{\bf(b)}} On suppose $\ell > 2a$. la cha\^{\i}nette r\'ealise le
minimum de $P(f)$ sous la contrainte $L(f) = \ell$, pour $f\in E$.
Il s'agit d'un probl\`eme \lguil\?d'extremas li\'es\?\rguil. Si l'infimum existe
et s'il est atteint en un point $f \in E \cap \cC^2([-a,a])$, la condition
$dL(f)(h)=0$ implique $dP(f)(h)=0$, donc les formes lin\'eaires associ\'ees
satisfont la condition de proportionnalit\'e 
$dP(f)(h)=\lambda dL(f)(h)$ ($\lambda\in\bR$), d'apr\`es le th\'eor\`eme 
des multiplicateurs de Lagrange. Comme cette \'egalit\'e doit
\^etre vraie pour toute fonction~$h$, on en d\'eduit que $f$ 
satisfait l'\'equation diff\'erentielle
$$
\sqrt{1+f'(x)^2}-
{d\over dx}\left({f(x)f'(x)\over\sqrt{1+f'(x)^2}}\right)=
c\,{d\over dx}\left({f'(x)\over\sqrt{1+f'(x)^2}}\right)
$$
avec $c={\lambda\over\mu g}$. En omettant la variable $x$ pour simplifier,
on trouve la condition
$$
\sqrt{1+f^{\prime 2}}-{ff''+f^{\prime 2}\over\sqrt{1+f^{\prime 2}}}
+{ff^{\prime 2}f''\over(1+f^{\prime 2})^{3/2}}=
c\left({f''\over \sqrt{1+f^{\prime 2}}}-
{f^{\prime 2}f''\over(1+f^{\prime 2})^{3/2}}\right).
$$
En  multipliant cette expression par $(1+f^{\prime 2})^{3/2}$, on obtient
la condition \'equivalente
$$
(1+f^{\prime 2})^2-(1+f^{\prime 2})(ff''+f^{\prime 2})+
ff^{\prime 2}f''=c\big(f''(1+f^{\prime 2})-f^{\prime 2}f''\big).
$$
Ceci se simplifie (quelque peu miraculeusement~!) en
$$
(1+f^{\prime 2})(1-ff'')+ff^{\prime 2}f''=cf'',
$$
soit encore
$$
cf''=1+f^{\prime 2}-ff''.
$$
La fonction $g=f+c$ v\'erifie alors $gg''=1+g^{\prime 2}$, ce qui peut s'\'ecrire
aussi
$$
{g'g''\over 1+g^{\prime 2}}={g'\over g}
$$
(on notera que $g$ ne peut s'annuler). Cette \'egalit\'e 
s'int\`egre sous la forme
$$
{1\over 2}\ln(1+g^{\prime 2})=\ln|g|+C,
$$
soit encore $1+g^{\prime 2}=b^2g^2$ o\`u $b=e^C$ est une 
constante positive. On a une solution constante $g(x)={1\over b}$,
$g'(x)=0$, tandis que pour $g'\ne 0$ il vient
$$
{g'\over\sqrt{b^2g^2-1}}=\varepsilon=\pm 1,
$$
par cons\'equent ${1\over b}\arg\cosh(b|g(x)|)=\pm x$.
Il s'ensuit que $g(x)=\pm{1\over b}\cosh(bx)$,~d'o\`u
$$
f(x)=\pm{1\over b}\cosh(bx)-c.
$$
On notera que la continuit\'e des fonctions mises en jeu impose que le signe
$\pm$ ne d\'epende pas de la valeur de la variable $x$).
La condition $f(a)=0$ implique 
$$
f(x)=\pm{1\over b}\big(\cosh(bx)-\cosh(ab)\big).
$$
La solution $g$ constante donnerait $f(x)=0$ identiquement, d'o\`u une
longueur $L(f)=2a$, alors qu'on a suppos\'e $\ell>2a$.
Le signe initial figurant dans $f(x)$ doit \'egalement \^etre $+1$, sinon 
on obtient une solution non physiquement
r\'ealiste $f(x)>0$ pour $x\in{}]-a,a[$, qui correspond en fait \`a un maximum
de $P(f)$ au lieu d'un minimum (noter que par changement de $f$ en $-f$, on a 
$P(-f)=-P(f)$).
Il reste \`a relier les constantes
$a,b,\ell$. Pour cela on calcule
$$
L(f)=\int_{-a}^a\!\sqrt{1+f'(x)^2}\,dx=
\int_{-a}^a\!\sqrt{1+\sinh(bx)^2}\,dx=
\int_{-a}^a\cosh(bx)\,dx={2\over b}\sinh(ab),
$$
et la contrainte $L(f)=\ell$ donne la solution cherch\'ee
$$
f(x)={1\over b}\big(\cosh(bx)-\cosh(ab)\big)\quad\hbox{o\`u $b$ est tel que}
\quad {2\over b}\arg\sinh(ab)=\ell.
$$
Il existe une valeur unique de $b$ v\'erifiant la derni\`ere condition, car
par convexit\'e de $\sinh$ sur $\bR_+$, la fonction $b\mapsto {2\over b}
\sinh(ab)$ est strictement croissante sur l'intervalle $]0,+\infty[$, 
et y varie de $2a$ \`a $+\infty$.
\medskip

\colorstate{5.12${}^{**}$} On appelle m\'etrique de Poincar\'e du disque 
unit\'e $D=\{|z|<1\}$ du plan complexe la m\'etrique riemannienne
$$
ds^2={|dz|^2\over (1-|z|^2)^2}={dx^2+dy^2\over(1-(x^2+y^2))^2},\qquad
z=x+i y.
$$
\maincolor{{\bf(a)}} Avec les notations du \S\?4.4, l'\'energie d'un chemin
$\gamma:[a,b]\to D$ s'exprime sous la forme
$$
E(\gamma)=\int_a^b\Vert\gamma'(t)\Vert_q^2\,dt
=\int_a^b{|\gamma'(t)|^2\over (1-|\gamma(t)|^2)^2}\,dt
=\int_a^b{\gamma'(t)\overline{\gamma'(t)}\over (1-\gamma(t)
\overline{\gamma(t)})^2}\,dt.
$$
En changeant $\gamma$ en $\gamma+h$, on voit que la diff\'erentielle de
l'\'energie est donn\'ee par
$$
E(\gamma)\cdot h=2\Re\int_a^b\left({\gamma'(t)\over
\big(1-\gamma(t)\overline{\gamma(t)}\big)^2}\overline{h'(t)}
+2{\gamma(t)\gamma'(t)\overline{\gamma'(t)}\over
\big(1-\gamma(t)\overline{\gamma(t)}\big)^3}\overline{h(t)}\right)dt,
$$
en effet les termes apparaissent par paires de termes conjugu\'es et 
$z+\overline z=2\,\Re z$. En int\'egrant
par parties le premier terme de l'int\'egrande et en utilisant le fait que
$h(a)=h(b)=0$, il~vient
$$
E(\gamma)\cdot h=2\Re\int_a^b\left(-{d\over dt}{\gamma'(t)\over
\big(1-\gamma(t)\overline{\gamma(t)}\big)^2}\overline{h(t)}
+2{\gamma(t)\gamma'(t)\overline{\gamma'(t)}\over
\big(1-\gamma(t)\overline{\gamma(t)}\big)^3}\overline{h(t)}\right)dt.
$$
Comme ceci doit \^etre vrai pour tout $h$, les g\'eod\'esiques, \`a savoir les
points critiques de l'\'energie, sont les solutions de l'\'equation 
d'Euler-Lagrange 
$$
-{d\over dt}{\gamma'(t)\over
\big(1-\gamma(t)\overline{\gamma(t)}\big)^2}
+2{\gamma(t)\gamma'(t)\overline{\gamma'(t)}\over
\big(1-\gamma(t)\overline{\gamma(t)}\big)^3}=0
$$
En d\'eveloppant la d\'eriv\'ee du premier terme et en multipliant par
$\big(1-\gamma(t)\overline{\gamma(t)}\big)^2$, on voit apr\`es simplification 
que l'\'equation des g\'eod\'esiques peut se r\'ecrire
$$
\gamma''(t)+{2\gamma'(t)^2\overline{\gamma(t)}\over
1-\gamma(t)\overline{\gamma(t)}}=0.\leqno(*)
$$

\maincolor{{\bf(b)}} Il est facile de voir que le chemin 
$\gamma(t)=\tanh(kt)$ (qui d\'ecrit le
diam\`etre $]-1,1[$ du disque) est solution de l'\'equation pour
tout $k\in\bR^*_+$, en effet $\gamma'(t)=k(1-\tanh(kt)^2)$ et
$\gamma''(t)=-2k^2\tanh(kt)(1-\tanh(kt)^2)$.
\medskip

\maincolor{{\bf(c)}} Si $t\mapsto \gamma(t)$ est solution, alors 
$t\mapsto\lambda\gamma(t)$ est encore solution pour tout nombre complexe 
$\lambda$ de module $1$. En effet par substitution $\gamma\mapsto
\lambda\gamma$ dans $(*)$, $\gamma''(t)$ est multipli\'e par $\lambda$,
$2\gamma'(t)^2\overline{\gamma(t)}$ est multipli\'e par $\lambda^2
\overline{\lambda}=\lambda$, et $1-\gamma(t)\overline{\gamma(t)}$ 
reste inchang\'e.

On v\'erifie maintenant que pour toute solution $\gamma$ de $(*)$
et toute homographie complexe  $h_a$ de la forme
$$h_a(z)={z+a\over 1+\overline a z},\qquad\hbox{$a\in D$},$$
la compos\'ee $h_a\circ\gamma$ est \'egalement solution de $(*)$.
Ceci d\'emontrera en fait que les g\'eod\'esiques sont invariantes
par les automorphismes holomorphes $z\mapsto\lambda{z+a\over 1+\overline a z}$
du disque unit\'e. On a d'abord
$$
\eqalign{
h_a'(z)&={(1+\overline a z)-\overline a(z+a)\over(1+\overline a z)^2}=
{1- |a|^2\over (1+\overline a z)^2},\quad
h_a''(z)=-2\overline a{1- |a|^2\over (1+\overline a z)^3}\cr
1-|h_a(z)|^2&={|1+\overline a z|^2-|z+a|^2\over |1+\overline a z|^2}=
{(1-|a|^2)(1-|z|^2)\over |1+\overline a z|^2},\cr}
$$
car les doubles produits $\overline a z+a\overline z$ figurant au 
num\'erateur s'auto-d\'etruisent. Ceci donne
$$
\eqalign{
(h_a\circ\gamma)'(t)&=h_a'(\gamma(t))\gamma'(t)={(1-|a|^2)\gamma'(t)\over
(1+\overline a\gamma(t))^2},\cr
(h_a\circ\gamma)''(t)&=h_a'(\gamma(t))\gamma''(t)+h_a''(\gamma(t))\gamma'(t)^2\cr
&={1- |a|^2\over (1+\overline a\gamma(t))^2}\gamma''(t)
-2\overline a{1- |a|^2\over (1+\overline a \gamma(t))^3}\gamma'(t)^2,\cr
{(h_a\circ\gamma)'(t)^2\overline{(h_a\circ\gamma)(t)}\over
1-|h_a\circ\gamma(t)|^2}
&={(1-|a|^2)^2\gamma'(t)^2(\overline{\gamma(t)}+\overline a)\over
(1+\overline a\gamma(t))^4(\overline{1+\overline a\gamma(t)})}
{|1+\overline a \gamma(t)|^2\over(1-|a|^2)(1-|\gamma(t)|^2)}\cr
&={(1-|a|^2)\gamma'(t)^2(\overline{\gamma(t)}+\overline a)\over
(1+\overline a\gamma(t))^3(1-|\gamma(t)|^2)},
\cr}
$$
de sorte qu'on a bien (apr\`es un calcul p\'enible !)
$$
(h_a\circ\gamma)''(t)+
{2(h_a\circ\gamma)'(t)^2\overline{(h_a\circ\gamma)(t)}\over
1-|h_a\circ\gamma(t)|^2}=0.
$$

\maincolor{{\bf(d)}} On rappelle que les automorphismes holomorphes du
disque unit\'e sont les homographies pr\'eservant le cercle unit\'e, et sont
donn\'es par l'une ou l'autre des formules
$$
\varphi(z)=h_a(\lambda z)=\lambda h_b(z)\quad\hbox{avec $a=\lambda b$,
$|\lambda|=1$}.
$$
Consid\'erons le chemin tel que
$\gamma(t)=h_a(\lambda\tanh(kt))$, $a\in D$, $|\lambda|=1$, $k\in\bR_+^*$.
D'apr\`es la question (c), c'est un arc de g\'eod\'esique. De plus
$\gamma(0)=h_a(0)=a$ et $\gamma'(0)=h_a'(0)\times\lambda k=\lambda k(1-|a|^2)$.
On voit ainsi que $\gamma(0)=a$ est un point quelconque de $D$ et que
$\gamma'(0)$ est un point quelconque de $\bC$, puisque $k$ et
$\lambda$ peuvent \^etre choisis pour ajuster respectivement le module
et l'argument. Le th\'eor\`eme d'unicit\'e de Cauchy-Lipschitz s'applique
car les coefficients de $(*)$ sont de classe $\cC^1$ (et m\^eme $\cC^\infty$).
On a donc bien ainsi trouv\'e l'unique g\'eod\'esique de donn\'ees initiales
$\gamma(0)=a\in D$ et $\gamma'(0)=b\in\bC$, il n'y en a pas d'autres.

\maincolor{{\bf(e)}} D'apr\`es (c) et (d) les trajectoires des
g\'eod\'esiques sont les images par les homographies $h_a$ des
diam\`etres du disque. Les diam\`etres du disque sont orthogonaux aux
cercle unit\'e $|z|=1$, et comme les automorphismes du disque sont des
transformations conformes, l'orthogonalit\'e est pr\'eserv\'ee. Toutes les
g\'eod\'esiques sont donc orthogonales au cercle unit\'e. Mais il est
classique (et on peut v\'erifier ais\'ement) que l'image d'une droite par
une homographie est une droite ou un cercle, les g\'eod\'esiques sont donc
des arcs de cercle orthogonaux au cercle unit\'e.


\end

% Local Variables:
% TeX-command-default: "uTeX"
% End:

