# NOT RUN {
x <- c(1:10, 20, 50, 100, 100000)
nu <- c(1, 10, 20, 50, 10^(2:10))
# }
# NOT RUN {
<!-- %% FIXME: do something more sensible -- compare with besselI() and *K() -->
# }
# NOT RUN {
sapply(0:4, function(k.)
sapply(nu, function(n.)
besselI.nuAsym(x, nu=n., k.max = k., log = TRUE)))
sapply(0:4, function(k.)
sapply(nu, function(n.)
besselK.nuAsym(x, nu=n., k.max = k., log = TRUE)))
# }
Run the code above in your browser using DataLab