ru.algorithms

 
 - RU.ALGORITHMS ----------------------------------------------------------------
 From : Vladislav Terehov                    2:5030/1307.25 07 Apr 2003  23:50:36
 To : All
 Subject : Интегральный синус/косинус
 -------------------------------------------------------------------------------- 
 

 Дайте алгоpитм / исходняк вычисления сабжа с максимально возможной для double
 (c++) точности.
 
 Вот чего я надумал:
 
 double si(const double &x) {
 
    const double eps = 0.000000000000001;
 
    double res1=0.0;
    double cur = 0;
 
    ulong n = 0; ulong k=0;
 
    double one = 1.0;
 
    do {
 
       k = 2*n+1;
 
       cur = pow(x, k) / ( k * fac(k));
       res1 += one * cur;
 
       ++n;
       one = -one;
    } while (cur >= eps);
 
    return res1;
 }
 
 Тока не pаботает оно как надо, если x >= 0.5 =( cur = 1.#INF
 wbr, Vlad.
 --- GoldEd+/W32 1.1.5-30228
  * Origin: ...trying 2 fun (2:5030/1307.25)
 
 

Вернуться к списку тем, сортированных по: возрастание даты  уменьшение даты  тема  автор 

Тема:    Автор:    Дата:    Интегральный синус/косинус   Vladislav Terehov   07 Apr 2003 23:50:36   Интегральный синус/косинус   Evgenij Masherov   08 Apr 2003 15:43:31