CRAN Package Check Results for Package DPQmpfr

Last updated on 2021-04-20 23:47:40 CEST.

Flavor Version Tinstall Tcheck Ttotal Status Flags
r-devel-linux-x86_64-debian-clang 0.3-0 4.33 167.94 172.27 OK
r-devel-linux-x86_64-debian-gcc 0.3-0 3.97 127.19 131.16 OK
r-devel-linux-x86_64-fedora-clang 0.3-0 209.19 NOTE
r-devel-linux-x86_64-fedora-gcc 0.3-0 200.11 NOTE
r-devel-windows-ix86+x86_64 0.3-0 16.00 168.00 184.00 OK
r-devel-windows-x86_64-gcc10-UCRT 0.3-0 NOTE
r-patched-linux-x86_64 0.3-0 6.45 159.18 165.63 WARN
r-patched-solaris-x86 0.3-0 229.50 NOTE
r-release-linux-x86_64 0.3-0 4.33 159.51 163.84 OK
r-release-macos-x86_64 0.3-0 NOTE
r-release-windows-ix86+x86_64 0.3-0 11.00 217.00 228.00 OK
r-oldrel-macos-x86_64 0.3-0 ERROR
r-oldrel-windows-ix86+x86_64 0.3-0 8.00 189.00 197.00 ERROR

Check Details

Version: 0.3-0
Check: dependencies in R code
Result: NOTE
    Namespace in Imports field not imported from: ‘graphics’
     All declared Imports should be used.
Flavors: r-devel-linux-x86_64-fedora-clang, r-devel-linux-x86_64-fedora-gcc, r-devel-windows-x86_64-gcc10-UCRT, r-patched-solaris-x86, r-release-macos-x86_64

Version: 0.3-0
Check: whether package can be installed
Result: WARN
    Found the following significant warnings:
     Rd warning: /tmp/RtmpQxImo6/R.INSTALL1b6d64b658e1/DPQmpfr/man/hyperGeomQ.Rd:9: file link ‘chooseZ’ in package ‘gmp’ does not exist and so has been treated as a topic
     Rd warning: /tmp/RtmpQxImo6/R.INSTALL1b6d64b658e1/DPQmpfr/man/hyperGeomQ.Rd:35: file link ‘chooseZ’ in package ‘gmp’ does not exist and so has been treated as a topic
Flavor: r-patched-linux-x86_64

Version: 0.3-0
Check: whether package can be installed
Result: ERROR
    Installation failed.
Flavor: r-oldrel-macos-x86_64

Version: 0.3-0
Check: tests
Result: ERROR
     Running 'beta-Ding94.R' [107s]
     Running 'dpqnorm-tst.R' [3s]
     Running 'hypergeom.R' [3s]
     Running 't-nonc_tst.R' [6s]
    Running the tests in 'tests/dpqnorm-tst.R' failed.
    Complete output:
     > require(DPQmpfr)
     Loading required package: DPQmpfr
     > require(Rmpfr)
     Loading required package: Rmpfr
     Loading required package: gmp
    
     Attaching package: 'gmp'
    
     The following objects are masked from 'package:base':
    
     %*%, apply, crossprod, matrix, tcrossprod
    
     C code of R package 'Rmpfr': GMP using 64 bits per limb
    
    
     Attaching package: 'Rmpfr'
    
     The following object is masked from 'package:gmp':
    
     outer
    
     The following objects are masked from 'package:stats':
    
     dbinom, dgamma, dnbinom, dnorm, dpois, pnorm
    
     The following objects are masked from 'package:base':
    
     cbind, pmax, pmin, rbind
    
     >
     > ## very accurate set of *large* abscissa values
     > x <- mpfr(2,1024)^ seq(3, 100, by=1/2)
     > xN <- asNumeric(x)
     >
     > ## Looking at pnormU_S53() = upper bound to pnorm(), for large x,
     > ## the term log(4 / (3 + sqrt(1 + 8/x^2))) early underflows to 0:
     > eps_S53 <- function(x, method = c("direct", "log1pExct", "deltaFull",
     + "delta1", "delta2", "delta3")) {
     + method <- match.arg(method)
     + delta <- function(x) { xm2 <- x^-2; (.5*(1 - sqrt(1 + 8*xm2)) + xm2) / (1 + xm2) }
     + switch(method,
     + "direct" = log(4 / (3 + sqrt(1 + 8*x^-2))),
     + "log1pExct" = log1p(delta(x)),
     + "deltaFull" = delta(x),
     + "delta1" = -1/(x^2+1) , ## 1 term for delta(x)
     + "delta2" = -(1 - (2/x)^2) / (x^2+1) , ## 2 terms for delta(x)
     + "delta3" = {tx2 <- (2/x)^2 ; -(1 - tx2*(1 - tx2)) / (x^2+1) }, ## 3 terms
     + stop("invalid 'method': ", method))
     + }
     >
     >
     > (S53meth <- eval(formals(eps_S53)$method))
     [1] "direct" "log1pExct" "deltaFull" "delta1" "delta2" "delta3"
     > ## "direct" "log1pExct" "deltaFull" "delta1" "delta2" "delta3"
     >
     > epsAllS53 <- function(x, meth = eval(formals(eps_S53)$method)) {
     + stopifnot(is.character(meth), is.numeric(x) || is(x,"mpfr"))
     + Rl <- lapply(meth, function(m) eps_S53(x, method=m))
     + new("mpfrMatrix", unlist(Rl),
     + Dim = c(length(x), length(meth)),
     + Dimnames = list(NULL, meth))
     + }
     >
     > EA <- epsAllS53(x)
     > asNumeric(EA) # looks +- ok
     direct log1pExct deltaFull delta1 delta2
     [1,] -1.505120e-02 -1.458470e-02 -1.447885e-02 -1.538462e-02 -1.442308e-02
     [2,] -7.664653e-03 -7.545372e-03 -7.516977e-03 -7.751938e-03 -7.509690e-03
     [3,] -3.868707e-03 -3.838543e-03 -3.831185e-03 -3.891051e-03 -3.830253e-03
     [4,] -1.943664e-03 -1.936080e-03 -1.934207e-03 -1.949318e-03 -1.934089e-03
     [5,] -9.741879e-04 -9.722861e-04 -9.718136e-04 -9.756098e-04 -9.717988e-04
     [6,] -4.876864e-04 -4.872103e-04 -4.870916e-04 -4.880429e-04 -4.870897e-04
     [7,] -2.439918e-04 -2.438726e-04 -2.438429e-04 -2.440810e-04 -2.438427e-04
     [8,] -1.220331e-04 -1.220033e-04 -1.219958e-04 -1.220554e-04 -1.219958e-04
     [9,] -6.102585e-05 -6.101840e-05 -6.101653e-05 -6.103143e-05 -6.101653e-05
     [10,] -3.051525e-05 -3.051339e-05 -3.051292e-05 -3.051665e-05 -3.051292e-05
     [11,] -1.525821e-05 -1.525774e-05 -1.525762e-05 -1.525856e-05 -1.525762e-05
     [12,] -7.629249e-06 -7.629133e-06 -7.629104e-06 -7.629336e-06 -7.629103e-06
     [13,] -3.814661e-06 -3.814632e-06 -3.814625e-06 -3.814683e-06 -3.814625e-06
     [14,] -1.907340e-06 -1.907332e-06 -1.907330e-06 -1.907345e-06 -1.907330e-06
     [15,] -9.536720e-07 -9.536702e-07 -9.536698e-07 -9.536734e-07 -9.536698e-07
     [16,] -4.768366e-07 -4.768361e-07 -4.768360e-07 -4.768369e-07 -4.768360e-07
     [17,] -2.384184e-07 -2.384183e-07 -2.384183e-07 -2.384185e-07 -2.384183e-07
     [18,] -1.192093e-07 -1.192092e-07 -1.192092e-07 -1.192093e-07 -1.192092e-07
     [19,] -5.960464e-08 -5.960463e-08 -5.960463e-08 -5.960464e-08 -5.960463e-08
     [20,] -2.980232e-08 -2.980232e-08 -2.980232e-08 -2.980232e-08 -2.980232e-08
     [21,] -1.490116e-08 -1.490116e-08 -1.490116e-08 -1.490116e-08 -1.490116e-08
     [22,] -7.450580e-09 -7.450580e-09 -7.450580e-09 -7.450581e-09 -7.450580e-09
     [23,] -3.725290e-09 -3.725290e-09 -3.725290e-09 -3.725290e-09 -3.725290e-09
     [24,] -1.862645e-09 -1.862645e-09 -1.862645e-09 -1.862645e-09 -1.862645e-09
     [25,] -9.313226e-10 -9.313226e-10 -9.313226e-10 -9.313226e-10 -9.313226e-10
     [26,] -4.656613e-10 -4.656613e-10 -4.656613e-10 -4.656613e-10 -4.656613e-10
     [27,] -2.328306e-10 -2.328306e-10 -2.328306e-10 -2.328306e-10 -2.328306e-10
     [28,] -1.164153e-10 -1.164153e-10 -1.164153e-10 -1.164153e-10 -1.164153e-10
     [29,] -5.820766e-11 -5.820766e-11 -5.820766e-11 -5.820766e-11 -5.820766e-11
     [30,] -2.910383e-11 -2.910383e-11 -2.910383e-11 -2.910383e-11 -2.910383e-11
     [31,] -1.455192e-11 -1.455192e-11 -1.455192e-11 -1.455192e-11 -1.455192e-11
     [32,] -7.275958e-12 -7.275958e-12 -7.275958e-12 -7.275958e-12 -7.275958e-12
     [33,] -3.637979e-12 -3.637979e-12 -3.637979e-12 -3.637979e-12 -3.637979e-12
     [34,] -1.818989e-12 -1.818989e-12 -1.818989e-12 -1.818989e-12 -1.818989e-12
     [35,] -9.094947e-13 -9.094947e-13 -9.094947e-13 -9.094947e-13 -9.094947e-13
     [36,] -4.547474e-13 -4.547474e-13 -4.547474e-13 -4.547474e-13 -4.547474e-13
     [37,] -2.273737e-13 -2.273737e-13 -2.273737e-13 -2.273737e-13 -2.273737e-13
     [38,] -1.136868e-13 -1.136868e-13 -1.136868e-13 -1.136868e-13 -1.136868e-13
     [39,] -5.684342e-14 -5.684342e-14 -5.684342e-14 -5.684342e-14 -5.684342e-14
     [40,] -2.842171e-14 -2.842171e-14 -2.842171e-14 -2.842171e-14 -2.842171e-14
     [41,] -1.421085e-14 -1.421085e-14 -1.421085e-14 -1.421085e-14 -1.421085e-14
     [42,] -7.105427e-15 -7.105427e-15 -7.105427e-15 -7.105427e-15 -7.105427e-15
     [43,] -3.552714e-15 -3.552714e-15 -3.552714e-15 -3.552714e-15 -3.552714e-15
     [44,] -1.776357e-15 -1.776357e-15 -1.776357e-15 -1.776357e-15 -1.776357e-15
     [45,] -8.881784e-16 -8.881784e-16 -8.881784e-16 -8.881784e-16 -8.881784e-16
     [46,] -4.440892e-16 -4.440892e-16 -4.440892e-16 -4.440892e-16 -4.440892e-16
     [47,] -2.220446e-16 -2.220446e-16 -2.220446e-16 -2.220446e-16 -2.220446e-16
     [48,] -1.110223e-16 -1.110223e-16 -1.110223e-16 -1.110223e-16 -1.110223e-16
     [49,] -5.551115e-17 -5.551115e-17 -5.551115e-17 -5.551115e-17 -5.551115e-17
     [50,] -2.775558e-17 -2.775558e-17 -2.775558e-17 -2.775558e-17 -2.775558e-17
     [51,] -1.387779e-17 -1.387779e-17 -1.387779e-17 -1.387779e-17 -1.387779e-17
     [52,] -6.938894e-18 -6.938894e-18 -6.938894e-18 -6.938894e-18 -6.938894e-18
     [53,] -3.469447e-18 -3.469447e-18 -3.469447e-18 -3.469447e-18 -3.469447e-18
     [54,] -1.734723e-18 -1.734723e-18 -1.734723e-18 -1.734723e-18 -1.734723e-18
     [55,] -8.673617e-19 -8.673617e-19 -8.673617e-19 -8.673617e-19 -8.673617e-19
     [56,] -4.336809e-19 -4.336809e-19 -4.336809e-19 -4.336809e-19 -4.336809e-19
     [57,] -2.168404e-19 -2.168404e-19 -2.168404e-19 -2.168404e-19 -2.168404e-19
     [58,] -1.084202e-19 -1.084202e-19 -1.084202e-19 -1.084202e-19 -1.084202e-19
     [59,] -5.421011e-20 -5.421011e-20 -5.421011e-20 -5.421011e-20 -5.421011e-20
     [60,] -2.710505e-20 -2.710505e-20 -2.710505e-20 -2.710505e-20 -2.710505e-20
     [61,] -1.355253e-20 -1.355253e-20 -1.355253e-20 -1.355253e-20 -1.355253e-20
     [62,] -6.776264e-21 -6.776264e-21 -6.776264e-21 -6.776264e-21 -6.776264e-21
     [63,] -3.388132e-21 -3.388132e-21 -3.388132e-21 -3.388132e-21 -3.388132e-21
     [64,] -1.694066e-21 -1.694066e-21 -1.694066e-21 -1.694066e-21 -1.694066e-21
     [65,] -8.470329e-22 -8.470329e-22 -8.470329e-22 -8.470329e-22 -8.470329e-22
     [66,] -4.235165e-22 -4.235165e-22 -4.235165e-22 -4.235165e-22 -4.235165e-22
     [67,] -2.117582e-22 -2.117582e-22 -2.117582e-22 -2.117582e-22 -2.117582e-22
     [68,] -1.058791e-22 -1.058791e-22 -1.058791e-22 -1.058791e-22 -1.058791e-22
     [69,] -5.293956e-23 -5.293956e-23 -5.293956e-23 -5.293956e-23 -5.293956e-23
     [70,] -2.646978e-23 -2.646978e-23 -2.646978e-23 -2.646978e-23 -2.646978e-23
     [71,] -1.323489e-23 -1.323489e-23 -1.323489e-23 -1.323489e-23 -1.323489e-23
     [72,] -6.617445e-24 -6.617445e-24 -6.617445e-24 -6.617445e-24 -6.617445e-24
     [73,] -3.308722e-24 -3.308722e-24 -3.308722e-24 -3.308722e-24 -3.308722e-24
     [74,] -1.654361e-24 -1.654361e-24 -1.654361e-24 -1.654361e-24 -1.654361e-24
     [75,] -8.271806e-25 -8.271806e-25 -8.271806e-25 -8.271806e-25 -8.271806e-25
     [76,] -4.135903e-25 -4.135903e-25 -4.135903e-25 -4.135903e-25 -4.135903e-25
     [77,] -2.067952e-25 -2.067952e-25 -2.067952e-25 -2.067952e-25 -2.067952e-25
     [78,] -1.033976e-25 -1.033976e-25 -1.033976e-25 -1.033976e-25 -1.033976e-25
     [79,] -5.169879e-26 -5.169879e-26 -5.169879e-26 -5.169879e-26 -5.169879e-26
     [80,] -2.584939e-26 -2.584939e-26 -2.584939e-26 -2.584939e-26 -2.584939e-26
     [81,] -1.292470e-26 -1.292470e-26 -1.292470e-26 -1.292470e-26 -1.292470e-26
     [82,] -6.462349e-27 -6.462349e-27 -6.462349e-27 -6.462349e-27 -6.462349e-27
     [83,] -3.231174e-27 -3.231174e-27 -3.231174e-27 -3.231174e-27 -3.231174e-27
     [84,] -1.615587e-27 -1.615587e-27 -1.615587e-27 -1.615587e-27 -1.615587e-27
     [85,] -8.077936e-28 -8.077936e-28 -8.077936e-28 -8.077936e-28 -8.077936e-28
     [86,] -4.038968e-28 -4.038968e-28 -4.038968e-28 -4.038968e-28 -4.038968e-28
     [87,] -2.019484e-28 -2.019484e-28 -2.019484e-28 -2.019484e-28 -2.019484e-28
     [88,] -1.009742e-28 -1.009742e-28 -1.009742e-28 -1.009742e-28 -1.009742e-28
     [89,] -5.048710e-29 -5.048710e-29 -5.048710e-29 -5.048710e-29 -5.048710e-29
     [90,] -2.524355e-29 -2.524355e-29 -2.524355e-29 -2.524355e-29 -2.524355e-29
     [91,] -1.262177e-29 -1.262177e-29 -1.262177e-29 -1.262177e-29 -1.262177e-29
     [92,] -6.310887e-30 -6.310887e-30 -6.310887e-30 -6.310887e-30 -6.310887e-30
     [93,] -3.155444e-30 -3.155444e-30 -3.155444e-30 -3.155444e-30 -3.155444e-30
     [94,] -1.577722e-30 -1.577722e-30 -1.577722e-30 -1.577722e-30 -1.577722e-30
     [95,] -7.888609e-31 -7.888609e-31 -7.888609e-31 -7.888609e-31 -7.888609e-31
     [96,] -3.944305e-31 -3.944305e-31 -3.944305e-31 -3.944305e-31 -3.944305e-31
     [97,] -1.972152e-31 -1.972152e-31 -1.972152e-31 -1.972152e-31 -1.972152e-31
     [98,] -9.860761e-32 -9.860761e-32 -9.860761e-32 -9.860761e-32 -9.860761e-32
     [99,] -4.930381e-32 -4.930381e-32 -4.930381e-32 -4.930381e-32 -4.930381e-32
     [100,] -2.465190e-32 -2.465190e-32 -2.465190e-32 -2.465190e-32 -2.465190e-32
     [101,] -1.232595e-32 -1.232595e-32 -1.232595e-32 -1.232595e-32 -1.232595e-32
     [102,] -6.162976e-33 -6.162976e-33 -6.162976e-33 -6.162976e-33 -6.162976e-33
     [103,] -3.081488e-33 -3.081488e-33 -3.081488e-33 -3.081488e-33 -3.081488e-33
     [104,] -1.540744e-33 -1.540744e-33 -1.540744e-33 -1.540744e-33 -1.540744e-33
     [105,] -7.703720e-34 -7.703720e-34 -7.703720e-34 -7.703720e-34 -7.703720e-34
     [106,] -3.851860e-34 -3.851860e-34 -3.851860e-34 -3.851860e-34 -3.851860e-34
     [107,] -1.925930e-34 -1.925930e-34 -1.925930e-34 -1.925930e-34 -1.925930e-34
     [108,] -9.629650e-35 -9.629650e-35 -9.629650e-35 -9.629650e-35 -9.629650e-35
     [109,] -4.814825e-35 -4.814825e-35 -4.814825e-35 -4.814825e-35 -4.814825e-35
     [110,] -2.407412e-35 -2.407412e-35 -2.407412e-35 -2.407412e-35 -2.407412e-35
     [111,] -1.203706e-35 -1.203706e-35 -1.203706e-35 -1.203706e-35 -1.203706e-35
     [112,] -6.018531e-36 -6.018531e-36 -6.018531e-36 -6.018531e-36 -6.018531e-36
     [113,] -3.009266e-36 -3.009266e-36 -3.009266e-36 -3.009266e-36 -3.009266e-36
     [114,] -1.504633e-36 -1.504633e-36 -1.504633e-36 -1.504633e-36 -1.504633e-36
     [115,] -7.523164e-37 -7.523164e-37 -7.523164e-37 -7.523164e-37 -7.523164e-37
     [116,] -3.761582e-37 -3.761582e-37 -3.761582e-37 -3.761582e-37 -3.761582e-37
     [117,] -1.880791e-37 -1.880791e-37 -1.880791e-37 -1.880791e-37 -1.880791e-37
     [118,] -9.403955e-38 -9.403955e-38 -9.403955e-38 -9.403955e-38 -9.403955e-38
     [119,] -4.701977e-38 -4.701977e-38 -4.701977e-38 -4.701977e-38 -4.701977e-38
     [120,] -2.350989e-38 -2.350989e-38 -2.350989e-38 -2.350989e-38 -2.350989e-38
     [121,] -1.175494e-38 -1.175494e-38 -1.175494e-38 -1.175494e-38 -1.175494e-38
     [122,] -5.877472e-39 -5.877472e-39 -5.877472e-39 -5.877472e-39 -5.877472e-39
     [123,] -2.938736e-39 -2.938736e-39 -2.938736e-39 -2.938736e-39 -2.938736e-39
     [124,] -1.469368e-39 -1.469368e-39 -1.469368e-39 -1.469368e-39 -1.469368e-39
     [125,] -7.346840e-40 -7.346840e-40 -7.346840e-40 -7.346840e-40 -7.346840e-40
     [126,] -3.673420e-40 -3.673420e-40 -3.673420e-40 -3.673420e-40 -3.673420e-40
     [127,] -1.836710e-40 -1.836710e-40 -1.836710e-40 -1.836710e-40 -1.836710e-40
     [128,] -9.183550e-41 -9.183550e-41 -9.183550e-41 -9.183550e-41 -9.183550e-41
     [129,] -4.591775e-41 -4.591775e-41 -4.591775e-41 -4.591775e-41 -4.591775e-41
     [130,] -2.295887e-41 -2.295887e-41 -2.295887e-41 -2.295887e-41 -2.295887e-41
     [131,] -1.147944e-41 -1.147944e-41 -1.147944e-41 -1.147944e-41 -1.147944e-41
     [132,] -5.739719e-42 -5.739719e-42 -5.739719e-42 -5.739719e-42 -5.739719e-42
     [133,] -2.869859e-42 -2.869859e-42 -2.869859e-42 -2.869859e-42 -2.869859e-42
     [134,] -1.434930e-42 -1.434930e-42 -1.434930e-42 -1.434930e-42 -1.434930e-42
     [135,] -7.174648e-43 -7.174648e-43 -7.174648e-43 -7.174648e-43 -7.174648e-43
     [136,] -3.587324e-43 -3.587324e-43 -3.587324e-43 -3.587324e-43 -3.587324e-43
     [137,] -1.793662e-43 -1.793662e-43 -1.793662e-43 -1.793662e-43 -1.793662e-43
     [138,] -8.968310e-44 -8.968310e-44 -8.968310e-44 -8.968310e-44 -8.968310e-44
     [139,] -4.484155e-44 -4.484155e-44 -4.484155e-44 -4.484155e-44 -4.484155e-44
     [140,] -2.242078e-44 -2.242078e-44 -2.242078e-44 -2.242078e-44 -2.242078e-44
     [141,] -1.121039e-44 -1.121039e-44 -1.121039e-44 -1.121039e-44 -1.121039e-44
     [142,] -5.605194e-45 -5.605194e-45 -5.605194e-45 -5.605194e-45 -5.605194e-45
     [143,] -2.802597e-45 -2.802597e-45 -2.802597e-45 -2.802597e-45 -2.802597e-45
     [144,] -1.401298e-45 -1.401298e-45 -1.401298e-45 -1.401298e-45 -1.401298e-45
     [145,] -7.006492e-46 -7.006492e-46 -7.006492e-46 -7.006492e-46 -7.006492e-46
     [146,] -3.503246e-46 -3.503246e-46 -3.503246e-46 -3.503246e-46 -3.503246e-46
     [147,] -1.751623e-46 -1.751623e-46 -1.751623e-46 -1.751623e-46 -1.751623e-46
     [148,] -8.758115e-47 -8.758115e-47 -8.758115e-47 -8.758115e-47 -8.758115e-47
     [149,] -4.379058e-47 -4.379058e-47 -4.379058e-47 -4.379058e-47 -4.379058e-47
     [150,] -2.189529e-47 -2.189529e-47 -2.189529e-47 -2.189529e-47 -2.189529e-47
     [151,] -1.094764e-47 -1.094764e-47 -1.094764e-47 -1.094764e-47 -1.094764e-47
     [152,] -5.473822e-48 -5.473822e-48 -5.473822e-48 -5.473822e-48 -5.473822e-48
     [153,] -2.736911e-48 -2.736911e-48 -2.736911e-48 -2.736911e-48 -2.736911e-48
     [154,] -1.368456e-48 -1.368456e-48 -1.368456e-48 -1.368456e-48 -1.368456e-48
     [155,] -6.842278e-49 -6.842278e-49 -6.842278e-49 -6.842278e-49 -6.842278e-49
     [156,] -3.421139e-49 -3.421139e-49 -3.421139e-49 -3.421139e-49 -3.421139e-49
     [157,] -1.710569e-49 -1.710569e-49 -1.710569e-49 -1.710569e-49 -1.710569e-49
     [158,] -8.552847e-50 -8.552847e-50 -8.552847e-50 -8.552847e-50 -8.552847e-50
     [159,] -4.276424e-50 -4.276424e-50 -4.276424e-50 -4.276424e-50 -4.276424e-50
     [160,] -2.138212e-50 -2.138212e-50 -2.138212e-50 -2.138212e-50 -2.138212e-50
     [161,] -1.069106e-50 -1.069106e-50 -1.069106e-50 -1.069106e-50 -1.069106e-50
     [162,] -5.345529e-51 -5.345529e-51 -5.345529e-51 -5.345529e-51 -5.345529e-51
     [163,] -2.672765e-51 -2.672765e-51 -2.672765e-51 -2.672765e-51 -2.672765e-51
     [164,] -1.336382e-51 -1.336382e-51 -1.336382e-51 -1.336382e-51 -1.336382e-51
     [165,] -6.681912e-52 -6.681912e-52 -6.681912e-52 -6.681912e-52 -6.681912e-52
     [166,] -3.340956e-52 -3.340956e-52 -3.340956e-52 -3.340956e-52 -3.340956e-52
     [167,] -1.670478e-52 -1.670478e-52 -1.670478e-52 -1.670478e-52 -1.670478e-52
     [168,] -8.352390e-53 -8.352390e-53 -8.352390e-53 -8.352390e-53 -8.352390e-53
     [169,] -4.176195e-53 -4.176195e-53 -4.176195e-53 -4.176195e-53 -4.176195e-53
     [170,] -2.088097e-53 -2.088097e-53 -2.088097e-53 -2.088097e-53 -2.088097e-53
     [171,] -1.044049e-53 -1.044049e-53 -1.044049e-53 -1.044049e-53 -1.044049e-53
     [172,] -5.220244e-54 -5.220244e-54 -5.220244e-54 -5.220244e-54 -5.220244e-54
     [173,] -2.610122e-54 -2.610122e-54 -2.610122e-54 -2.610122e-54 -2.610122e-54
     [174,] -1.305061e-54 -1.305061e-54 -1.305061e-54 -1.305061e-54 -1.305061e-54
     [175,] -6.525304e-55 -6.525304e-55 -6.525304e-55 -6.525304e-55 -6.525304e-55
     [176,] -3.262652e-55 -3.262652e-55 -3.262652e-55 -3.262652e-55 -3.262652e-55
     [177,] -1.631326e-55 -1.631326e-55 -1.631326e-55 -1.631326e-55 -1.631326e-55
     [178,] -8.156631e-56 -8.156631e-56 -8.156631e-56 -8.156631e-56 -8.156631e-56
     [179,] -4.078315e-56 -4.078315e-56 -4.078315e-56 -4.078315e-56 -4.078315e-56
     [180,] -2.039158e-56 -2.039158e-56 -2.039158e-56 -2.039158e-56 -2.039158e-56
     [181,] -1.019579e-56 -1.019579e-56 -1.019579e-56 -1.019579e-56 -1.019579e-56
     [182,] -5.097894e-57 -5.097894e-57 -5.097894e-57 -5.097894e-57 -5.097894e-57
     [183,] -2.548947e-57 -2.548947e-57 -2.548947e-57 -2.548947e-57 -2.548947e-57
     [184,] -1.274474e-57 -1.274474e-57 -1.274474e-57 -1.274474e-57 -1.274474e-57
     [185,] -6.372368e-58 -6.372368e-58 -6.372368e-58 -6.372368e-58 -6.372368e-58
     [186,] -3.186184e-58 -3.186184e-58 -3.186184e-58 -3.186184e-58 -3.186184e-58
     [187,] -1.593092e-58 -1.593092e-58 -1.593092e-58 -1.593092e-58 -1.593092e-58
     [188,] -7.965460e-59 -7.965460e-59 -7.965460e-59 -7.965460e-59 -7.965460e-59
     [189,] -3.982730e-59 -3.982730e-59 -3.982730e-59 -3.982730e-59 -3.982730e-59
     [190,] -1.991365e-59 -1.991365e-59 -1.991365e-59 -1.991365e-59 -1.991365e-59
     [191,] -9.956824e-60 -9.956824e-60 -9.956824e-60 -9.956824e-60 -9.956824e-60
     [192,] -4.978412e-60 -4.978412e-60 -4.978412e-60 -4.978412e-60 -4.978412e-60
     [193,] -2.489206e-60 -2.489206e-60 -2.489206e-60 -2.489206e-60 -2.489206e-60
     [194,] -1.244603e-60 -1.244603e-60 -1.244603e-60 -1.244603e-60 -1.244603e-60
     [195,] -6.223015e-61 -6.223015e-61 -6.223015e-61 -6.223015e-61 -6.223015e-61
     delta3
     [1,] -1.448317e-02
     [2,] -7.517260e-03
     [3,] -3.831203e-03
     [4,] -1.934208e-03
     [5,] -9.718137e-04
     [6,] -4.870916e-04
     [7,] -2.438429e-04
     [8,] -1.219958e-04
     [9,] -6.101653e-05
     [10,] -3.051292e-05
     [11,] -1.525762e-05
     [12,] -7.629104e-06
     [13,] -3.814625e-06
     [14,] -1.907330e-06
     [15,] -9.536698e-07
     [16,] -4.768360e-07
     [17,] -2.384183e-07
     [18,] -1.192092e-07
     [19,] -5.960463e-08
     [20,] -2.980232e-08
     [21,] -1.490116e-08
     [22,] -7.450580e-09
     [23,] -3.725290e-09
     [24,] -1.862645e-09
     [25,] -9.313226e-10
     [26,] -4.656613e-10
     [27,] -2.328306e-10
     [28,] -1.164153e-10
     [29,] -5.820766e-11
     [30,] -2.910383e-11
     [31,] -1.455192e-11
     [32,] -7.275958e-12
     [33,] -3.637979e-12
     [34,] -1.818989e-12
     [35,] -9.094947e-13
     [36,] -4.547474e-13
     [37,] -2.273737e-13
     [38,] -1.136868e-13
     [39,] -5.684342e-14
     [40,] -2.842171e-14
     [41,] -1.421085e-14
     [42,] -7.105427e-15
     [43,] -3.552714e-15
     [44,] -1.776357e-15
     [45,] -8.881784e-16
     [46,] -4.440892e-16
     [47,] -2.220446e-16
     [48,] -1.110223e-16
     [49,] -5.551115e-17
     [50,] -2.775558e-17
     [51,] -1.387779e-17
     [52,] -6.938894e-18
     [53,] -3.469447e-18
     [54,] -1.734723e-18
     [55,] -8.673617e-19
     [56,] -4.336809e-19
     [57,] -2.168404e-19
     [58,] -1.084202e-19
     [59,] -5.421011e-20
     [60,] -2.710505e-20
     [61,] -1.355253e-20
     [62,] -6.776264e-21
     [63,] -3.388132e-21
     [64,] -1.694066e-21
     [65,] -8.470329e-22
     [66,] -4.235165e-22
     [67,] -2.117582e-22
     [68,] -1.058791e-22
     [69,] -5.293956e-23
     [70,] -2.646978e-23
     [71,] -1.323489e-23
     [72,] -6.617445e-24
     [73,] -3.308722e-24
     [74,] -1.654361e-24
     [75,] -8.271806e-25
     [76,] -4.135903e-25
     [77,] -2.067952e-25
     [78,] -1.033976e-25
     [79,] -5.169879e-26
     [80,] -2.584939e-26
     [81,] -1.292470e-26
     [82,] -6.462349e-27
     [83,] -3.231174e-27
     [84,] -1.615587e-27
     [85,] -8.077936e-28
     [86,] -4.038968e-28
     [87,] -2.019484e-28
     [88,] -1.009742e-28
     [89,] -5.048710e-29
     [90,] -2.524355e-29
     [91,] -1.262177e-29
     [92,] -6.310887e-30
     [93,] -3.155444e-30
     [94,] -1.577722e-30
     [95,] -7.888609e-31
     [96,] -3.944305e-31
     [97,] -1.972152e-31
     [98,] -9.860761e-32
     [99,] -4.930381e-32
     [100,] -2.465190e-32
     [101,] -1.232595e-32
     [102,] -6.162976e-33
     [103,] -3.081488e-33
     [104,] -1.540744e-33
     [105,] -7.703720e-34
     [106,] -3.851860e-34
     [107,] -1.925930e-34
     [108,] -9.629650e-35
     [109,] -4.814825e-35
     [110,] -2.407412e-35
     [111,] -1.203706e-35
     [112,] -6.018531e-36
     [113,] -3.009266e-36
     [114,] -1.504633e-36
     [115,] -7.523164e-37
     [116,] -3.761582e-37
     [117,] -1.880791e-37
     [118,] -9.403955e-38
     [119,] -4.701977e-38
     [120,] -2.350989e-38
     [121,] -1.175494e-38
     [122,] -5.877472e-39
     [123,] -2.938736e-39
     [124,] -1.469368e-39
     [125,] -7.346840e-40
     [126,] -3.673420e-40
     [127,] -1.836710e-40
     [128,] -9.183550e-41
     [129,] -4.591775e-41
     [130,] -2.295887e-41
     [131,] -1.147944e-41
     [132,] -5.739719e-42
     [133,] -2.869859e-42
     [134,] -1.434930e-42
     [135,] -7.174648e-43
     [136,] -3.587324e-43
     [137,] -1.793662e-43
     [138,] -8.968310e-44
     [139,] -4.484155e-44
     [140,] -2.242078e-44
     [141,] -1.121039e-44
     [142,] -5.605194e-45
     [143,] -2.802597e-45
     [144,] -1.401298e-45
     [145,] -7.006492e-46
     [146,] -3.503246e-46
     [147,] -1.751623e-46
     [148,] -8.758115e-47
     [149,] -4.379058e-47
     [150,] -2.189529e-47
     [151,] -1.094764e-47
     [152,] -5.473822e-48
     [153,] -2.736911e-48
     [154,] -1.368456e-48
     [155,] -6.842278e-49
     [156,] -3.421139e-49
     [157,] -1.710569e-49
     [158,] -8.552847e-50
     [159,] -4.276424e-50
     [160,] -2.138212e-50
     [161,] -1.069106e-50
     [162,] -5.345529e-51
     [163,] -2.672765e-51
     [164,] -1.336382e-51
     [165,] -6.681912e-52
     [166,] -3.340956e-52
     [167,] -1.670478e-52
     [168,] -8.352390e-53
     [169,] -4.176195e-53
     [170,] -2.088097e-53
     [171,] -1.044049e-53
     [172,] -5.220244e-54
     [173,] -2.610122e-54
     [174,] -1.305061e-54
     [175,] -6.525304e-55
     [176,] -3.262652e-55
     [177,] -1.631326e-55
     [178,] -8.156631e-56
     [179,] -4.078315e-56
     [180,] -2.039158e-56
     [181,] -1.019579e-56
     [182,] -5.097894e-57
     [183,] -2.548947e-57
     [184,] -1.274474e-57
     [185,] -6.372368e-58
     [186,] -3.186184e-58
     [187,] -1.593092e-58
     [188,] -7.965460e-59
     [189,] -3.982730e-59
     [190,] -1.991365e-59
     [191,] -9.956824e-60
     [192,] -4.978412e-60
     [193,] -2.489206e-60
     [194,] -1.244603e-60
     [195,] -6.223015e-61
     >
     > ## testing cbind() i.e., the Rmpfr-method a bit here !
     > chkDn <- function(r) stopifnot(identical(
     + dimnames(r),
     + list(NULL, c("x", "log2x", eval(formals(eps_S53)$method)))))
     > chkDn2 <- function(r) stopifnot(identical(
     + colnames(r)[-(1:2)], eval(formals(eps_S53)$method)))
     > chkDn ( r1 <- cbind(x, log2x = log2(x), EA))
     > chkDn2( r2 <- selectMethod("cbind","mpfr")(x, log2x = log2(x), EA) )
     > ## are the same if names are *made* to match:
     > stopifnot(identical(r1, local({t <- r2; dimnames(t) <- dimnames(r1); t})))
     > chkDn2( r3 <- selectMethod("cbind","mpfr")(x, log2(x), epsAllS53(x)) )
     > chkDn2( r4 <- cbind(x, log2(x), epsAllS53(x)) )
     > stopifnot(identical(r1, local({t <- r4; dimnames(t) <- dimnames(r1); t})))
     >
     > U_S53mat <- function(x) cbind(x, log2x = log2(x), epsAllS53(x))
     > chkDn(Rx <- U_S53mat(x))
     > dimnames(Rx) # .. all empty, i.e., list(NULL, NULL) ... even though epsAll..() gives nice colnames
     [[1]]
     NULL
    
     [[2]]
     [1] "x" "log2x" "direct" "log1pExct" "deltaFull" "delta1"
     [7] "delta2" "delta3"
    
     > ##
     > if(FALSE) ## no longer needed !
     + ## A version with correct colnames :
     + U_S53mat <- function(x) {
     + eA <- epsAllS53(x)
     + r <- cbind(x, log2(x), eA)
     + colnames(r) <- c("x", "log2.x", colnames(eA))
     + r
     + }
     >
     > stopifnot(Rx[,-(1:2)] < 0)
     >
     > op <- options(width = 100)
     > print(asNumeric(Rx), digits=5) ; options(op)
     x log2x direct log1pExct deltaFull delta1 delta2 delta3
     [1,] 8.0000e+00 3.0 -1.5051e-02 -1.4585e-02 -1.4479e-02 -1.5385e-02 -1.4423e-02 -1.4483e-02
     [2,] 1.1314e+01 3.5 -7.6647e-03 -7.5454e-03 -7.5170e-03 -7.7519e-03 -7.5097e-03 -7.5173e-03
     [3,] 1.6000e+01 4.0 -3.8687e-03 -3.8385e-03 -3.8312e-03 -3.8911e-03 -3.8303e-03 -3.8312e-03
     [4,] 2.2627e+01 4.5 -1.9437e-03 -1.9361e-03 -1.9342e-03 -1.9493e-03 -1.9341e-03 -1.9342e-03
     [5,] 3.2000e+01 5.0 -9.7419e-04 -9.7229e-04 -9.7181e-04 -9.7561e-04 -9.7180e-04 -9.7181e-04
     [6,] 4.5255e+01 5.5 -4.8769e-04 -4.8721e-04 -4.8709e-04 -4.8804e-04 -4.8709e-04 -4.8709e-04
     [7,] 6.4000e+01 6.0 -2.4399e-04 -2.4387e-04 -2.4384e-04 -2.4408e-04 -2.4384e-04 -2.4384e-04
     [8,] 9.0510e+01 6.5 -1.2203e-04 -1.2200e-04 -1.2200e-04 -1.2206e-04 -1.2200e-04 -1.2200e-04
     [9,] 1.2800e+02 7.0 -6.1026e-05 -6.1018e-05 -6.1017e-05 -6.1031e-05 -6.1017e-05 -6.1017e-05
     [10,] 1.8102e+02 7.5 -3.0515e-05 -3.0513e-05 -3.0513e-05 -3.0517e-05 -3.0513e-05 -3.0513e-05
     [11,] 2.5600e+02 8.0 -1.5258e-05 -1.5258e-05 -1.5258e-05 -1.5259e-05 -1.5258e-05 -1.5258e-05
     [12,] 3.6204e+02 8.5 -7.6292e-06 -7.6291e-06 -7.6291e-06 -7.6293e-06 -7.6291e-06 -7.6291e-06
     [13,] 5.1200e+02 9.0 -3.8147e-06 -3.8146e-06 -3.8146e-06 -3.8147e-06 -3.8146e-06 -3.8146e-06
     [14,] 7.2408e+02 9.5 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06
     [15,] 1.0240e+03 10.0 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07
     [16,] 1.4482e+03 10.5 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07
     [17,] 2.0480e+03 11.0 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07
     [18,] 2.8963e+03 11.5 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07
     [19,] 4.0960e+03 12.0 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08
     [20,] 5.7926e+03 12.5 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08
     [21,] 8.1920e+03 13.0 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08
     [22,] 1.1585e+04 13.5 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09
     [23,] 1.6384e+04 14.0 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09
     [24,] 2.3170e+04 14.5 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09
     [25,] 3.2768e+04 15.0 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10
     [26,] 4.6341e+04 15.5 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10
     [27,] 6.5536e+04 16.0 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10
     [28,] 9.2682e+04 16.5 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10
     [29,] 1.3107e+05 17.0 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11
     [30,] 1.8536e+05 17.5 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11
     [31,] 2.6214e+05 18.0 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11
     [32,] 3.7073e+05 18.5 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12
     [33,] 5.2429e+05 19.0 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12
     [34,] 7.4146e+05 19.5 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12
     [35,] 1.0486e+06 20.0 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13
     [36,] 1.4829e+06 20.5 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13
     [37,] 2.0972e+06 21.0 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13
     [38,] 2.9658e+06 21.5 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13
     [39,] 4.1943e+06 22.0 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14
     [40,] 5.9316e+06 22.5 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14
     [41,] 8.3886e+06 23.0 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14
     [42,] 1.1863e+07 23.5 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15
     [43,] 1.6777e+07 24.0 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15
     [44,] 2.3727e+07 24.5 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15
     [45,] 3.3554e+07 25.0 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16
     [46,] 4.7453e+07 25.5 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16
     [47,] 6.7109e+07 26.0 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16
     [48,] 9.4906e+07 26.5 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16
     [49,] 1.3422e+08 27.0 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17
     [50,] 1.8981e+08 27.5 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17
     [51,] 2.6844e+08 28.0 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17
     [52,] 3.7963e+08 28.5 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18
     [53,] 5.3687e+08 29.0 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18
     [54,] 7.5925e+08 29.5 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18
     [55,] 1.0737e+09 30.0 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19
     [56,] 1.5185e+09 30.5 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19
     [57,] 2.1475e+09 31.0 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19
     [58,] 3.0370e+09 31.5 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19
     [59,] 4.2950e+09 32.0 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20
     [60,] 6.0740e+09 32.5 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20
     [61,] 8.5899e+09 33.0 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20
     [62,] 1.2148e+10 33.5 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21
     [63,] 1.7180e+10 34.0 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21
     [64,] 2.4296e+10 34.5 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21
     [65,] 3.4360e+10 35.0 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22
     [66,] 4.8592e+10 35.5 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22
     [67,] 6.8719e+10 36.0 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22
     [68,] 9.7184e+10 36.5 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22
     [69,] 1.3744e+11 37.0 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23
     [70,] 1.9437e+11 37.5 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23
     [71,] 2.7488e+11 38.0 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23
     [72,] 3.8874e+11 38.5 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24
     [73,] 5.4976e+11 39.0 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24
     [74,] 7.7747e+11 39.5 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24
     [75,] 1.0995e+12 40.0 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25
     [76,] 1.5549e+12 40.5 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25
     [77,] 2.1990e+12 41.0 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25
     [78,] 3.1099e+12 41.5 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25
     [79,] 4.3980e+12 42.0 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26
     [80,] 6.2198e+12 42.5 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26
     [81,] 8.7961e+12 43.0 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26
     [82,] 1.2440e+13 43.5 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27
     [83,] 1.7592e+13 44.0 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27
     [84,] 2.4879e+13 44.5 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27
     [85,] 3.5184e+13 45.0 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28
     [86,] 4.9758e+13 45.5 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28
     [87,] 7.0369e+13 46.0 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28
     [88,] 9.9516e+13 46.5 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28
     [89,] 1.4074e+14 47.0 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29
     [90,] 1.9903e+14 47.5 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29
     [91,] 2.8147e+14 48.0 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29
     [92,] 3.9807e+14 48.5 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30
     [93,] 5.6295e+14 49.0 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30
     [94,] 7.9613e+14 49.5 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30
     [95,] 1.1259e+15 50.0 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31
     [96,] 1.5923e+15 50.5 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31
     [97,] 2.2518e+15 51.0 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31
     [98,] 3.1845e+15 51.5 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32
     [99,] 4.5036e+15 52.0 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32
     [100,] 6.3691e+15 52.5 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32
     [101,] 9.0072e+15 53.0 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32
     [102,] 1.2738e+16 53.5 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33
     [103,] 1.8014e+16 54.0 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33
     [104,] 2.5476e+16 54.5 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33
     [105,] 3.6029e+16 55.0 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34
     [106,] 5.0952e+16 55.5 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34
     [107,] 7.2058e+16 56.0 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34
     [108,] 1.0190e+17 56.5 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35
     [109,] 1.4412e+17 57.0 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35
     [110,] 2.0381e+17 57.5 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35
     [111,] 2.8823e+17 58.0 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35
     [112,] 4.0762e+17 58.5 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36
     [113,] 5.7646e+17 59.0 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36
     [114,] 8.1524e+17 59.5 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36
     [115,] 1.1529e+18 60.0 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37
     [116,] 1.6305e+18 60.5 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37
     [117,] 2.3058e+18 61.0 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37
     [118,] 3.2610e+18 61.5 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38
     [119,] 4.6117e+18 62.0 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38
     [120,] 6.5219e+18 62.5 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38
     [121,] 9.2234e+18 63.0 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38
     [122,] 1.3044e+19 63.5 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39
     [123,] 1.8447e+19 64.0 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39
     [124,] 2.6088e+19 64.5 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39
     [125,] 3.6893e+19 65.0 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40
     [126,] 5.2175e+19 65.5 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40
     [127,] 7.3787e+19 66.0 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40
     [128,] 1.0435e+20 66.5 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41
     [129,] 1.4757e+20 67.0 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41
     [130,] 2.0870e+20 67.5 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41
     [131,] 2.9515e+20 68.0 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41
     [132,] 4.1740e+20 68.5 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42
     [133,] 5.9030e+20 69.0 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42
     [134,] 8.3480e+20 69.5 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42
     [135,] 1.1806e+21 70.0 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43
     [136,] 1.6696e+21 70.5 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43
     [137,] 2.3612e+21 71.0 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43
     [138,] 3.3392e+21 71.5 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44
     [139,] 4.7224e+21 72.0 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44
     [140,] 6.6784e+21 72.5 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44
     [141,] 9.4447e+21 73.0 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44
     [142,] 1.3357e+22 73.5 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45
     [143,] 1.8889e+22 74.0 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45
     [144,] 2.6714e+22 74.5 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45
     [145,] 3.7779e+22 75.0 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46
     [146,] 5.3427e+22 75.5 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46
     [147,] 7.5558e+22 76.0 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46
     [148,] 1.0685e+23 76.5 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47
     [149,] 1.5112e+23 77.0 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47
     [150,] 2.1371e+23 77.5 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47
     [151,] 3.0223e+23 78.0 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47
     [152,] 4.2742e+23 78.5 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48
     [153,] 6.0446e+23 79.0 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48
     [154,] 8.5484e+23 79.5 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48
     [155,] 1.2089e+24 80.0 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49
     [156,] 1.7097e+24 80.5 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49
     [157,] 2.4179e+24 81.0 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49
     [158,] 3.4194e+24 81.5 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50
     [159,] 4.8357e+24 82.0 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50
     [160,] 6.8387e+24 82.5 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50
     [161,] 9.6714e+24 83.0 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50
     [162,] 1.3677e+25 83.5 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51
     [163,] 1.9343e+25 84.0 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51
     [164,] 2.7355e+25 84.5 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51
     [165,] 3.8686e+25 85.0 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52
     [166,] 5.4710e+25 85.5 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52
     [167,] 7.7371e+25 86.0 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52
     [168,] 1.0942e+26 86.5 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53
     [169,] 1.5474e+26 87.0 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53
     [170,] 2.1884e+26 87.5 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53
     [171,] 3.0949e+26 88.0 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53
     [172,] 4.3768e+26 88.5 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54
     [173,] 6.1897e+26 89.0 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54
     [174,] 8.7536e+26 89.5 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54
     [175,] 1.2379e+27 90.0 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55
     [176,] 1.7507e+27 90.5 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55
     [177,] 2.4759e+27 91.0 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55
     [178,] 3.5014e+27 91.5 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56
     [179,] 4.9518e+27 92.0 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56
     [180,] 7.0028e+27 92.5 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56
     [181,] 9.9035e+27 93.0 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56
     [182,] 1.4006e+28 93.5 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57
     [183,] 1.9807e+28 94.0 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57
     [184,] 2.8011e+28 94.5 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57
     [185,] 3.9614e+28 95.0 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58
     [186,] 5.6023e+28 95.5 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58
     [187,] 7.9228e+28 96.0 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58
     [188,] 1.1205e+29 96.5 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59
     [189,] 1.5846e+29 97.0 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59
     [190,] 2.2409e+29 97.5 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59
     [191,] 3.1691e+29 98.0 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60
     [192,] 4.4818e+29 98.5 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60
     [193,] 6.3383e+29 99.0 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60
     [194,] 8.9636e+29 99.5 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60
     [195,] 1.2677e+30 100.0 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61
     >
     > matplot(Rx[,"x"], - Rx[,-(1:2)], type = "b", log="xy")
     Error in xy.coords(x, y, xlabel, ylabel, log = log) :
     'list' object cannot be coerced to type 'double'
     Calls: matplot -> xy.coords
     Execution halted
Flavor: r-oldrel-windows-ix86+x86_64