# NOT RUN {
x <- c(1:10, 20, 50, 100^(2:10))
nu <- c(1, 10, 20, 50, 100)
# }
# NOT RUN {
<!-- %% FIXME: do something more sensible -- compare with besselI() and *K() -->
# }
# NOT RUN {
r <- lapply(c(0:4,10,20), function(k.)
sapply(nu, function(n.)
besselIasym(x, nu=n., k.max = k., log = TRUE)))
warnings()
# }
# NOT RUN {
<!-- %% NaN's ... seem a bit scattered. FIXME ?? -->
# }
# NOT RUN {
try( # needs improvement in R [or a local workaround]
besselIasym(10000*(1+1i), nu=200, k.max=20, log=TRUE)
) # Error in log1p(-d) : unimplemented complex function
# }
Run the code above in your browser using DataLab