#include <iostream>
#include "bessely.h"
extern "C" {
#include "nr.h"
}
/* 
% 
% Written by:
% -- 
% John L. Weatherwax                2005-07-30
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----
*/

double bessely( double nu, double z ){
  float rj,ry,rjp,ryp;

  bessjy( (float) z, (float) nu, &rj, &ry, &rjp, &ryp ); 

  return ((double) ry);
}

#ifdef _BUILD_TESTS_
int main(){
  double nu, z;

  nu = 1; z = 1; 
  std::cout << "bessely( " << nu << ", " << z << ") = " << bessely( nu, z ) << std::endl; 

  nu = 1.5; z = 2; 
  std::cout << "bessely( " << nu << ", " << z << ") = " << bessely( nu, z ) << std::endl; 
}
#endif