\begin{center}
\LARGE
he Use of Mathematica for Computation of Stability
              Regions with Guaranteed Accuracy

\end{center}

\begin{center}
\Large
Victor G. Ganzha, K. Schaub, Evgenii V. Vorozhtsov

\end{center}

\noindent
Date:  July 18th (Thursday)\\
Time:  10:20-10:45\\

\begin{center}
\large
Abstract
\end{center}

We present a symbolic-numerical algorithm for the stability investigation of difference schemes
approximating the partial differential equations of hyperbolic or parabolic type. The use of a
combination of rational arithmetic and the arithmetic of floating-point numbers available in the
Mathematica system has enabled us to compute the stability region boundaries with guaranteed
accuracy. Computational examples are given, in which we analyze a number of difference schemes
for one- and two-dimensional partial differential equations. 

