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

Решение уравнения методом Ньютона. Локализация корня - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Определить структуру ModelWindow для работы с моделями экранных окон http://www.cyberforum.ru/cpp-beginners/thread1050386.html
Определить структуру ModelWindow для работы с моделями экранных окон. В качестве полей задаются: заголовок окна, координаты левого верхнего угла, размер по горизонтали, размер по вертикали, состояние "видимое\невидимое", состояние "с рамкой\без рамки". Координаты и размеры указываются в целых числах. Кроме функций read и print реализовать следующие функции: передвижение окна по горизонтали, по...
C++ необходимо изменить программу #include <iostream> int main() { const int N = 5; using namespace std; double n; int count = N; float S = 0; std::cout << "vvedi 5 chisel\n"; http://www.cyberforum.ru/cpp-beginners/thread1050377.html
C++ ЗАполнение
struct person { char*name; char *adres; int age; }; void input (person *c, int kol) { char str; int p; for (int i =0;i<kol;i++) //Иду по массиву указателей
C++ Дан текстовый файл. Напечатать в алфавитном порядке
Дан текстовый файл. Напечатать в алфавитном порядке все слова из заданного файла, имеющие длину n.
C++ Строковые массивы http://www.cyberforum.ru/cpp-beginners/thread1050365.html
Здравствуйте! Подскажите пожалуйста, что есть такое строковые массивы в С++? Если можно с примером. Заранее благодарен :)
C++ В чем ошибка? #include <iostream.h> #include <time.h> #include <string.h> #include <stdio.h> #include <stdlib.h> int n, col_razr=3; int velich_razr(int chislo,int razr) { while(razr>1) подробнее

Показать сообщение отдельно
Chiakiii
0 / 0 / 0
Регистрация: 27.11.2011
Сообщений: 38

Решение уравнения методом Ньютона. Локализация корня - C++

22.12.2013, 03:19. Просмотров 1465. Ответов 10
Метки (Все метки)

Собственно, программа имеется, все довольно таки просто, но необходимо локализовать корень и найти именно положительный корень уравнения. Я так понял, нужно найти начальное приближение Xn и Xm. Приближение необходимо искать методом деления пополам. Как мне объяснили, необходимо выставить пределы для деления пополам от 0 до 100 и потом, когда метод будет сходится уже что-то принимать за Xn и Xm. Я в этом не силен. Необходима помощь, желательно программно.
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
62
63
64
65
66
67
68
69
70
71
72
73
#include <cstdlib>
#include <iostream>
#include <math.h>
#include <iomanip>
 
using namespace std;
double e = 0.001; 
 
 
double funct(double x)
{
    double f = x*x+4*sin(5*x)+exp(x+1)-3.2;
    return f;
}
 
 
double diff(double x)
{
    double f = 2*x+exp(x+1)+20*cos(5*x);
    return f;
}
 
 
int task_2()
{
    int n = 0;
    double xm = 0.95;
    double xn = 0;
    double xnprev;
    double acc = fabs(xn - xm);
    //int a=0,b=2;
    /*
    double c;
    while ((xn-xm)>=e)
          {
               c=(xn+xm)/2;
               if (funct(xm)*funct(c)<0)
               xn=c;
               else xm=c;                         
           }*/
    cout << "Решение методом Ньютона:" << endl << endl;
    while (acc > e)
    {
        n++;
        xnprev=xn;
        xn = xm - funct(xm)/diff(xm);
        if(xn<0){
            //xn=xnprev;
            //break;
        }
        acc = fabs(xn - xm);
        cout << setw(2) << n << " " << xm << " " << xn << " " << acc << endl;
        xm = xn;
    }
    cout << endl << "Корень уравнения x = " << xn << endl;
    return 0;
}
 
int main()
{
    
          
    setlocale(LC_ALL,"Rus");
    int ch;
   
    cout.precision(4);
    cout.setf(ios::fixed);
    
    task_2();
        
    system("PAUSE");
    return EXIT_SUCCESS;
}
Добавлено через 9 часов 38 минут
Метод деления пополам закомменчен в task_2()

Добавлено через 17 часов 36 минут
Кхм, не понимаю. Код дан, все работает. Мне нужно узнать буквально теоретический ответ на вопрос.
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 16:43. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru