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

Метод хорд - C++

Восстановить пароль Регистрация
 
longcat
 Аватар для longcat
3 / 3 / 1
Регистрация: 15.10.2012
Сообщений: 18
16.06.2013, 14:11     Метод хорд #1
Всем привет! В общем, при решении курсовой по сабжу, столкнулся со следующей проблемой. Первое уравнение программа решает правильно, а второе — выводит на экран "nan". (Примечание: фиксурем точку b, т.к. f(b)*f''(b)>0 для обеих функций.) Собственно, весь код:
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
53
54
55
56
57
58
59
60
61
#include <iostream>
#include <clocale>
#include <conio.h>
#include <cmath>
using namespace std;
 
double f1(double x){
    return 10*x-exp(x);
}
double f2(double x){
    return tan(3.43*x)-5.07*x;
}
int bCheck (const int a, const int b){
    if(f1(a)*(-exp(a)>0)) return a;
        else return b;
}
int fCheck (const int a, const int b){
    if (f2(a)*(tan(3.47*a)*((2*3.47*3.47)/(cos(3.47*a)*cos(3.47*a))))>0)return a;
        else return b;
}
 
double aSolveFind (const int a, const int b, double eps, double (*f) (double x)){
    double next=b;
    double next1=a-((f(a)*(next-a))/(f(next)-f(a)));
    while (abs(next1-next)>=eps){
        next=next1;
        next1=a-((f(a)*(next-a))/(f(next)-f(a)));
    }
    return next1;
}
double bSolveFind (const int a, const int b, double eps, double (*f)(double x)){
    double next=a;
    double next1=next-((f(next)*(b-next))/(f(b)-f(next)));
    while(abs(next1-next)>=eps){
        next=next1;
        next1=next-((f(next)*(b-next))/(f(b)-f(next)));
    }
    return next1;
}
int main(){
    setlocale (0,"");
    double a,b;
    double eps=0.00001;
    int SwOp;
    cout<<"Выберите уравнение:\n";
    cout<<"1. f(x)=10x-e^x; a=3, b=4.\n";
    cout<<"2. f(x)=tg(3.43x)-5.07x; a=0.1, b=0.3\n";
    cout<<"Ваш выбор...........";
    cin>>SwOp;
    switch(SwOp){
        case 1:{a=3; b=4;
            if (bCheck(a,b)==b)
                cout<<"Корни данного уравнения на промежутке ["<<a<<","<<b<<"] : "<<bSolveFind(a,b,eps,f1)<<endl;
            else cout<<"Корни данного уравнения на промежутке ["<<a<<","<<b<<"] : "<<aSolveFind(a,b,eps,f1)<<endl; break;}
        case 2:{a=0.1; b=0.3;
            //if(fCheck(a,b)==a)
               // cout<<"Корни данного уравнения на промежутке ["<<a<<","<<b<<"] : "<<aSolveFind(a,b,eps,f2)<<endl;
                   /* else*/ cout<<"Корни данного уравнения на промежутке ["<<a<<","<<b<<"] : "<<bSolveFind(a,b,eps,f2)<<endl; break;}
    }
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
16.06.2013, 14:11     Метод хорд
Посмотрите здесь:

МЕТОД ХОРД И КАСАТЕЛЬНЫХ C++
метод хорд C++
Комбинированный метод хорд и касательных C++
C++ МЕтод хорд
метод хорд C++
C++ Метод хорд
C++ Метод хорд
C++ Метод хорд в С++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
longcat
 Аватар для longcat
3 / 3 / 1
Регистрация: 15.10.2012
Сообщений: 18
16.06.2013, 14:39  [ТС]     Метод хорд #2
Метод хорд
Вот... ┴──┴ ︵╰(°□°╰)
longcat
 Аватар для longcat
3 / 3 / 1
Регистрация: 15.10.2012
Сообщений: 18
16.06.2013, 16:35  [ТС]     Метод хорд #3
Upd: Разобрался, во 2 функции мы же передаём неправильную дробь...
Yandex
Объявления
16.06.2013, 16:35     Метод хорд
Ответ Создать тему
Опции темы

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