The Use of Mathematica for Computation of Stability Regions with Guaranteed Accuracy
Date: July 18th (Thursday)
Time: 10:20-10:45
Abstract
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.