19 lines
400 B
C
19 lines
400 B
C
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) );
|
|
}
|