# NOT RUN {
x <- (0:100)/8 # (have exact binary representation)
stopifnot( all.equal(besselY(x, 0), bY0 <- y0(x))
, all.equal(besselJ(x, 1), bJ1 <- j1(x))
, all.equal(yn(0,x), bY0)
, all.equal(jn(1,x), bJ1)
)
# }
Run the code above in your browser using DataLab