############################################################
# Output from 'catprec' illustrating catastrophic precision 
# loss due to subtraction of nearly-equal floating point 
# values
############################################################

% catprec

       h      d(sin) approx   d(sin) exact    d(sin) err

   5.000E-01  8.91096713E-01  9.80066578E-01  -8.897E-02
   6.250E-02  9.73222242E-01  9.80066578E-01  -6.844E-03
   7.812E-03  9.79280560E-01  9.80066578E-01  -7.860E-04
   9.766E-04  9.79969416E-01  9.80066578E-01  -9.716E-05
   1.221E-04  9.80054450E-01  9.80066578E-01  -1.213E-05
   1.526E-05  9.80065062E-01  9.80066578E-01  -1.516E-06
   1.907E-06  9.80066388E-01  9.80066578E-01  -1.895E-07
   2.384E-07  9.80066554E-01  9.80066578E-01  -2.369E-08
   2.980E-08  9.80066574E-01  9.80066578E-01  -3.663E-09
   3.725E-09  9.80066575E-01  9.80066578E-01  -2.731E-09
   4.657E-10  9.80066538E-01  9.80066578E-01  -3.998E-08
   5.821E-11  9.80066299E-01  9.80066578E-01  -2.784E-07
   7.276E-12  9.80064392E-01  9.80066578E-01  -2.186E-06
   9.095E-13  9.80041504E-01  9.80066578E-01  -2.507E-05
   1.137E-13  9.79980469E-01  9.80066578E-01  -8.611E-05
   1.421E-14  9.78515625E-01  9.80066578E-01  -1.551E-03