예제 #1
0
/// Sine integral.
double
Si(double x)
{
  gsl_sf_result result;
  int stat = gsl_sf_Si_e(x, &result);
  if (stat != GSL_SUCCESS)
    {
      std::ostringstream msg("Error in Si:");
      msg << " x=" << x;
      throw std::runtime_error(msg.str());
    }
  else
    return result.val;
}
예제 #2
0
파일: sinint.c 프로젝트: tommyliu/visionPJ1
double gsl_sf_Si(const double x)
{
  EVAL_RESULT(gsl_sf_Si_e(x, &result));
}