similar to: besselK

Assistance, Im a student intrested in using R in my learning and research work in option pricing however i have a problem with besselK function In R. Would you assit me in computing the besselK of third kind of a complex number in R. Any code or suggestion will be highly appriceiated eg besselK(2,10) works well.. but besselK(2,10i) doesnt work !! im supprised it works in MATLAB but NOT

The besselK() function knows about overflows/underflows internally; there is a constant xmax_BESS_K in src/nmath/bessel.h (and referred to only in bessel_k.c), equal to 705.342, which is checked if expon.scaled is FALSE. (The equivalent number for bessel_i.c is 709, defined as exparg_BESS in bessel.h.) However, besselK(x) silently returns +Inf if x>705.342. This behavior is reasonable for

hello all i searched the archives and couldn't get a solution to the following question. i have the following function: F=function(z,v) { if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1)) } and while testing whether it works (ie F(z,v) is approx = besselK(z,1+v)/besselK(z,V). the recursion formula allows one to calculate besselK(z,1+v)/besselK(z,V) for large values of z )

[cc'ing to R-devel] It's moderately obscure to me too, but ... There are two sets of bessel code in the package, one in C++ from Chris Bond and the other in FORTRAN from Netlib (I think). The FORTRAN code is the one that's giving trouble, you might try just removing the FORTRAN code from the src directory and trying again.

I am looking for a method of dealing with the modified Bessel function K_\nu(x) for large \nu. The besselK function implementation of this allows for dealing with large values of x by allowing for exponential scaling, but there is no facility for dealing with large \nu. What would work for me would be an lbesselK function in the manner of lgamma which returned the log of K_\nu(x) for large

Hi: There is a Modified Bessel function (third kind, real order) besselK in R. Is there a C version of this function? I can only find a Fortran function rkbesel on http://www.netlib.org/specfun/. However, this rkbesel fortran function is the second kind. Thank you! liu _______________________________ Declare Yourself - Register online to vote today!

hello all i searched the archives and couldn't get a solution to the following question. i have the following function: F=function(z,v) { if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1)) } and while testing whether it works (ie F(z,v) is approx = besselK(z,1+v)/besselK(z,V). the recursion formula allows one to calculate besselK(z,1+v)/besselK(z,V) for large values of z )

Hello, I am searching for code in C++ or fortran for Modified Bessel function of third kind (fractional or real order). Can someone help me? Thank you --------------------------------- Découvrez un nouveau moyen de poser toutes vos questions quelque soit le sujet ! Yahoo! Questions/Réponses pour partager vos connaissances, vos opinions et vos expériences. Cliquez ici.

Hi all, Running R 0.90.1 on a RH 6.1 system. Installation of the integrate_2.1-2 package went smoothly. My code contains a loop in which integrate() is called several times in each pass. I get a segmentation fault after what seems to be a random number of calls to integrate(). Debug output shows: Program received signal SIGSEGV, Segmentation fault. promiseArgs (el=0x40276414,

In order to fit a probability distribution proposed by Sichel [Journal of the Royal Statistical Society. Series A (General), Vol. 137, No. 1. (1974), pp. 25-34], I need a modified Bessel function of the 2nd kind. I notice that the base package of "R" only has modified Bessel functions of the 1st and 3rd kind. Does a modified Bessel function of the 2nd kind exist anywhere? Many

I've rolled up R-1.6.1.tgz a short while ago. This is a minor upgrade, fixing an assortment of minor bugs, and one major one: 1.6.0 introduced a memory leak in deparsing, causing much trouble for people running simulations. You can get it from the developer site at http://cvs.r-project.org/pub/CRAN/src/base/R-1.6.1.tgz or wait for it to be mirrored at a CRAN site near you. Binaries for

I've rolled up R-1.6.1.tgz a short while ago. This is a minor upgrade, fixing an assortment of minor bugs, and one major one: 1.6.0 introduced a memory leak in deparsing, causing much trouble for people running simulations. You can get it from the developer site at http://cvs.r-project.org/pub/CRAN/src/base/R-1.6.1.tgz or wait for it to be mirrored at a CRAN site near you. Binaries for

Dear R-Group: Recently, I ran into a problem. I was using a function called "I.1", which evaluates the first-order modified Bessel function of the first kind, in the package "CircStats". This function is not vectorized, since it uses a couple of "if" conditions. However, when I called this function with a vector argument, I got no error/warning messages in

I've put up R-0.63.0.tgz up for FTP from Auckland some minutes ago. As usual, don't get it from there unless you are desperate, but wait for it to be mirrored at a CRAN site near you within a day or two. For those who *are* desperate, I've left a copy in ftp://blueberry.kubism.ku.dk/pub/R-devel/R-0.63.tgz (Be gentle, that's my desktop PC!) Here is the relevant part of the NEWS

I've put up R-0.63.0.tgz up for FTP from Auckland some minutes ago. As usual, don't get it from there unless you are desperate, but wait for it to be mirrored at a CRAN site near you within a day or two. For those who *are* desperate, I've left a copy in ftp://blueberry.kubism.ku.dk/pub/R-devel/R-0.63.tgz (Be gentle, that's my desktop PC!) Here is the relevant part of the NEWS

Example run and stack trace: wazor /s/src/stat/R-1.8.0/tests/Examples $ ../../bin/R --no-save < base-Ex.R R : Copyright 2003, The R Development Core Team Version 1.8.0 (2003-10-08) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a

Hi all I seem to have this problem whereby if I delete an email, then the deleted email does not appear in the Trash folder. I tried setting dovecot mail_debug = yes, but I find in syslog, mail.log etc that there is very limited information. If anyone could assit, I would be most greatful. Kind Regards Brent Clark P.s. I use Thunderbird as the mail client

R-1.7.0 built on NetBSD 1.6, but the validation test suite failed: Machinetype: Intel Pentium III (600 MHz); NetBSD 1.6 (GENERIC) Remote gcc version: gcc (GCC) 3.2.2 Remote g++ version: g++ (GCC) 3.2.2 Configure environment: CC=gcc CXX=g++ LDFLAGS=-Wl,-rpath,/usr/local/lib make[5]: Entering directory `/local/build/R-1.7.0/src/library' >>> Building/Updating