Ejemplo n.º 1
0
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) );
}
Ejemplo n.º 2
0
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) );

}