+++ /dev/null
-extern float A[], B[];
-extern float MAXNUMF;
-float chbevlf(float, float *, int);
-float expf(float), i1f(float), logf(float), sqrtf(float);
-
-float k1f(float xx)
-{
- float x, y;
-
- x = xx;
- if( x <= 2.0 )
- {
- y = x * x - 2.0;
- y = logf( 0.5f * x ) * i1f(x) + chbevlf( y, A, 7 ) / x;
- return( y );
- }
- return( expf(-x) * chbevlf( (float)(8.0/x - 2.0), B, 10 ) / sqrtf(x) );
-}