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

Ошибка при использовании метода хорд - C++

Восстановить пароль Регистрация
 
korsarqa
 Аватар для korsarqa
7 / 7 / 0
Регистрация: 06.12.2011
Сообщений: 170
25.11.2012, 18:16     Ошибка при использовании метода хорд #1
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
// chord method.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include "conio.h"
#include <math.h>
#include "iostream"
#define e 2.71
 
//метод хорд
using namespace std;
 
 
double func(double x)
{
    return pow(e,sqrt(x))-1/x;
}
 
double find(double infinum, double supremum, double epsilon)
{
    while (fabs(supremum - infinum) > epsilon)
    {
        infinum = supremum - (supremum - infinum) * func(supremum) / (func(
                supremum) - func(infinum));
        supremum = infinum - (infinum - supremum) * func(infinum) / (func(
                infinum) - func(supremum));
    }
 
    return supremum;
}
 
int main()
{
    double a = -5, b = 5;
    cout << find(a, b, 0.001)<<endl;
    getch();
    return 0;
}
Использовал код, подставил свое уравнение, и ответ нормальным не получается...
Выходит -1.#IND

Что означает данная надпись? В задании у меня говорится что отделить корни выполнить надо самостоятельно
Мб мне корень из x заменить просто на x? И в ответе потом сделать обратное для нормального ответа

Добавлено через 4 минуты
Заменил sqrt(x) на просто x. Ответ так же не получается =(

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
// chord method.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include "conio.h"
#include <math.h>
#include "iostream"
#define e 2.71
 
//метод хорд
using namespace std;
 
 
double func(double x)
{
    //return pow(e,sqrt(x))-1/x;
    return pow(e,x)-1/(x*x);
}
 
double find(double infinum, double supremum, double epsilon)
{
    while (fabs(supremum - infinum) > epsilon)
    {
        infinum = supremum - (supremum - infinum) * func(supremum) / (func(
                supremum) - func(infinum));
        supremum = infinum - (infinum - supremum) * func(infinum) / (func(
                infinum) - func(supremum));
        supremum=sqrt(supremum);
    }
 
    return supremum;
}
 
int main()
{
    double a = -5, b = 5;
    cout << find(a, b, 0.001)<<endl;
    getch();
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
25.11.2012, 18:16     Ошибка при использовании метода хорд
Посмотрите здесь:

C++ ошибка при использовании delete
Ошибка при использовании конструктора C++ Builder
Ошибка при использовании STL C++
Ошибка при использовании вектора C++
Ошибка при использовании new[] / delete[] C++
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Kuzia domovenok
 Аватар для Kuzia domovenok
1882 / 1737 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
25.11.2012, 18:46     Ошибка при использовании метода хорд #2
В определённый момент
supremum становится очень близок infinum это недопустимо! Нельзя уменьшать супремум, увеличивать инфинум
Ведь тогда
func(supremum) == func(infinum) или почти равен
Следовательно делить на их разность нельзя
korsarqa
 Аватар для korsarqa
7 / 7 / 0
Регистрация: 06.12.2011
Сообщений: 170
25.11.2012, 18:59  [ТС]     Ошибка при использовании метода хорд #3
Цитата Сообщение от Kuzia domovenok Посмотреть сообщение
В определённый момент
func(supremum) == func(infinum) или почти равен
Следовательно делить на их разность нельзя
Тогда надо поставить в цикле в конце условие ?
if func(supremum) == func(infinum) break?

Добавлено через 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
38
39
#include <iostream>
#include <math.h>
#include <conio.h>
#include <stdio.h>
#include <conio.h>
#define e 2.71
 
using namespace std;
 
void main() {
    double a, b, x, eps, fc, fx, c; //a,b - концы отрезка, x - корень, e - эпсилон, fx/fc - значения функции в точках x и c
    int count=0;
    a=1;    
    b=2;
 
    eps=0.001;
 
    x=b;
    c=a;
 
    while(1) {
        count++;
        
        fc=pow(e,sqrt(x))-1/x;
        fx=pow(e,sqrt(x))-1/x;
 
 
        x=x-fx*(x-c)/(fx-fc);
 
        cout << "x: " << x << "\n";
        
        if(abs(fx)<eps) 
            break;
 
    }
    cout << "\nSteps: " << count;
    _getch();
 
}
Kuzia domovenok
 Аватар для Kuzia domovenok
1882 / 1737 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
25.11.2012, 19:04     Ошибка при использовании метода хорд #4
а правильный ответ какой?

Добавлено через 3 минуты
Цитата Сообщение от korsarqa Посмотреть сообщение
может как то преобразовать уравнение надо?
Дык, конечно надо!! зачем ты вообще в области отрицательных икс корень ищешь!?!? У тебя же корень из икс в уравнении стоит? А ты от МИНУС пяти до пяти работаешь.
korsarqa
 Аватар для korsarqa
7 / 7 / 0
Регистрация: 06.12.2011
Сообщений: 170
25.11.2012, 19:08  [ТС]     Ошибка при использовании метода хорд #5
Цитата Сообщение от Kuzia domovenok Посмотреть сообщение
а правильный ответ какой?

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

Дык, конечно надо!! зачем ты вообще в области отрицательных икс корень ищешь!?!? У тебя же корень из икс в уравнении стоит? А ты от МИНУС пяти до пяти работаешь.
точно не заметил, ty
Yandex
Объявления
25.11.2012, 19:08     Ошибка при использовании метода хорд
Ответ Создать тему
Опции темы

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