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) ); }
float k1f(float xx) { float x, y; x = xx; if( x <= MINNUMF ) /* { mtherr( "k1f", DOMAIN ); return( MAXNUMF ); } */ assert(0); if( x <= 2.0f ) { y = x * x - 2.0f; y = logf( 0.5f * x ) * i1f(x) + chbevlf( y, A_k1f, 7 ) / x; return( y ); } return( expf(-x) * chbevlf( 8.0f/x - 2.0f, B_k1f, 10 ) / sqrtf(x) ); }