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

найти ошибку. циклы - C++

Восстановить пароль Регистрация
 
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
12.02.2011, 13:11     найти ошибку. циклы #1
вычислить значение, соответствующее значению функции для указанного значения x0 (при указанных ограничениях на аргумент функции) с точностью ε (например, ε=0.001). Точность вычисления считается выполненной, если последний по модулю член в сумме меньше ε. Полученное значение сравнить с истинным значением данной функции, найденное при помощи математических встроенных функций, найдя абсолютную погрешность вычисления. Указать количество просуммированных членов такого ряда.

смотреть картинку.
мой код.
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
38
39
40
41
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <clocale>
using namespace std;
 
int main(int argc, char *argv[])
{
        setlocale( LC_ALL, ".1251" );
    int n;
    double e,delta,f,a,s,x,k;
    cout << "Введите переменную удовлетворяющую условию |x|<1: " << endl;
    cin >> x;
    cout << "Введите точность вычисления: " << endl;
    cin >> e;
     a = 1;
     s = 1;
     n = 1;
    f = (3*(exp(log(1+x)/3)))-3;
    cout << "Истинное значение f(" << x << ")= " << f << endl;
    do
    {
         k =(2*5/6*9)*pow(x,3);
         if (n % 2 !=0)         
  k = k*(-1);
  a =a*((x-(2/6)*pow(x,2))+k-((2*5*(3*n-4))/(6*9*(3*n))*pow(x,n)));
  s =s+a;
  n =n+1;
  
    } while (a<e);
    delta = f-s;
    cout << " При функции f(" << x << ")= " << s << endl;
    cout << " Погрешность = " << delta << endl;
    cout << " Количество просуммированых членов = " << n << endl;
    
    
    
    
    system("PAUSE");
    return EXIT_SUCCESS;
}

код моего друга на паскале
Pascal
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
program project1;
var
  n:integer;
  e,delta,f,a,s,x,k:double;
begin
 writeln('Vvedite x ydovletvor ysloviy |x|<1');
 read(x);
 writeln('Vvedite tochnost vychisleniy');
 read(e);
 f:=(3*(exp(ln(1+x)/3)))-3;
 writeln('Istinnoe znacheniy f(',x,')=',f);
 a:=1;
 s:=1;
 n:=1;
 repeat
  k:=(2*5/6*9)*sqr(x)*x;
  if n mod 2<>0 then
   k:=k*(-1);
  a:=a*((x-(2/6)*sqr(x))+(k)-((2*5*(3*n-4))/(6*9*(3*n))*exp(n*ln(x))));
  s:=s+a;
  n:=n+1;
 until a<e ;
 delta:=f-s;
 writeln('Pri function f(',x,')= ',s);
 writeln('Pogreshnost = ',delta);
 writeln('Kolochestvo prosumirovannyh chlenov = ',n);
 readln;
 readln;
end.
Миниатюры
найти ошибку. циклы  
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Jupiter
Каратель
Эксперт C++
6543 / 3963 / 226
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
12.02.2011, 14:25     найти ошибку. циклы #2
И в чем проблема? В паскале для возведения числа в степень используют логарифмическое тождество, потому ваш друг использовал exp, в С/С++ для этого есть функция pow
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
12.02.2011, 14:56  [ТС]     найти ошибку. циклы #3
проблема в неправильном вычислении, либо у него либо у меня либо у нас обоих
Jupiter
Каратель
Эксперт C++
6543 / 3963 / 226
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
12.02.2011, 15:23     найти ошибку. циклы #4
Да у вас ошибка в вычислениях , внутри цикла нужно считать общий член ряда потом проверять на точность и если он подходит по точности суммировать его а если нет то выходит из цикла. Вот что у меня получилось, вариант не идеален так как много повторяющихся операций, в идеале нужно испльзовать массив/вектор для взятия предыдущего члена
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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <clocale>
using namespace std;
 
int fact( int n )
{
    return n == 0 ? 1 : n * fact(n-1);
}
 
double P( double pow, int n )
{
    double res(1);
    for( int i = 0; i < n; i++ )
    {
        res *= pow - i;
    }
    return res;
}
 
int main()
{
    setlocale( LC_ALL, ".1251" );
    
    cout << "Введите переменную удовлетворяющую условию |x|<1: " << endl;
    double x;
    cin >> x;
    double e;
    cout << "Введите точность вычисления: " << endl;
    cin >> e;
    double f = (3*(exp(log(1+x)/3)))-3;
    cout << "Истинное значение f(" << x << ")= " << f << endl;
    double member(0), total_sum(0);
    int n(1);
    while(true)
    {
        member = 3 * P(1./3, n) * pow( x, n ) / fact(n);
        if( abs(member) > e )
        {
            total_sum += member;
            n++;
        }
        else break;
    }
    double delta = f-total_sum;
    cout << " При функции f(" << x << ")= " << total_sum << endl;
    cout << " Погрешность = " << delta << endl;
    cout << " Количество просуммированых членов = " << n << endl;
    system("PAUSE");
    return EXIT_SUCCESS;
}
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
12.02.2011, 15:38  [ТС]     найти ошибку. циклы #5
Цитата Сообщение от Maxwe11 Посмотреть сообщение
нужно испльзовать массив/вектор для взятия предыдущего члена
цель лабы научится использывать циклы, все остальное не использывать...
разбираю что вы посоветывали

Добавлено через 8 минут
ваша программа рассадник вопросов для меня
ну для чего там пи?
верхняя часть для чего, а конкретно фактр?

можно прокомментировать как то строки?
слишком все сложно... не помню для чего "?"
эти возрасщяющие значения return вообще пугают, нельзя ли проще?
Jupiter
Каратель
Эксперт C++
6543 / 3963 / 226
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
12.02.2011, 16:00     найти ошибку. циклы #6
Данный ряд является биномиальным рядом, общий член которого
http://www.cyberforum.ru/cgi-bin/latex.cgi?{\left(1 + x \right)}^{\alpha } = 1 + \sum_{n=1}^{\propto }\frac{\alpha *\left(\alpha - 1 \right)*{}_{....}*\left(\alpha  - n + 1 \right)}{n!}*{x}^{n}
fact - рекурсивное вычисление факториала, надеюсь знаете что такое рекурсия?
C++
1
2
3
4
return n == 0 ? 1 : n * fact(n-1);
// это эквивалентно
if(n==0) return 1;
else return fact(n-1)
P - вычисление произведения http://www.cyberforum.ru/cgi-bin/latex.cgi?\alpha *\left(\alpha - 1 \right)*{}_{....}*\left(\alpha  - n + 1 \right)

Не по теме:

PS если не секрет где и на каком курсе учитесь? а то странно - здание вам задали, а по вышке судя по всему ряды еще не изучали

IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
12.02.2011, 16:47  [ТС]     найти ошибку. циклы #7
1 курс 2 семестр, только началось программирование, в первом была информатика.
лабы дали - еще лекций не было) что умею - написал.
Учусь в Белгородском технологе. Техническая кибернетика.
Я один с потока с++ учу, тяжеловато. (остальной народ паскаль-делфи)

по теме:
ткните носом в строку где в моем коде начался косяк.
Jupiter
Каратель
Эксперт C++
6543 / 3963 / 226
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
12.02.2011, 17:19     найти ошибку. циклы #8
21 - 30
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
12.02.2011, 18:19  [ТС]     найти ошибку. циклы #9
как поправить этот участок, при том сохранить ду - вайл?
Jupiter
Каратель
Эксперт C++
6543 / 3963 / 226
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
12.02.2011, 18:25     найти ошибку. циклы #10
в моем коде поменять
C++
1
while(true){...}
на
C++
1
do{...}while(true)
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
13.02.2011, 11:49  [ТС]     найти ошибку. циклы #11
Цитата Сообщение от Maxwe11 Посмотреть сообщение
в моем коде поменять
а в моем)))
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
19.02.2011, 18:49  [ТС]     найти ошибку. циклы #12
Не правильный код, полная неразбериха с вычислениями..
Помогите с правкой..

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
38
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <clocale>
using namespace std;
 
int main(int argc, char *argv[])
{
        setlocale( LC_ALL, ".1251" );
    int n;
    double e,delta,f,a,s,x,k;
    cout << "Введите переменную удовлетворяющую условию |x|<1: " << endl;
    cin >> x;
    cout << "Введите точность вычисления: " << endl;
    cin >> e;
     a = 1;
     s = 1;
     n = 1;
    f = (3*(exp(log(1+x)/3)))-3;
    cout << "Истинное значение f(" << x << ")= " << f << endl;
    do
    {
         k =(2*5/6*9)*pow(x,3);
         if (n % 2 !=0)         
  k = k*(-1);
  a =a*((x-(2/6)*pow(x,2))+k-((2*5*(3*n-4))/(6*9*(3*n))*pow(x,n)));
  s =s+a;
  n =n+1;
  
    } while (a<e);
    delta = f-s;
    cout << " При функции f(" << x << ")= " << s << endl;
    cout << " Погрешность = " << delta << endl;
    cout << " Количество просуммированых членов = " << n << endl;
 
    system("PAUSE");
    return EXIT_SUCCESS;
}
мой альтернативный код:

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
38
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <clocale>
using namespace std;
 
int main(int argc, char *argv[])
{
        setlocale( LC_ALL, ".1251" );
    int n;
    double e,delta,f,a,s,x,k;
    cout << "Введите переменную удовлетворяющую условию |x|<1: " << endl;
    cin >> x;
    cout << "Введите точность вычисления: " << endl;
    cin >> e;
     a = 1;
     s = 1;
     n = 1;
    f = (3*(exp(log(1+x)/3)))-3;
    cout << "Истинное значение f(" << x << ")= " << f << endl;
    do
    {
         k =(2*5/6*9)*pow(x,3);
         if (n % 2 !=0)         
  k = k*(-1);
  a =a*((x-(2/6)*pow(x,2))+k-((2*5*(3*n-4))/(6*9*(3*n))*pow(x,n)));
  s =s+a;
  n =n+1;
  
    } while (a<e);
    delta = abs(f-s);
    cout << " При функции f(" << x << ")= " << s << endl;
    cout << " Погрешность = " << delta << endl;
    cout << " Количество просуммированых членов = " << n << endl;
 
    system("PAUSE");
    return EXIT_SUCCESS;
}
timchuchok
15 / 15 / 0
Регистрация: 21.12.2010
Сообщений: 55
20.02.2011, 00:44     найти ошибку. циклы #13
Тебе в программе можно обойтись без функции pow, ведь у тебя каждый следующий член равен предадущему умноженому на коефициэнт и икс ! ! !
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
20.02.2011, 00:56  [ТС]     найти ошибку. циклы #14
Цитата Сообщение от timchuchok Посмотреть сообщение
можно обойтись
как бы да... но смысл программы не правильный не из за этого....
Jupiter
Каратель
Эксперт C++
6543 / 3963 / 226
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
20.02.2011, 02:20     найти ошибку. циклы #15
IIIa66uMEM6eP, откройте учебник по высшей математике и почитайте 1-ю главу по теме ряды
timchuchok
15 / 15 / 0
Регистрация: 21.12.2010
Сообщений: 55
20.02.2011, 03:33     найти ошибку. циклы #16
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
double f( double  x) 
{
    double  e = 0.000001, sum, a;
    int n  = 2;
    a = x;
    sum = x;
    while ( a > e )
    {
        a *= -1.*(x * ( 3 * n - 4)) / (3 * n);
        sum += a;
        n++;
    }
    return sum;
}
это для 0< х <1
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.02.2011, 10:37     найти ошибку. циклы
Еще ссылки по теме:

C++ Найти сумму первых n-четных натуральных чисел (найти ошибку)
Циклы с параметром: найти сумму ряда C++
C++ Циклы. Найти сумму ряда

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

Или воспользуйтесь поиском по форуму:
IIIa66uMEM6eP
заставил Бендера
 Аватар для IIIa66uMEM6eP
432 / 288 / 10
Регистрация: 05.12.2010
Сообщений: 1,642
Записей в блоге: 6
20.02.2011, 10:37  [ТС]     найти ошибку. циклы #17
Цитата Сообщение от Maxwe11 Посмотреть сообщение
откройте учебник
хм... и вас в коде ошибка и у меня, но учебник должен только я открывать - странно
Yandex
Объявления
20.02.2011, 10:37     найти ошибку. циклы
Ответ Создать тему
Опции темы

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