Метод Ньютона - программа уходит в бесконечный цикл

@dfg · Регистрация: 29.07.2014

Студворк — интернет-сервис помощи студентам

Подскажите где туплю, уже недели две над этой прогой сижу, а оно всё не работает хоть убей.... получается в бесконечный цикл заводит, потому что к корням никак не подойдёт...

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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
#include <iostream>
#include <cstring>
#include <cctype>
#include <conio.h>
#include <cmath>
 
using std::cout;
using std::cin;
using std::endl;
 
using namespace std;
 
struct member  //структура членов 
{
    char *name; //имя 
    double k; //коэффицент
    int st;//степень
};
 
void Input(member *s) // функция заполнения структуры
    {
        cout << "Введите имя переменной:";
        s -> name = new char;
        cin.get();
        cin.getline(s->name,2);
        cout << "Введите коэффициент:";
        cin >> s-> k;
        cout << "Введите степень:";
        cin >> s-> st;
        cout << endl;
    }
 
int main()
{
    setlocale(LC_ALL , "Russian");
    double e;
    cout << "Задайте относительную погрешность вычисления аргументов:" << endl;
    cin >> e;
    int N, M;
    cout << "Введите количество неизвестных:" << endl;
    cin >> M;
    cout << "Введите количество уравнений:" << endl;
    cin >> N;
 
    member **p = new member*[N];
    for(int i(0); i < N; i++)
        p[i] = new member[M];
    
    for(int i(0); i < N; i++)
    {
        for(int j(0); j < M; j++)
        {
            Input(&p[i][j]);
        }
        cout << endl;
    }
    
    double *t = new double[N];
    for(int i(0); i<N; i++)
    {
        cout << "Ввведите свободный член уравнения №" << i+1 << ':' << endl;
        cin >> t[i];
    }
 
    double *f = new double[M];
    for(int i(0); i<M; i++)
    {
        cout << "Введите значение начального приближения " << p[0][i].name << ':' << endl;
        cin >> f[i];
    }
    int n;
    //
        point:
 
 
    //
 
    double **jak = new double*[N]; // формирование матрицы Якоби
    for(int i(0); i < N; i++)
        jak[i] = new double[M];
 
    for(int i(0); i < N; i++)
    {
        for(int j(0); j < M; j++)
        {
                double c=0;
                for(int n(0); n < M; n++)
                        c=c+(p[i][n].k)*pow((float)f[n],(p[i][n].st));   
                double h=0;
                for(int n1(0); n1 < M; n1++)
                {
                    if (n1 == j)
                        h=h+(p[i][n1].k)*pow((float)(f[n1]+f[n1]*e),(p[i][n1].st));      
                    else
                        h=h+(p[i][n1].k)*pow((float)f[n1],(p[i][n1].st));
                }
            jak[i][j]=(h-c)/(f[j]*e);
        }
    }
 
    // Формирование расширенной матрицы системы линейных алгебраических уравнений:
 
    double **mass = new double*[N]; 
    for(int i(0); i < N; i++)
        mass[i] = new double[(M+1)];
 
    for(int i(0); i < N; i++)
    {
        for(int j(0); j < M; j++)
            mass[i][j]=jak[i][j];
    }
        for(int i(0); i < N; i++)
    {
            for(int j(M); j < M+1; j++)
            mass[i][j]=t[i];
    }
 
    // Решение системы линейных алгебраических уравнений методом Гауса:
    
    // Прямой ход:
    for(int i(0); i < N; i++) 
    {
        double a=mass[i][i];
        for(int j=i+1; j < N; j++)  
        {
            double b=mass[j][i];
            for(int k=i; k < M+1; k++) 
            mass[j][k]=mass[i][k]*b-mass[j][k]*a;   
        }
            
    }
 
    // Обратный ход:
    double *t1 = new double[M];
    for(int i=N-1; i>-1; i--)
    {
        double summ=0;
        for(int j=i+1; j<M+1; j++)
            summ=summ+mass[i][j]*t1[j];
        summ=mass[i][M]-summ;
        if(mass[i][i]==0)
            cout << "Нет решения" << endl;
        t1[i]=summ/mass[i][i];
    }
    //Проверка условия:
 
    n=n+1;
        for(int i(0); i < M; i++)
        if (fabs(t1[i]/f[i])>e || n>20)
        {
            for(int i(0); i<N; i++)
            cout << p[0][i].name << '=' << t1[i] << endl;
            cout << "Количество итераций:" << n << endl;
        }
        else
        {
            for(int j(0); j < M; j++)
                f[j]=t1[j];     
            goto point;
        }
        
    for(int i(0); i < N; i++)
        delete [] p[i];
    delete [] p;
 
    for(int i(0); i < N; i++)
        delete [] p[i];
    delete [] p;
 
    for(int i(0); i < N; i++)
        delete [] p[i];
    delete [] p;
    
    _getch();
    return 0;
}

@Donnerwetter · 17.08.2014, 21:49

C++
1
n=n+1;

переменная не инициализирована, какое значение инкрементируется?

C++
1
2
3
4
5
6
7
8
9
10
11
for (int i(0); i < N; i++)
    delete[] p[i];
delete[] p;
 
for (int i(0); i < N; i++)
    delete[] p[i];
delete[] p;
 
for (int i(0); i < N; i++)
    delete[] p[i];
delete[] p;

зачем эти три одинаковые действия? память удаляется уже в первом случае, а потом идет обращение к несуществующему участку

@-=ЮрА=- · 17.08.2014, 22:24

dfg, Методом Гаусса решить систему n линейных алгебраических уравнений

@dfg · 18.08.2014, 00:56 **[ТС]**

ну вот так вроде правильно, но сходится к ответу всё-равно не хочет

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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
#include <iostream>
#include <cstring>
#include <cctype>
#include <conio.h>
#include <cmath>
 
using std::cout;
using std::cin;
using std::endl;
 
using namespace std;
 
struct member  //структура членов 
{
    char *name; //имя 
    double k; //коэффицент
    int st;//степень
};
 
void Input(member *s) // функция заполнения структуры
    {
        cout << "Введите имя переменной:";
        s -> name = new char;
        cin.get();
        cin.getline(s->name,2);
        cout << "Введите коэффициент:";
        cin >> s-> k;
        cout << "Введите степень:";
        cin >> s-> st;
        cout << endl;
    }
 
int main()
{
    setlocale(LC_ALL , "Russian");
    double e;
    cout << "Задайте относительную погрешность вычисления аргументов:" << endl;
    cin >> e;
    int N, M;
    cout << "Введите количество неизвестных:" << endl;
    cin >> M;
    cout << "Введите количество уравнений:" << endl;
    cin >> N;
 
    member **p = new member*[N];
    for(int i(0); i < N; i++)
        p[i] = new member[M];
    
    for(int i(0); i < N; i++)
    {
        for(int j(0); j < M; j++)
        {
            Input(&p[i][j]);
        }
        cout << endl;
    }
    
    double *t = new double[N];
    for(int i(0); i<N; i++)
    {
        cout << "Ввведите свободный член уравнения №" << i+1 << ':' << endl;
        cin >> t[i];
    }
 
    double *f = new double[M];
    for(int i(0); i<M; i++)
    {
        cout << "Введите значение начального приближения " << p[0][i].name << ':' << endl;
        cin >> f[i];
    }
    int n=0;
    //
        point:
 
 
    //
 
    double **jak = new double*[N]; // формирование матрицы Якоби
    for(int i(0); i < N; i++)
        jak[i] = new double[M];
 
    for(int i(0); i < N; i++)
    {
        for(int j(0); j < M; j++)
        {
                double c=0;
                for(int n(0); n < M; n++)
                        c=c+(p[i][n].k)*pow((float)f[n],(p[i][n].st));   
                double h=0;
                for(int n1(0); n1 < M; n1++)
                {
                    if (n1 == j)
                        h=h+(p[i][n1].k)*pow((float)(f[n1]+f[n1]*e),(p[i][n1].st));      
                    else
                        h=h+(p[i][n1].k)*pow((float)f[n1],(p[i][n1].st));
                }
            jak[i][j]=(h-c)/(f[j]*e);
        }
    }
 
    // Формирование расширенной матрицы системы линейных алгебраических уравнений:
 
    double **mass = new double*[N]; 
    for(int i(0); i < N; i++)
        mass[i] = new double[(M+1)];
 
    for(int i(0); i < N; i++)
    {
        for(int j(0); j < M; j++)
            mass[i][j]=jak[i][j];
    }
        for(int i(0); i < N; i++)
    {
            for(int j(M); j < M+1; j++)
            mass[i][j]=t[i];
    }
 
    // Решение системы линейных алгебраических уравнений методом Гауса:
    
    // Прямой ход:
    for(int i(0); i < N; i++) 
    {
        double a=mass[i][i];
        for(int j=i+1; j < N; j++)  
        {
            double b=mass[j][i];
            for(int k=i; k < M+1; k++) 
            mass[j][k]=mass[i][k]*b-mass[j][k]*a;   
        }
            
    }
 
    // Обратный ход:
    double *t1 = new double[M];
    for(int i=N-1; i>-1; i--)
    {
        double summ=0;
        for(int j=i+1; j<M+1; j++)
            summ=summ+mass[i][j]*t1[j];
        summ=mass[i][M]-summ;
        if(mass[i][i]==0)
            cout << "Нет решения" << endl;
        t1[i]=summ/mass[i][i];
    }
    //Проверка условия:
 
    n=n+1;
        //for(int i(0); i < M; i++)
        if (fabs(t1[1]/f[1])<e || n>=100)
        {
            for(int i(0); i<N; i++)
            cout << p[0][i].name << '=' << t1[i] << endl;
            cout << "Количество итераций:" << n << endl;
        }
        else
        {
            for(int j(0); j < M; j++)
                f[j]=t1[j];     
            goto point;
        }
        
    for(int i(0); i < N; i++)
        delete [] p[i];
    delete [] p;
    
    _getch();
    return 0;
}

Новые блоги и статьи Все статьи Все блоги /
Загрузка PNG-файла с альфа-каналом с помощью библиотеки SDL3_image на Android 8Observer8 27.01.2026 Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .	влияние грибов на сукцессию anaschu 26.01.2026 Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .	Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .
Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .	Использование SDL3-callbacks вместо функции main() на Android, Desktop и WebAssembly 8Observer8 24.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .	моя боль iceja 24.01.2026 Выложила интерполяцию кубическими сплайнами www. iceja. net REST сервисы временно не работают, только через Web. Написала за 56 рабочих часов этот сайт с нуля. При помощи perplexity. ai PRO , при. . .	Модель сукцессии микоризы anaschu 24.01.2026 Решили писать научную статью с неким РОманом

@-=ЮрА=- 6614 / 4256 / 401 Регистрация: 08.08.2009 Сообщений: 10,325 Записей в блоге: 24
	17.08.2014, 22:24
	dfg, Методом Гаусса решить систему n линейных алгебраических уравнений 0