Проблемы с соиЬ (C++)

Question

Я имею трудное время выяснить, что здесь не так:

#include <iostream> 
#include <cmath> 
#include <iomanip> 

using namespace std; 

double fact(double); 
double sinTaylor(double); 
double cosTaylor(double); 

int main() 
{ 
    double number, sineOfnumber, cosineOfnumber; 

    cout << "Enter a number, then I will calculate the sine and cosine of this number" << endl; 

    cin >> number; 

    sineOfnumber = sinTaylor(number); 
    cosineOfnumber = cosTaylor(number); 

    cout << fixed << endl; 
    cout << cosineOfnumber << endl; 
    cout << sineOfnumber << endl; 

    return 0; 
} 

double fact(double n) 
{ 
    double product = 1; 
    while(n > 1) 
    product *= n--; 
    return product; 
} 

double sinTaylor(double x) 
{ 
    double currentIteration, sumSine; 

    for(double n = 0; n < 5; n++) 
    { 
     currentIteration = pow(-1, n)*pow(x, 2*n+1)/fact(2*n+1); 
     sumSine += currentIteration; 
    } 
    return sumSine; 
} 

double cosTaylor(double y) 
{ 
    double currentIteration, sumCosine; 

    for(double n = 0; n < 5; n++) 
    { 
     double currentIteration = pow(-1, n)*pow(y, 2*n)/fact(2*n); 
     sumCosine += currentIteration; 
    } 
    return sumCosine; 
} 

Итак, вот мой код. Я довольно доволен этим. За исключением одного: sineOfnumber и cosOfnumber после вызова sinTaylor и cosTaylor добавят друг друга в следующую строку cout, которая будет печатать друг друга. Другими словами, если число равно разрешению, то в строке, предназначенной для печати cosineOfnumber, будет напечатано 1.14, а в строке sineOfnumber будет напечатан результат. Может ли кто-нибудь помочь мне определить, почему это так? Спасибо огромное!

Answer 1

Вы когда-нибудь инициализировали переменные sumSine и sumCosine в своих функциях? Они не гарантированно начинаются с нуля, поэтому, когда вы вызываете + = внутри цикла, вы можете добавлять вычисленные значения в мусор.

Попробуйте инициализировать эти две переменные до нуля и посмотреть, что происходит, как и ничего, что код выглядит нормально.

Answer 2

Серия для синуса (извините за LaTeX):

sin(x) = \sum_{n \ge 0} \frac{x^{2 n + 1}}{(2 n + 1)!} 

Если вы посмотрите, данный термин t_ {2 п + 1} можно вычислить срок T_ {2 п + 3} как

t_{2 n + 3} = t_{2 n + 1} * \frac{x^2}{(2 n + 2)(2 n + 3)} 

Так что, учитывая термин, вы можете легко вычислить следующий. Если вы посмотрите на серию для косинуса, это похоже. Полученная программа более эффективна (без пересчета факториалов) и может быть более точным. При добавлении чисел с плавающей запятой точнее добавлять их от самых маленьких до самых больших, но я сомневаюсь, что это будет иметь значение здесь.