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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
amfisat
67 / 67 / 1
Регистрация: 16.06.2009
Сообщений: 240
#1

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

06.03.2011, 19:18. Просмотров 695. Ответов 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;
}
Спасибо!
0
Лучшие ответы (1)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.03.2011, 19:18
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Методом Ньютона решить линейное уравнение (C++):

Решить уравнение по заданной формуле методом хорд и модифицированным методом Ньютона - C++
8. Решить уравнение 5sinx - x = 0 на отрезке с точностью E методом хорд и модифицированным методом Ньютона.

Решить линейное уравнение ax = b - C++
Помогите пожалуйста решить задачу на си++ Решите линейное уравнение ax = b.

Решить линейное уравнение - C++
Решить линейное уравнение вида ax=b. Коэффициенты а и b читаются с экрана. Значение х вывести на экран с тремя знаками после запятой. Можно...

Решить линейное уравнение A*x + B = 0 - C++
Скоро экзамен , кто чем может ) Помоги ближнему , ибо смех продлевает жизнь...) 2. Решить линейное уравнение A*x + B = 0,...

Решить функцию методом Ньютона - C++
Нужно решить функцию методом Ньютона 0.1 x^2 - x \ln \,x=0. Отрезок содержащий корень . Точное значение должно выдать 1,1183. Но моя прога...

Уравнение методом Ньютона - C++
Доброго времени суток, помогите пожалуйста, справиться с уравнением, которое нужно реализовать методом Ньютона. Есть уравнение, скрин...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Евгений М.
1035 / 976 / 54
Регистрация: 28.02.2010
Сообщений: 2,829
Завершенные тесты: 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?
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
07.03.2011, 09:52
Привет! Вот еще темы с ответами:

Решить нелинейное уравнение методом хорд и методом простой итерации - C++
Help, please, послезавтра сдавать лабу, а я не знаю как делать=( Нужно решить нелинейное уравнение методом хорд и методом простой...

Решить уравнение методом Гаусса - C++
надо посчитать методом гаусса ур-е A \cdot X = B 3-размерность A = \begin{pmatrix} 2.74 &amp; -1.18 &amp; 3.17 \\ 1.12 &amp; 0.83 &amp;...

Решить уравнение методом итераций - C++
Здравствуйте,помогите пожалуйста разобраться,4 уравнение,с методом дихотомии разобрался,а с итерацией никак не выходит

Решить уравнение методом простых итераций - C++
Помогите решить уравнение х^4+2*x^3-x-1 на интервале #include &lt;iostream&gt; #include &lt;math.h&gt; using namespace std; //Функция ...


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

Или воспользуйтесь поиском по форуму:
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru