> y = sin ( mp (x ) ) % compute sine values in quadruple precision > f = sin (x ) % call built-in function for double precision sine % Check accuracy of sine using 1M random points on (0, 16] Having extended precision routines, accuracy check of any function needs only few commands: Along the way, we will also check accuracy of the commonly used open source libraries. Today we will investigate how accurately MATLAB computes Bessel functions of the first and second kind Y n(x) and J n(x) in double precision. Please refer to the pages for more details. In spite of the fact that modified Bessel functions are easy to compute (they are monotonous and do not cross x-axis) we saw that MATLAB provides accuracy much lower than expected for double precision. In previous posts we studied accuracy of computation of modified Bessel functions: K 1(x), K 0(x), I 0(x) and I 1(x).
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |