vpa_test

vpa_test, a MATLAB code which calls vpa(), which implements MATLAB's variable precision arithmetic (VPA) feature.

MATLAB's variable precision arithmetic is only available through the Symbolic Math Toolbox. To find out if you have this toolbox installed, you can type "ver" within your MATLAB session.

In MATLAB, evaluating a symbolic variable means creating a representation of its value, to a given number of digits, as a floating point number. The vpa() command is used to create the floating point representation. Typical uses are:

        r = vpa ( s )
      

        r = vpa ( s, d )
      

The current default number of digits can be found by using the digits() command:

        dignum = digits;
      

        dignum = digits ( );
      

        digits ( 50 );
      

To my mind, a significant issue with the symbolic math toolbox is that it is not apparent how to control the printing of the value of a given variable. If I compute something to 80 digits, I can always display it by naming it. But what if I prefer a labeled printout? Various alternatives involving "disp()" and "fprintf()" and string concatenation seem unsatisfactory.

Licensing:

The computer code and data files made available on this web page are distributed under the BSD license.

Languages:

vpa_test is available in a MATLAB version and an Octave version.

Reference:

1. David Bailey, Jonathan Borwein,
   Experimental Mathematics: Examples, Methods and Implications,
   Notices of the American Mathematical Society,
   Volume 52, Number 5, pages 502-514.
2. David Bailey, Karthik Jeyabalan, Xiaoye Li,
   A Comparison of Three High-Precision Quadrature Schemes,
   Experimental Mathematics,
   Volume 14, Number 3, pages 317-329.
3. Cleve Moler,
   Cleve's Corner: Computing Pi,
   MATLAB News and Notes,
   February 2012.

Source Code:

• vpa_test.m, calls all the tests.
• vpa_test.sh, runs all the tests.
• vpa_test.txt, the output file.

computing_pi is a collection of three functions discussed by Cleve Moler for computing pi to many digits.

• agm_pi.m, uses the arithmetic-geometric mean algorithm.
• agm_pi_test.m, tests agm_pi.
• bbp_pi.m, uses the Bailey-Borwein-Plouffe algorithm.
• bbp_pi_test.m, tests bbp_pi.
• chud_pi.m, uses the Chudnovsky algorithm.
• chud_pi_test.m, tests chud_pi.