Помогите найти ошибку

@Sinko · Регистрация: 27.10.2011

Author24 — интернет-сервис помощи студентам

Проблема , по моему мнению, заключается в том что выдает неправильную полную сумму ряда.
Для ряда, члены которого вычисляются по формуле, подсчитать сумму членов ряда с точностью до 0.000001 и сумму первых 10 членов ряда. Если Вы считаете это необходимым, можете упростить или преобразовать выражение.
https://www.cyberforum.ru/atta... 1319727768

C++

#include <stdio.h>
#include <math.h>
#include <conio.h>
int main(void) { long n;              /* параметр ряда */
 double dbln;         /* параметр ряда в форме с плавающей точкой */
 double sum=0;        /* сумма членов ряда */
 double term;         /* значение текущего члена */
 const double eps=0.000001;  /* предел точности */
 long k2=1;           /* pow(2,n)*/
 short k1=1;          /* pow(-1,n)*/
 long k5=1;           /* pow(5,n)*/
   /* основной цикл; в модификациях вычисляются 
      следующие значения pow(2,n) и pow(-1,n)*/
   for (n=0; ; n++, k2*=2, k1=-k1, k5*=5) {
     /* преобразование n в форму с плавающей точкой */
     dbln=n;
     /* вычисление очередного члена */                          
     term=k1*(k2)/(dbln+k5);
     /* проверка достижения предела точности */       
     if (fabs(term)>=eps) 
        /* если не достигнут - накопление суммы */ 
        sum+=term;     /* если достигнут - выход из цикла */ 
     else break;      /* если 10 членов - вывод суммы */               
     if (n==9)       
        printf("summa 10 chlenov rada= %10.7lf\n",sum);     }
   /* конец основного цикла */
   /* вывод окончательной суммы */
   printf("polnay summa rada = %10.7lf\n",sum);     
    getch();
   return 0;
} /* конец программы */

@~~-=ЮрА=-~~ · 28.10.2011, 10:07

Сообщение от Sinko

Для ряда, члены которого вычисляются по формуле, подсчитать сумму членов ряда с точностью до 0.000001 и сумму первых 10 членов ряда. Если Вы считаете это необходимым, можете упростить или преобразовать выражение.

C

#include <stdio.h>
#include <conio.h>
 
double my_abs(double val)
{
    if(val < 0)
        val = -1*val;
    return val;
}
 
int main()
{
    double eps = 0;
    printf("Enter accuracy e : ");
    scanf("%lf",&eps);
    long m, n = 1;
    printf("Enter num of first elements : ");
    scanf("%u",&m);
    double pn   = 5;
    double sum  = 0;//Ñóììà ýëåìåíòîâ ðÿäà
    double summ = 0;//Ñóììà ïåðâûõ m ýëåìåíòîâ
    double an = (-1)*2.0/(1 + pn);
    while(eps <= my_abs(an))
    {
        n   += 1;
        sum += an;
        if(n <= m)
            summ += an;
        pn *= 5;
        an *= (-1)*2.0*(n - 1 + pn/5)/(n + pn);
    }
    printf("Sum of m-elements : %lf\n",summ);
    printf("Sum of sequence   : %lf\n",sum);
    printf("Press any key to continue\n");
    getch();
    return 0;
}

@~~-=ЮрА=-~~ · 28.10.2011, 10:18

Сообщение от -=ЮрА=-

while(eps <= my_abs(an))

Условиями окончания подобных задачь является fabs(f(x + dx) - f(x)) < e (я не буду вдаваться в дебри математике, можете смело поверить). Т.е в задаче нас интересуют две суммы SUM(a(1..n)) и
SUM(a(1..n - 1)), эти суммы отличаются на аn, т.е fabs(SUM(a(1..n)) - SUM(a(1..n - 1))) = fabs(an)
Это и будет условием стоп-итерационного процесса. В задаче решил обойтись без стандартной библиотеки math.h и ввёл свою функцию по расчёту модулю числа

Сообщение от -=ЮрА=-

double my_abs(double val)
{
* * * * if(val < 0)
* * * * * * * * val = -1*val;
* * * * return val;
}

- ну право подключать хедер всего из за одной простенькой функции это несерьёзно

Добавлено через 3 минуты

Сообщение от -=ЮрА=-

an *= (-1)*2.0*(n - 1 + pn/5)/(n + pn);

- а теперь это. Возьмём два смежных члена ряда an и an+1:
Найдём an+1/an
an = (-1)^n*2^n/(n + 5^n)
an+1 = (-1)^(n+1)*2^(n+1)/(n + 1 + 5^(n+1))

an+1/an = (-1)*2*(n + 5^n)/(n + 1 + 5^(n+1)) - вводим доп переменную отвечающую за 5^n

Сообщение от -=ЮрА=-

pn *= 5;

и свободно считаем элементы ряда без pow.
Думаю после пояснений всё стало предельно ясно

@~~-=ЮрА=-~~ Заблокирован
	28.10.2011, 10:18	3
	Сообщение от -=ЮрА=- while(eps <= my_abs(an)) Условиями окончания подобных задачь является fabs(f(x + dx) - f(x)) < e (я не буду вдаваться в дебри математике, можете смело поверить). Т.е в задаче нас интересуют две суммы SUM(a(1..n)) и SUM(a(1..n - 1)), эти суммы отличаются на аn, т.е fabs(SUM(a(1..n)) - SUM(a(1..n - 1))) = fabs(an) Это и будет условием стоп-итерационного процесса. В задаче решил обойтись без стандартной библиотеки math.h и ввёл свою функцию по расчёту модулю числа Сообщение от -=ЮрА=- double my_abs(double val) { * * * * if(val < 0) * * * * * * * * val = -1val; * * * return val; } - ну право подключать хедер всего из за одной простенькой функции это несерьёзно Добавлено через 3 минуты Сообщение от -=ЮрА=- an = (-1)2.0(n - 1 + pn/5)/(n + pn); - а теперь это. Возьмём два смежных члена ряда an и an+1: Найдём an+1/an an = (-1)^n2^n/(n + 5^n) an+1 = (-1)^(n+1)2^(n+1)/(n + 1 + 5^(n+1)) an+1/an = (-1)2(n + 5^n)/(n + 1 + 5^(n+1)) - вводим доп переменную отвечающую за 5^n Сообщение от -=ЮрА=-* pn *= 5; и свободно считаем элементы ряда без pow. Думаю после пояснений всё стало предельно ясно 0