Форум программистов, компьютерный форум CyberForum.ru

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

Восстановить пароль Регистрация
 
Sinko
0 / 0 / 0
Регистрация: 27.10.2011
Сообщений: 7
27.10.2011, 19:06     Помогите найти ошибку #1
Проблема , по моему мнению, заключается в том что выдает неправильную полную сумму ряда.
Для ряда, члены которого вычисляются по формуле, подсчитать сумму членов ряда с точностью до 0.000001 и сумму первых 10 членов ряда. Если Вы считаете это необходимым, можете упростить или преобразовать выражение.
http://www.cyberforum.ru/attachment....1&d=1319727768

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
#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;
} /* конец программы */
Изображения
 
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
27.10.2011, 19:06     Помогите найти ошибку
Посмотрите здесь:

C++ Помогите найти ошибку
Помогите найти ошибку C++
C++ помогите найти ошибку.
C++ помогите найти ошибку
C++ Помогите найти ошибку
Помогите найти ошибку C++
Помогите найти ошибку C++
Помогите найти ошибку C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
-=ЮрА=-
Заблокирован
Автор FAQ
28.10.2011, 10:07     Помогите найти ошибку #2
Цитата Сообщение от Sinko Посмотреть сообщение
Для ряда, члены которого вычисляются по формуле, подсчитать сумму членов ряда с точностью до 0.000001 и сумму первых 10 членов ряда. Если Вы считаете это необходимым, можете упростить или преобразовать выражение.
C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
#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;
}
Миниатюры
Помогите найти ошибку  
-=ЮрА=-
Заблокирован
Автор FAQ
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 = -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.
Думаю после пояснений всё стало предельно ясно
Yandex
Объявления
28.10.2011, 10:18     Помогите найти ошибку
Ответ Создать тему
Опции темы

Текущее время: 15:49. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru