3709
 3710
 3711
 3712
 3713
 3714
 3715
 3716      DOUBLE PRECISION   ERRBND, PREC, PUSCLR
 3717      COMPLEX*16         X
 3718
 3719
 3720
 3721
 3722
 3723
 3724
 3725
 3726
 3727
 3728
 3729
 3730
 3731
 3732
 3733
 3734
 3735
 3736
 3737
 3738
 3739
 3740
 3741
 3742
 3743
 3744
 3745
 3746
 3747
 3748
 3749
 3750
 3751
 3752
 3753
 3754
 3755
 3756
 3757
 3758
 3759      DOUBLE PRECISION   TWO
 3760      parameter( two = 2.0d+0 )
 3761
 3762
 3763      INTRINSIC          abs, dble, dcmplx, dimag
 3764
 3765
 3766
 3767      x = dcmplx( pusclr * dble( x ), pusclr * dimag( x ) )
 3768
 3769      errbnd = ( two * prec ) * abs( x )
 3770
 3771      RETURN
 3772
 3773
 3774