Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: RE: RE: Bessel


From   "Nick Cox" <n.j.cox@durham.ac.uk>
To   <statalist@hsphsun2.harvard.edu>
Subject   st: RE: RE: Bessel
Date   Sat, 1 Oct 2005 17:43:21 +0100

Bobby Gutierrez needed some Bessel function 
for some purpose and wrote an undocumented helper
program. I can't remember which and can't find 
the code. 

I once wrote a rather specific Stata program for I_0, 
not what you specify here. 

1.5708 means _pi/2 here. 

Odds are you need to transliterate formulae
from Abramowitz and Stegun's handbook to Stata or Mata. 

Nick 
n.j.cox@durham.ac.uk 

Jhilbe@aol.com
 
> To be more specific, has anyone  made a Bessel of the 2nd 
> kind function in 
> Stata? An approximation is  sqrt(1.5708/x) * 
> sin(x-0.7854(2y+1)) for values 
> closer to the max values  and
> -[ (2^n)(n-1)! /pi ] * x^(-n) for values closer to 0. If 
> someone has a  
> better approximation, or the "exact" function in Stata, I'd 
> love to see  it.  

*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index