search for: besselk

Assistance, besselK- complex number problem 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 appriceiate...

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.) Howeve...

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 in R. rgds ivivi mwaniki student...

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 ) i got the following error : "evaluation nested too deeply: infinite recursion / options(expressions=)?" i.e. zz=1 nn=1000 R=array(0,dim=nn) for (i in 1:nn) { print(...

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 \nu. Does...

...essel' >> ** Removing >> '/home/baptiste/R/i486-pc-linux-gnu-library/2.7/complexbessel' > > > > ________________________________________ > From: Ben Bolker [bolker at ufl.edu] > Sent: 24 November 2008 15:14 > To: Auguie, Baptiste > Subject: Re: [R] besselK > > I took a little time to build a "complexbessel" package. > It's barely documented ... but see if this helps. > It would be nice to test this for speed and accuracy, > document it further, etc ... > > Ben Bolker > > > > baptiste auguie wrote...

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 ) i got the following error : "evaluation nested too deeply: infinite recursion / options(expressions=)?" i.e. zz=1 nn=1000 R=array(0,dim=nn) for (i in 1:nn) { print(...

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, 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

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

...#39; (smoothing and local MODern REGression methods) contributed by B. D. Ripley. o par() has a `no.readonly = FALSE' argument which allows more sensible op <- par(no.readonly = TRUE); on.exit(par(op)). o Real Bessel functions of 1st to 3rd kind, of arbitrary order: besselI(), besselK(), besselJ(), besselY() are the I(), K(), J(), and Y() Bessel functions. o New conflicts() function from B.D.Ripley o uniroot has a new `maxiter' argument and returns #{iter} and precision. o New option `show.coef.Pvalue' (default: TRUE). If FALSE, print.summary.[g]lm does...

...#39; (smoothing and local MODern REGression methods) contributed by B. D. Ripley. o par() has a `no.readonly = FALSE' argument which allows more sensible op <- par(no.readonly = TRUE); on.exit(par(op)). o Real Bessel functions of 1st to 3rd kind, of arbitrary order: besselI(), besselK(), besselJ(), besselY() are the I(), K(), J(), and Y() Bessel functions. o New conflicts() function from B.D.Ripley o uniroot has a new `maxiter' argument and returns #{iter} and precision. o New option `show.coef.Pvalue' (default: TRUE). If FALSE, print.summary.[g]lm does...

...t, as documented. Previously it was read unless R_ENVIRON pointed to an actual file. o Startup.Rd described the proceesing under Unix-alikes but incorrectly implied it happened that way on the Windows and MacOS ports. Neither use Renviron.site, for example. o besselK(x,*) now returns 0 instead of Inf for large x. (PR#2179) o The Tcl console code didn't work with Tcl/TK 8.0, and has been #ifdef'd out. (PR#2090) o format.AsIs() was not handling matrices. o sd() was not passing na.rm to var() for matrices and data frames....

...t, as documented. Previously it was read unless R_ENVIRON pointed to an actual file. o Startup.Rd described the proceesing under Unix-alikes but incorrectly implied it happened that way on the Windows and MacOS ports. Neither use Renviron.site, for example. o besselK(x,*) now returns 0 instead of Inf for large x. (PR#2179) o The Tcl console code didn't work with Tcl/TK 8.0, and has been #ifdef'd out. (PR#2090) o format.AsIs() was not handling matrices. o sd() was not passing na.rm to var() for matrices and data frames....

...iodic [1] 1 0 1 0 1 0 1 0 1 0 1 0 1 0 > x %/% 5 [1] -1 0 0 0 0 0 1 1 1 1 1 2 2 2 > > > > cleanEx(); ..nameEx <- "Bessel" > > ### * Bessel > > ### Name: Bessel > ### Title: Bessel Functions > ### Aliases: bessel Bessel besselI besselJ besselK besselY gammaCody > ### Keywords: math > > ### ** Examples > > nus <- c(0:5,10,20) > > x <- seq(0,4, len= 501) > plot(x,x, ylim = c(0,6), ylab="",type='n', main = "Bessel Functions I_nu(x)") Bus Error(coredump) wazor /s/src/stat/R-1.8....

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