diff --git a/README b/README
index f121543fd4bb550dea713867f895c684c53acb5d..bb08f37598870b341df5c5b60fbef8b80f8b2b6e 100644
--- a/README
+++ b/README
@@ -17,3 +17,7 @@ is that the GSL library must be loaded and accessible.
 access to the functions.  Defining the functions you need yourself is
 very easy.
 
+======================
+Needs work
+bessel-sequence-Jnu in bessel.lisp.
+Why do sequence-returns crash on declination?
diff --git a/general/mathematical.lisp b/general/mathematical.lisp
index 46d7d19b972569748c938ab40527d8d6a7346708..161bc98d3290b46af4bdec708c5577969dc5ccde 100644
--- a/general/mathematical.lisp
+++ b/general/mathematical.lisp
@@ -3,13 +3,12 @@
 ; description: Mathematical functions                    
 ; date:        Wed Mar  8 2006 - 22:09                   
 ; author:      Liam M. Healy
-; modified:    Thu Mar  9 2006 - 23:10
+; modified:    Sat Mar 18 2006 - 15:13
 ;********************************************************
 
 (in-package :gsl)
 
-(eval-when (:compile-toplevel :load-toplevel :execute)
-  (export '(nanp infinityp finitep log+1 exp-1 hypotenuse approximately=)))
+(export '(nanp infinityp finitep log+1 exp-1 hypotenuse approximately=))
 
 ;;;;****************************************************************************
 ;;; Infinities and Not-a-number
diff --git a/special-functions/airy.lisp b/special-functions/airy.lisp
index 82fb920162b143ae6de23b3647a2f77a66ec5352..94fb4fcb1729af6f66bc56a5e6797fde5dfb573e 100644
--- a/special-functions/airy.lisp
+++ b/special-functions/airy.lisp
@@ -3,20 +3,16 @@
 ; description: Airy functions                            
 ; date:        Fri Mar 17 2006 - 18:41                   
 ; author:      Liam M. Healy
-; modified:    Sat Mar 18 2006 - 00:21
+; modified:    Sat Mar 18 2006 - 15:21
 ;********************************************************
 
 (in-package :gsl)
 
-(eval-when (:compile-toplevel :load-toplevel :execute)
-  (export
-   '()))
-
 ;;;;****************************************************************************
 ;;;; Airy functions
 ;;;;****************************************************************************
 
-(defun-sf airy-ai ((x :double))
+(defun-sf airy-Ai ((x :double))
   "gsl_sf_airy_Ai_e"
   :return (sf-result)
   :documentation "The Airy function Ai(x)."
@@ -26,70 +22,69 @@
 ;;; 0.24308135437540998d0
 ;;; 6.288670879282072d-17
 
-(defun-sf airy-bi ((x :double))
+(defun-sf airy-Bi ((x :double))
   "gsl_sf_airy_Bi_e"
   :return (sf-result)
   :documentation "The Airy function Bi(x)."
   :mode t)
 
-(defun-sf airy-ai-scaled ((x :double))
+(defun-sf airy-Ai-scaled ((x :double))
   "gsl_sf_airy_Ai_scaled_e"
   :return (sf-result)
   :documentation "The scaled Airy function @math{S_A(x) Ai(x)}.  For @math{x>0} the scaling factor @math{S_A(x)} is @c{$\exp(+(2/3) x^{3/2})$} @math{\exp(+(2/3) x^(3/2))}, and is 1 for @math{x<0}."
   :mode t)
 
-(defun-sf airy-bi-scaled ((x :double))
+(defun-sf airy-Bi-scaled ((x :double))
   "gsl_sf_airy_Bi_scaled_e"
   :return (sf-result)
   :documentation "The scaled Airy function @math{S_B(x) Bi(x)}.  For @math{x>0} the scaling factor @math{S_B(x)} is @c{$\exp(-(2/3) x^{3/2})$} @math{exp(-(2/3) x^(3/2))}, and is 1 for @math{x<0}."
   :mode t)
 
-(defun-sf airy-ai-deriv ((x :double))
+(defun-sf airy-Ai-deriv ((x :double))
   "gsl_sf_airy_Ai_deriv_e"
   :return (sf-result)
   :documentation "The Airy function derivative Ai'(x)."
   :mode t)
 
-(defun-sf airy-bi-deriv ((x :double))
+(defun-sf airy-Bi-deriv ((x :double))
   "gsl_sf_airy_Bi_deriv_e"
   :return (sf-result)
   :documentation "The Airy function derivative Bi'(x)."
   :mode t)
 
-(defun-sf airy-ai-deriv-scaled ((x :double))
+(defun-sf airy-Ai-deriv-scaled ((x :double))
   "gsl_sf_airy_Ai_deriv_scaled_e"
   :return (sf-result)
   :documentation "The scaled Airy function derivative S_A(x) Ai'(x).  For @math{x>0} the scaling factor @math{S_A(x)} is @c{$\exp(+(2/3) x^{3/2})$} @math{\exp(+(2/3) x^(3/2))}, and is 1 for @math{x<0}."
   :mode t)
 
-(defun-sf airy-bi-deriv-scaled ((x :double))
+(defun-sf airy-Bi-deriv-scaled ((x :double))
   "gsl_sf_airy_Bi_deriv_scaled_e"
   :return (sf-result)
   :documentation "The scaled Airy function derivative S_B(x) Bi'(x).  For @math{x>0} the scaling factor @math{S_B(x)} is @c{$\exp(-(2/3) x^{3/2})$} @math{exp(-(2/3) x^(3/2))}, and is 1 for @math{x<0}."
   :mode t)
 
-(defun-sf airy-zero-ai ((s :unsigned-int))
+(defun-sf airy-zero-Ai ((s :unsigned-int))
   "gsl_sf_airy_zero_Ai_e"
   :return (sf-result)
   :documentation "The location of the @var{s}-th zero of the Airy function @math{Ai(x)}."
   :mode t)
 
-(defun-sf airy-zero-bi ((s :unsigned-int))
+(defun-sf airy-zero-Bi ((s :unsigned-int))
   "gsl_sf_airy_zero_Bi_e"
   :return (sf-result)
   :documentation "The location of the @var{s}-th zero of the Airy function @math{Bi(x)}."
   :mode t)
 
-(defun-sf airy-zero-ai-deriv ((s :unsigned-int))
+(defun-sf airy-zero-Ai-deriv ((s :unsigned-int))
   "gsl_sf_airy_zero_Ai_deriv_e"
   :return (sf-result)
   :documentation "The location of the @var{s}-th zero of the Airy
 function derivative @math{Ai'(x)}."
   :mode t)
 
-(defun-sf airy-zero-bi-deriv ((s :unsigned-int))
+(defun-sf airy-zero-Bi-deriv ((s :unsigned-int))
   "gsl_sf_airy_zero_Bi_deriv_e"
   :return (sf-result)
   :documentation "The location of the @var{s}-th zero of the Airy function derivative @math{Bi'(x)}."
   :mode t)
-