\batchmode

\documentclass[a4j,12pt]{jarticle}
\RequirePackage{ifthen}


\usepackage{Dennou6}
\usepackage{ascmac}




\usepackage[dvips]{color}


\pagecolor[gray]{.7}

\usepackage[]{inputenc}



\makeatletter

\makeatletter
\count@=\the\catcode`\_ \catcode`\_=8 
\newenvironment{tex2html_wrap}{}{}%
\catcode`\<=12\catcode`\_=\count@
\newcommand{\providedcommand}[1]{\expandafter\providecommand\csname #1\endcsname}%
\newcommand{\renewedcommand}[1]{\expandafter\providecommand\csname #1\endcsname{}%
  \expandafter\renewcommand\csname #1\endcsname}%
\newcommand{\newedenvironment}[1]{\newenvironment{#1}{}{}\renewenvironment{#1}}%
\let\newedcommand\renewedcommand
\let\renewedenvironment\newedenvironment
\makeatother
\let\mathon=$
\let\mathoff=$
\ifx\AtBeginDocument\undefined \newcommand{\AtBeginDocument}[1]{}\fi
\newbox\sizebox
\setlength{\hoffset}{0pt}\setlength{\voffset}{0pt}
\addtolength{\textheight}{\footskip}\setlength{\footskip}{0pt}
\addtolength{\textheight}{\topmargin}\setlength{\topmargin}{0pt}
\addtolength{\textheight}{\headheight}\setlength{\headheight}{0pt}
\addtolength{\textheight}{\headsep}\setlength{\headsep}{0pt}
\setlength{\textwidth}{349pt}
\newwrite\lthtmlwrite
\makeatletter
\let\realnormalsize=\normalsize
\global\topskip=2sp
\def\preveqno{}\let\real@float=\@float \let\realend@float=\end@float
\def\@float{\let\@savefreelist\@freelist\real@float}
\def\liih@math{\ifmmode$\else\bad@math\fi}
\def\end@float{\realend@float\global\let\@freelist\@savefreelist}
\let\real@dbflt=\@dbflt \let\end@dblfloat=\end@float
\let\@largefloatcheck=\relax
\let\if@boxedmulticols=\iftrue
\def\@dbflt{\let\@savefreelist\@freelist\real@dbflt}
\def\adjustnormalsize{\def\normalsize{\mathsurround=0pt \realnormalsize
 \parindent=0pt\abovedisplayskip=0pt\belowdisplayskip=0pt}%
 \def\phantompar{\csname par\endcsname}\normalsize}%
\def\lthtmltypeout#1{{\let\protect\string \immediate\write\lthtmlwrite{#1}}}%
\newcommand\lthtmlhboxmathA{\adjustnormalsize\setbox\sizebox=\hbox\bgroup\kern.05em }%
\newcommand\lthtmlhboxmathB{\adjustnormalsize\setbox\sizebox=\hbox to\hsize\bgroup\hfill }%
\newcommand\lthtmlvboxmathA{\adjustnormalsize\setbox\sizebox=\vbox\bgroup %
 \let\ifinner=\iffalse \let\)\liih@math }%
\newcommand\lthtmlboxmathZ{\@next\next\@currlist{}{\def\next{\voidb@x}}%
 \expandafter\box\next\egroup}%
\newcommand\lthtmlmathtype[1]{\gdef\lthtmlmathenv{#1}}%
\newcommand\lthtmllogmath{\lthtmltypeout{l2hSize %
:\lthtmlmathenv:\the\ht\sizebox::\the\dp\sizebox::\the\wd\sizebox.\preveqno}}%
\newcommand\lthtmlfigureA[1]{\let\@savefreelist\@freelist
       \lthtmlmathtype{#1}\lthtmlvboxmathA}%
\newcommand\lthtmlpictureA{\bgroup\catcode`\_=8 \lthtmlpictureB}%
\newcommand\lthtmlpictureB[1]{\lthtmlmathtype{#1}\egroup
       \let\@savefreelist\@freelist \lthtmlhboxmathB}%
\newcommand\lthtmlpictureZ[1]{\hfill\lthtmlfigureZ}%
\newcommand\lthtmlfigureZ{\lthtmlboxmathZ\lthtmllogmath\copy\sizebox
       \global\let\@freelist\@savefreelist}%
\newcommand\lthtmldisplayA{\bgroup\catcode`\_=8 \lthtmldisplayAi}%
\newcommand\lthtmldisplayAi[1]{\lthtmlmathtype{#1}\egroup\lthtmlvboxmathA}%
\newcommand\lthtmldisplayB[1]{\edef\preveqno{(\theequation)}%
  \lthtmldisplayA{#1}\let\@eqnnum\relax}%
\newcommand\lthtmldisplayZ{\lthtmlboxmathZ\lthtmllogmath\lthtmlsetmath}%
\newcommand\lthtmlinlinemathA{\bgroup\catcode`\_=8 \lthtmlinlinemathB}
\newcommand\lthtmlinlinemathB[1]{\lthtmlmathtype{#1}\egroup\lthtmlhboxmathA
  \vrule height1.5ex width0pt }%
\newcommand\lthtmlinlineA{\bgroup\catcode`\_=8 \lthtmlinlineB}%
\newcommand\lthtmlinlineB[1]{\lthtmlmathtype{#1}\egroup\lthtmlhboxmathA}%
\newcommand\lthtmlinlineZ{\egroup\expandafter\ifdim\dp\sizebox>0pt %
  \expandafter\centerinlinemath\fi\lthtmllogmath\lthtmlsetinline}
\newcommand\lthtmlinlinemathZ{\egroup\expandafter\ifdim\dp\sizebox>0pt %
  \expandafter\centerinlinemath\fi\lthtmllogmath\lthtmlsetmath}
\newcommand\lthtmlindisplaymathZ{\egroup %
  \centerinlinemath\lthtmllogmath\lthtmlsetmath}
\def\lthtmlsetinline{\hbox{\vrule width.1em \vtop{\vbox{%
  \kern.1em\copy\sizebox}\ifdim\dp\sizebox>0pt\kern.1em\else\kern.3pt\fi
  \ifdim\hsize>\wd\sizebox \hrule depth1pt\fi}}}
\def\lthtmlsetmath{\hbox{\vrule width.1em\kern-.05em\vtop{\vbox{%
  \kern.1em\kern0.8 pt\hbox{\hglue.17em\copy\sizebox\hglue0.8 pt}}\kern.3pt%
  \ifdim\dp\sizebox>0pt\kern.1em\fi \kern0.8 pt%
  \ifdim\hsize>\wd\sizebox \hrule depth1pt\fi}}}
\def\centerinlinemath{%
  \dimen1=\ifdim\ht\sizebox<\dp\sizebox \dp\sizebox\else\ht\sizebox\fi
  \advance\dimen1by.5pt \vrule width0pt height\dimen1 depth\dimen1 
 \dp\sizebox=\dimen1\ht\sizebox=\dimen1\relax}

\def\lthtmlcheckvsize{\ifdim\ht\sizebox<\vsize 
  \ifdim\wd\sizebox<\hsize\expandafter\hfill\fi \expandafter\vfill
  \else\expandafter\vss\fi}%
\providecommand{\selectlanguage}[1]{}%
\makeatletter \tracingstats = 1 


\begin{document}
\pagestyle{empty}\thispagestyle{empty}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength hsize=\the\hsize}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength vsize=\the\vsize}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength hoffset=\the\hoffset}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength voffset=\the\voffset}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength topmargin=\the\topmargin}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength topskip=\the\topskip}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength headheight=\the\headheight}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength headsep=\the\headsep}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength parskip=\the\parskip}\lthtmltypeout{}%
\lthtmltypeout{latex2htmlLength oddsidemargin=\the\oddsidemargin}\lthtmltypeout{}%
\makeatletter
\if@twoside\lthtmltypeout{latex2htmlLength evensidemargin=\the\evensidemargin}%
\else\lthtmltypeout{latex2htmlLength evensidemargin=\the\oddsidemargin}\fi%
\lthtmltypeout{}%
\makeatother
\setcounter{page}{1}
\onecolumn

% !!! IMAGES START HERE !!!



\setlength{\parskip}{4mm}%

\setlength{\parskip}{4mm}
\stepcounter{section}
\stepcounter{section}
\stepcounter{subsection}
\stepcounter{subsection}
\stepcounter{subsection}
\stepcounter{subsection}
\stepcounter{section}
\stepcounter{subsection}
{\newpage\clearpage
\lthtmldisplayA{displaymath1103}%
\begin{displaymath}
          \tilde{\rho }_{i} = \rho _{i}^{n}-\frac{\Delta t}{\Delta x}
          [\rho _{i}^{n}u_{i}^{n}-\rho _{i-1}^{n}u_{i-1}^{n}].
        \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1104}%
\begin{displaymath}
          \Delta _{i+\frac{1}{2}} = \tilde{\rho }_{i+1}-\tilde{\rho }_{i}.
        \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1106}%
\begin{displaymath}
        S = \left\{
          \begin{array}{lcl} 
            +1 & \mbox{if} & \Delta _{i+\frac{1}{2}} \geq 0, \\
            -1 & \mbox{if} & \Delta f_{i+\frac{1}{2}} < 0.
          \end{array}
        \right. 
        \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1107}%
\begin{displaymath}
        \eta _{i+\frac{1}{2}}=\frac{1}{2}\epsilon _{i+\frac{1}{2}}
        (1-\epsilon _{i+\frac{1}{2}}),
        \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1108}%
\begin{displaymath}
        \epsilon _{i+\frac{1}{2}} = \frac{\Delta t}{\Delta x}
        \mbox{max}[u_{i}^{n}, u_{i+1}^{n}].
        \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1109}%
\begin{displaymath}
        \rho _{i}^{n+1}=\tilde{\rho }_{i} - (f_{i+\frac{1}{2}}^{c}-
        f_{i-\frac{1}{2}}^{c}).
        \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1117}%
$\eta _{i+\frac{1}{2}}|\Delta _
    {i+\frac{1}{2}}|$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

\stepcounter{subsection}
{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1131}%
$\Delta _
    {i-\frac{1}{2}}, \Delta _{i+\frac{1}{2}}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlfigureA{figure1063}%
\begin{figure}      \begin{center}
        \Depsf[][]{ps-fig/fig5.ps}
              \end{center}
    \end{figure}%
\lthtmlfigureZ
\lthtmlcheckvsize\clearpage}

\stepcounter{subsection}
{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1135}%
$(x,y)=(i\Delta x, j\Delta y)$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1137}%
$y$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1139}%
$x$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1141}%
$(i\Delta x, j\Delta y)$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlfigureA{figure1074}%
\begin{figure}      \begin{center}
        \Depsf[][]{ps-fig/fig6.ps}
              \end{center}
    \end{figure}%
\lthtmlfigureZ
\lthtmlcheckvsize\clearpage}

\stepcounter{section}
{\newpage\clearpage
\lthtmldisplayA{displaymath1193}%
\begin{displaymath}
      \rho _{i}^{n+1}=\rho _{i}^{n} - \frac{1}{\Delta x}[F_{i+\frac{1}{2}}
      -F_{i-\frac{1}{2}}],
    \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

\stepcounter{subsection}
{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1569}%
$F_{i+\frac{1}{2}}^{L}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1571}%
$F_{i+\frac{1}{2}}^{H}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1214}%
\begin{displaymath}
        A_{i+\frac{1}{2}}\equiv F_{i+\frac{1}{2}}^{H}-F_{i+\frac{1}{2}}^{L}.
      \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1225}%
\begin{displaymath}
        \rho _{i}^{td}=\rho _{i}^{n} - \frac{1}{\Delta x}[F_{i+\frac{1}{2}}^{L}
        -F_{i-\frac{1}{2}}^{L}].
      \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1240}%
\begin{displaymath}
        A_{i+\frac{1}{2}}^{c} = C_{i+\frac{1}{2}}A_{i+\frac{1}{2}}, \;
        0 \le C_{i+\frac{1}{2}}\le 1.
      \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1573}%
$C_{i+\frac{1}{2}}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1254}%
\begin{displaymath}
        \rho _{i}^{td}=\rho _{i}^{n} - \frac{1}{\Delta x}[A_{i+\frac{1}{2}}^{c}
        -A_{i-\frac{1}{2}}^{c}].
      \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1575}%
$F_{i+\frac{1}{2}}^{L}, F_{i+\frac{1}{2}}^{H}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1279}%
\begin{displaymath}
      \rho _{i}^{n+1}=\rho _{i}^{n}
       - \frac{1}{\Delta x}\left\{[C_{i+\frac{1}{2}}F_{i+\frac{1}{2}}^{H}+
      (1-C_{i+\frac{1}{2}})F_{i+\frac{1}{2}}^{L}]
      - [C_{i-\frac{1}{2}}F_{i-\frac{1}{2}}^{H}+
      (1-C_{i-\frac{1}{2}})F_{i-\frac{1}{2}}^{L}]\right\},
    \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

\stepcounter{subsection}
{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1585}%
$\rho _{i}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1587}%
$\rho _
    {i}^{max}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1589}%
$\rho _{i}^{min}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3007}%
$\displaystyle A_{i+\frac{1}{2}} =0$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3008}%
$\textstyle \mbox{if}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3009}%
$\displaystyle A_{i+\frac{1}{2}}(\rho _
{i+1}^{td} - \rho _{i}^{td}) < 0,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3010}%
$\textstyle \mbox{and either}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3011}%
$\displaystyle A_{i+\frac{1}{2}}(\rho _
{i+2}^{td} - \rho _{i+1}^{td}) < 0,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3013}%
$\textstyle \mbox{or}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3014}%
$\displaystyle A_{i+\frac{1}{2}}(\rho _
{i}^{td} - \rho _{i-1}^{td}) < 0.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1593}%
$ A_{i+\frac{1}{2}}(\rho _ {i+1}^{td} - \rho _
        {i}^{td}) < 0 $%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1595}%
$\rho _{i}^{max}, \rho _{i}^{min}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3019}%
$\displaystyle \rho _{i}^{max}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3021}%
$\displaystyle \mbox{max}(\rho _{i-1}^{td}, \; \rho _
{i}^{td}, \; \rho _{i+1}^{td}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3023}%
$\displaystyle \rho _{i}^{min}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3025}%
$\displaystyle \mbox{min}(\rho _{i-1}^{td}, \; \rho _
{i}^{td}, \; \rho _{i+1}^{td}).$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3028}%
$\displaystyle \rho _{i}^{a}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3030}%
$\displaystyle \mbox{max}(\rho _{i}^{n}, \; \rho _{i}^{td}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3034}%
$\displaystyle \mbox{max}(\rho _{i-1}^{a}, \; \rho _
{i}^{a}, \; \rho _{i+1}^{a}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3036}%
$\displaystyle \rho _{i}^{b}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3038}%
$\displaystyle \mbox{min}(\rho _{i}^{n}, \; \rho _{i}^{td}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3042}%
$\displaystyle \mbox{min}(\rho _{i-1}^{b}, \; \rho _
{i}^{b}, \; \rho _{i+1}^{b}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3044}%
$\displaystyle P_{i}^{+}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3046}%
$\displaystyle \mbox{格子点 $i$\  に入る全ての antidiffusive フ
ラックスの和}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3049}%
$\displaystyle \mbox{max}(0, A_{i-\frac{1}{2}}) - \mbox{min}(0,
A_{i+\frac{1}{2}}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3051}%
$\displaystyle Q_{i}^{+}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3053}%
$\displaystyle (\rho _{i}^{max}-\rho _{i})\Delta x,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3055}%
$\displaystyle R_{i}^{+}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3057}%
$\displaystyle \left\{
\begin{array}{lcl}
\mbox{min}(1, Q_{i}^{+}/P_{i}^{+}) & \mbox{if} & P_{i}^{+}>0, \\
0 & \mbox{if} & P_{i}^{+}=0.
\end{array}
\right.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1599}%
$R_{i}^{+}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3061}%
$\displaystyle P_{i}^{-}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3063}%
$\displaystyle \mbox{格子点 $i$\  から出て行く全ての antidiffusive フ
ラックスの和}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3066}%
$\displaystyle \mbox{max}(0, A_{i+\frac{1}{2}}) - \mbox{min}(0,
A_{i-\frac{1}{2}}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3068}%
$\displaystyle Q_{i}^{-}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3070}%
$\displaystyle (\rho _{i}-\rho _{i}^{min})\Delta x,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3072}%
$\displaystyle R_{i}^{-}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3074}%
$\displaystyle \left\{
\begin{array}{lcl}
\mbox{min}(1, Q_{i}^{-}/P_{i}^{-}) & \mbox{if} & P_{i}^{-}>0, \\
0 & \mbox{if} & P_{i}^{-}=0.
\end{array}
\right.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1605}%
$R_{i}^{-}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline1609}%
$R_{i}^{+}, R_{i}^{-}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{displaymath1502}%
\begin{displaymath}
          C_{i+\frac{1}{2}} = \left\{
            \begin{array}{lcl}
               \mbox{min}(R_{i+1}^{+}, R_{i}^{-}) & \mbox{if} &
              A_{i+\frac{1}{2}}\ge 0, \\
               \mbox{min}(R_{i}^{+}, R_{i+1}^{-}) & \mbox{if} &
              A_{i+\frac{1}{2}} < 0.
            \end{array}
            \right.
        \end{displaymath}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

\stepcounter{section}
\stepcounter{subsection}
{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3083}%
$\displaystyle \rho _{i.j}^{n+1}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3085}%
$\displaystyle \rho _{i.j}^{n} - \frac{1}{\Delta x\Delta
y}\{[C_{i+\frac{1}{2},j}F_{i+\frac{1}{2},j}^{H} +
(1-C_{i+\frac{1}{2},j})F_{i+\frac{1}{2},j}^{L}]$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3086}%
$\displaystyle - [C_{i-\frac{1}{2},j}F_{i-\frac{1}{2},j}^{H}+
(1-C_{i-\frac{1}{2},j})F_{i-\frac{1}{2},j}^{L}],$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3087}%
$\displaystyle + [C_{i,j+\frac{1}{2}}G_{i,j+\frac{1}{2}}^{H}+
(1-C_{i,j+\frac{1}{2}})G_{i,j+\frac{1}{2}}^{L}]$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3089}%
$\displaystyle - [C_{i,j-\frac{1}{2}}G_{i,j-\frac{1}{2}}^{H}+
(1-C_{i,j-\frac{1}{2}})G_{i,j-\frac{1}{2}}^{L}]\},$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2213}%
$C_{i\pm \frac{1}{2},j}, C_{i,j\pm \frac{1}{2}}$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

\setcounter{equation}{18}

\renewcommand{\theequation}{\arabic{equation}'}
{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3093}%
$\displaystyle A_{i+\frac{1}{2},j} =0$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3095}%
$\displaystyle A_{i+\frac{1}{2},j}(\rho _
{i+1,j}^{td} - \rho _{i,j}^{td}) < 0,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3097}%
$\displaystyle A_{i+\frac{1}{2},j}(\rho _
{i+2,j}^{td} - \rho _{i+1,j}^{td}) < 0,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3099}%
$\displaystyle A_{i+\frac{1}{2},j}(\rho _
{i,j}^{td} - \rho _{i-1,j}^{td}) < 0,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3100}%
$\displaystyle A_{i,j+\frac{1}{2}} =0$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3102}%
$\displaystyle A_{i,j+\frac{1}{2}}(\rho _
{i,j+1}^{td} - \rho _{i,j}^{td}) < 0,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3104}%
$\displaystyle A_{i,j+\frac{1}{2}}(\rho _
{i,j+2}^{td} - \rho _{i,j+1}^{td}) < 0,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3107}%
$\displaystyle A_{i,j+\frac{1}{2}}(\rho _
{i,j}^{td} - \rho _{i,j-1}^{td}) < 0.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3111}%
$\displaystyle \rho _{i,j}^{max}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3113}%
$\displaystyle \mbox{max}(\rho _{i-1,j}^{td}, \; \rho _
{i,j}^{td}, \; \rho _{i+1,j}^{td}, \; \rho _{i,j-1}^{td}, \;
\rho _{i,j+1}^{td}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3115}%
$\displaystyle \rho _{i,j}^{min}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3117}%
$\displaystyle \mbox{min}(\rho _{i-1,j}^{td}, \; \rho _
{i,j}^{td}, \; \rho _{i+1,j}^{td}, \; \rho _{i,j-1}^{td}, \;
\rho _{i,j+1}^{td}).$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3120}%
$\displaystyle \rho _{i,j}^{a}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3122}%
$\displaystyle \mbox{max}(\rho _{i,j}^{n}, \; \rho _
{i,j}^{td}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3126}%
$\displaystyle \mbox{max}(\rho _{i-1,j}^{a}, \; \rho _ {i,j}^{a}
, \; \rho _{i+1,j}^{a}, \; \rho _{i,j-1}^{a}, \;
\rho _{i,j+1}^{a}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3128}%
$\displaystyle \rho _{i,j}^{b}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3130}%
$\displaystyle \mbox{min}(\rho _{i,j}^{n}, \; \rho _
{i,j}^{td}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3134}%
$\displaystyle \mbox{min}(\rho _{i-1,j}^{b}, \;
\rho _ {i,j}^{b}, \; \rho _{i+1,j}^{b}, \; \rho _{i,j-1}^{b}, \;
\rho _{i,j-1}^{b}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3136}%
$\displaystyle P_{i,j}^{+}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3138}%
$\displaystyle \mbox{格子点 $i,j$\  に入る全ての antidiffusive フ
ラックスの和}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3140}%
$\displaystyle \mbox{max}(0, A_{i-\frac{1}{2},j}) - \mbox{min}(0,
A_{i+\frac{1}{2},j})$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3142}%
$\displaystyle + \mbox{max}(0, A_{i,j-\frac{1}{2}}) -
\mbox{min}(0, A_{i,j+\frac{1}{2}}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3144}%
$\displaystyle Q_{i,j}^{+}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3146}%
$\displaystyle (\rho _{i,j}^{max}-\rho _{i,j})\Delta x \Delta y,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3148}%
$\displaystyle R_{i,j}^{+}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3150}%
$\displaystyle \left\{
\begin{array}{lcl}
\mbox{min}(1, Q_{i.j}^{+}/P_{i,j}^{+}) & \mbox{if} & P_{i,j}^{+}>0,
\\
0 & \mbox{if} & P_{i,j}^{+}=0.
\end{array}
\right.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3152}%
$\displaystyle P_{i,j}^{-}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3154}%
$\displaystyle \mbox{格子点 $i,j$\  から出て行く全ての antidiffusive
フラックスの和}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3156}%
$\displaystyle \mbox{max}(0, A_{i+\frac{1}{2},j}) - \mbox{min}(0,
A_{i-\frac{1}{2},j})$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3158}%
$\displaystyle + \mbox{max}(0, A_{i,j+\frac{1}{2}}) -
\mbox{min}(0,A_{i,j-\frac{1}{2}}),$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3160}%
$\displaystyle Q_{i,j}^{-}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3162}%
$\displaystyle (\rho _{i,j}-\rho _{i,j}^{min})\Delta x \Delta y
,$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3164}%
$\displaystyle R_{i,j}^{-}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3166}%
$\displaystyle \left\{
\begin{array}{lcl}
\mbox{min}(1, Q_{i,j}^{-}/P_{i,j}^{-}) & \mbox{if} & P_{i,j}^{-}>0,
\\
0 & \mbox{if} & P_{i,j}^{-}=0.
\end{array}
\right.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3168}%
$\displaystyle C_{i+\frac{1}{2},j}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3170}%
$\displaystyle \left\{
\begin{array}{lcl}
\mbox{min}(R_{i+1,j}^{+}, R_{i,j}^{-}) & \mbox{if} &
A_{i+\frac{1}{2},j}\ge 0, \\
\mbox{min}(R_{i,j}^{+}, R_{i+1,j}^{-}) & \mbox{if} &
A_{i+\frac{1}{2},j} < 0.
\end{array}
\right.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3172}%
$\displaystyle C_{i,j+\frac{1}{2}}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3174}%
$\displaystyle \left\{
\begin{array}{lcl}
\mbox{min}(R_{i,j+1}^{+}, R_{i,j}^{-}) & \mbox{if} &
A_{i,j+\frac{1}{2}}\ge 0, \\
\mbox{min}(R_{i,j}^{+}, R_{i,j+1}^{-}) & \mbox{if} &
A_{i,j+\frac{1}{2}} < 0.
\end{array}
\right.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

\addtocounter{equation}{1}
\stepcounter{subsection}
{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2225}%
$0\leq x \leq 1, 0\leq y \leq 1$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2227}%
$\Delta x = 
        \Delta y = 0.01$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2229}%
$(x_{0}, y_{0})=(0.5, 0.5)$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2231}%
$\omega = 0.1$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2233}%
$(u, v) = ( - \omega (y-y_{0}), 
        \omega (x-x_{0}))$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2237}%
$(x_{m}, y_{m})=(0.75, 0.5)$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_inline2241}%
$(u\Delta t/\Delta x)$%
\lthtmlinlinemathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlfigureA{figure2006}%
\begin{figure}      \begin{center}
        \Depsf[][80mm]{ps-fig/init-cone.ps}
              \end{center}
    \end{figure}%
\lthtmlfigureZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3190}%
$\displaystyle \rho _{i,j}^{n+1}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3192}%
$\displaystyle \rho _{i,j}^{n} -
\{F_{x}(\rho _{i,j}^{n},\rho _{i+1,j}^{n},u_{i+\frac{1}{2},j}^{n}) -
F_{x}(\rho _{i-1,j}^{n},\rho _{i,j}^{n},u_{i-\frac{1}{2},j}^{n})\}$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlinlinemathA{tex2html_wrap_indisplay3194}%
$\displaystyle - \{F_{y}(\rho _{i,j}^{n},\rho _{i,j+1}^{n},v_{i,j+\frac{1}{2}}^{n})
- F_{y}(\rho _{i,j-1}^{n},\rho _{i,j}^{n},v_{i,j-\frac{1}{2}}^{n})\}.$%
\lthtmlindisplaymathZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmldisplayA{eqnarraystar2052}%
\begin{eqnarray*}
      F_{x}(\rho _{i,j}^{n},\rho _{i+1,j}^{n},u_{i+\frac{1}{2},j}^{n})&=& 
      [(u_{i+\frac{1}{2},j}^{n}+|u_{i+\frac{1}{2},j}^{n}|)\rho _{i,j}^{n} + 
      (u_{i+\frac{1}{2},j}^{n}-|u_{i+\frac{1}{2},j}^{n}|)\rho _{i+1,j}^{n}]
      \frac{\Delta t}{2\Delta x}, \\
      F_{y}(\rho _{i.j}^{n},\rho _{i,j+1}^{n},v_{i,j+\frac{1}{2}}^{n})&=&
      [(v_{i,j+\frac{1}{2}}^{n}+|v_{i,j+\frac{1}{2}}^{n}|)\rho _{i,j}^{n} + 
      (v_{i,j+\frac{1}{2}}^{n}-|v_{i,j+\frac{1}{2}}^{n}|)\rho _{i,j+1}^{n}]
      \frac{\Delta t}{2\Delta y}, 
    \end{eqnarray*}%
\lthtmldisplayZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlfigureA{figure2108}%
\begin{figure}      \begin{center}
        \Depsf[][80mm]{ps-fig/upstream1.ps}
                \Depsf[][80mm]{ps-fig/upstream2.ps}
              \end{center}
    \end{figure}%
\lthtmlfigureZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlfigureA{figure2120}%
\begin{figure}      \begin{center}
        \Depsf[][80mm]{ps-fig/fct-exp1.ps}
                \Depsf[][80mm]{ps-fig/fct-exp2.ps}
              \end{center}
    \end{figure}%
\lthtmlfigureZ
\lthtmlcheckvsize\clearpage}

{\newpage\clearpage
\lthtmlfigureA{figure2133}%
\begin{figure}      \begin{center}
        \Depsf[][80mm]{ps-fig/init-cylind.ps}
                \Depsf[][80mm]{ps-fig/fct-exp3.ps}
              \end{center}
    \end{figure}%
\lthtmlfigureZ
\lthtmlcheckvsize\clearpage}

\stepcounter{section}

\end{document}
