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

Методом Ньютона решить линейное уравнение - C++

Восстановить пароль Регистрация
 
amfisat
 Аватар для amfisat
67 / 67 / 1
Регистрация: 16.06.2009
Сообщений: 235
06.03.2011, 19:18     Методом Ньютона решить линейное уравнение #1
Всем привет и хорошего настроения!
Решаю вычмат. Методом Ньютона надо решить линейное уравненьице. Сам метод запроган. Хочу запрогать процесс нахождения начального приближения x0.
Мои рассуждения: с помощью графика определяю границы промежутка, где сидят все корни уравнения. Далее иду с шагом step=0.1 от левого конца этого промежутка, пока не встречу такое х0, что f(x0)*fdx(x0)<=0, где fdx - производная функции. Функция здесь меняет знак=> на этом небольшом промежутке есть корень=> уменьшаю шаг :step*=0.1 и повторяю тоже самое, пока f(x0) по модулю не станет <=0.1 - после этого прекращаю вычисления и беру этот x0.
У меня всё зацикливается. Помогите, пжлст, понять, в чем проблема и советом.
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
#include<iostream>
#include<cmath>
using namespace std;
 
double fx (double n) // вычисление самой функции
{
    return (n*n*n-2.25*n*n-0.25*n+1.5);
}
double fdx (double n) // вычисление производной
{
    return (3*n*n*-2.25*2*n-0.25);
}
 
double X0(double step, int deg) // поиск нач. приближения, границы корней: [-1, 2.1]
{
    double x0=-1;
    while (deg!=0) // deg - степень полинома
    {
        if (fx(x0)*fdx(x0)<=0)
        {
            do
            {
                return X0(0.1*step,deg);
            }
            while (abs(fx(x0))>=0.1);
 
            if (abs(fx(x0))<0.1)
                return x0;
            
        }
        deg--;
        x0+=step;
    }
    
}
int main()
{
    double st=0.1; 
    cout<<X0(st, 3)<<endl;
    
    return 0;
}
Спасибо!
Лучшие ответы (1)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.03.2011, 19:18     Методом Ньютона решить линейное уравнение
Посмотрите здесь:

C++ Методом Эйлера решить дифференциальное уравнение
Решить уравнение x^3+8x+10=0 методом хорд и касательных C++
Решить уравнение методом простых итераций C++
Решить уравнение методом Гаусса C++
C++ Решить функцию методом Ньютона
Решить линейное уравнение A*x + B = 0 C++
Решить линейное уравнение ax = b C++
Решить линейное уравнение C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Евгений М.
1033 / 974 / 53
Регистрация: 28.02.2010
Сообщений: 2,817
Завершенные тесты: 2
07.03.2011, 09:52     Методом Ньютона решить линейное уравнение #2
Сообщение было отмечено автором темы, экспертом или модератором как ответ
amfisat, это:
C++
1
2
3
4
5
do
{
       return X0(0.1*step,deg);
}
while (abs(fx(x0))>=0.1);
и это:
C++
1
return X0(0.1*step,deg);
одно и тоже. Этот return завершает выполнение функции X0 и возвращает значение вызываемой функции.

Вычмат я не успел пройти в универе. Есть такой вопрос: почему вызывая функцию X0(...) всегда корни начинаются находится с x=-1? Почему не начинается с того значения x, при котором f(x)/f'(x)>0 и f(x+step)/f'(x+step)<0?
Yandex
Объявления
07.03.2011, 09:52     Методом Ньютона решить линейное уравнение
Ответ Создать тему
Опции темы

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