Вывод вектора ax+b

@apachan · Регистрация: 07.12.2010

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

Программа решает систему методом зейделя, помогите пожалуйста в завершающей точке вывести вектор ах-б

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
//#include "stdafx.h"
    #include <iostream>//i/o
    #include <conio.h> //getch
    #include <cmath>   //fabs, abs
  //  #include "stdafx.h"
    #include <math.h>
    #include <tchar.h> 
 
using namespace std;
 
 
 
     const int zn=2000;
     double zeps = 0.01;
     int zi;
     double zx[zn];
     int zr=4, N=2;
     //double zb[zn];
     //double zla[zn][zn];
     double zla[zn][zn]={{25,2,-2,1},{0,-20-0.3*N,3,0.1*N},{2,5,25,-2},{-4,-0.1*N,1,16+0.3*N}},
zb[zn]={{25*N+11},{-49-0.5*N},{2*N-67},{-3.9*N+13}};
   int vibs;
    
   double Zeydel (double zla[zn][zn], double zb[zn], double zx[zn], double ze);
    
 
    int main () {
    char vibormetoda;
    setlocale (LC_ALL, "Russian"); 
    
                cout<<"Zeidel:"<<endl;
                cout << '\n'; 
 
               // std::cout << "Aaaaeoa NUM o?aaiaiee:";
                //std::cin >>  zr;
            
        // for (int zi = 0; zi < zr; zi++)
          //  for (int zj = 0; zj < zr; zj++)
           // {{
              //   std::cout << "Caiieieoa iao?eoo A [" << zi+1 << "][" << zj+1 << "] = ";
             //   std::cin >> zla[zi][zj];
            
        //    }}
            //for ( zi = 0; zi < zr; zi++) {
             //   std::cout << "Aaeoi? B [" << zi << "] = ";
               // std::cin >> zb[zi];
          //  }
 
            // Auaiaei iao?eoa e aaeoi? 
 
                cout<<("Matrica A:\n")<<endl;
    
        for (int i = 0; i<zr; i++)
        {
            for(int j = 0;j<zr; j++)
                cout<<zla[i][j]<<"\t";
            
            cout<<"\n"<<endl;
        }
        
        printf("Matrica B:\n\n");
           for ( i = 0; i<zr; i++){
        cout<<zb[i]<<"\n";
           }
 
        
            cout << '\n'<<endl; 
            
 
     for ( zi=0;zi<zr;zi++)
     {
      zx[zi]=0;
     }
 
     if (Zeydel (zla,zb,zx,zeps))
     {
         cout << '\n'; cout << '\n'; 
         std::cout<<"Rezultat:";
         cout << '\n'; cout << '\n'; 
      for (int zi=0;zi<zr;zi++)
      {
       std::cout<<zx[zi]<<endl;
      
      }
     cout<<endl;}   
 
    
            
   system ("Pause");
    
 
        }
        
 
    double Zeydel (double zla[zn][zn], double zb[zn], double zx[zn], double ze)
    {
     int zi,zj;
     double zs1, zs2, zs, zv, zm;
    double zii=0;
     do
     {
      zm = 0;
      for (zi=0; zi< zr; zi++)
      {
       zs1=0;
       zs2=0;
 
       for (zj=0;zj<zi;zj++)  zs1+= zla[zi][zj]*zx[zj];
       for (zj=zi;zj<zr;zj++)  zs2+= zla[zi][zj]*zx[zj];
        zv=zx[zi];
        zx[zi]-=(zs1+zs2-zb[zi])/zla[zi][zi];
 
        if (fabs(zv-zx[zi])>zm)
         zm=fabs(zv-zx[zi]);
       }
        zii++;
       }
      while (zm>=ze);
 
        std::cout<<"Chislo Iteracii: "<<zii;
       return 1;
     }

Добавлено через 9 минут
в общем, перемножить матрицу а на х вычесть из этого матрицу б

@apachan · 03.06.2012, 17:37 **[ТС]**

Вот до чего дошел, выдает только последнее значение вместо всего столбца ах-б почему-то

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
    #include <iostream>//i/o
    #include <conio.h> //getch
    #include <cmath>   //fabs, abs
    #include <math.h>
    #include <tchar.h> 
 
using namespace std;
 
 
 
     const int zn=2000, q1=1000;
     double zeps = 0.01;
     int zi;
     double zx[zn];
     int zr=4, N=2;
     double re[q1];
     double zla[zn][zn]={{25,2,-2,1},{0,-20-0.3*N,3,0.1*N},{2,5,25,-2},{-4,-0.1*N,1,16+0.3*N}},
zb[zn]={{25*N+11},{-49-0.5*N},{2*N-67},{-3.9*N+13}};
   int vibs;
    
   double Zeydel (double zla[zn][zn], double zb[zn], double zx[zn], double ze);
    
 
    int main () {
    char vibormetoda;
    setlocale (LC_ALL, "Russian"); 
    
                cout<<"Zeidel:"<<endl;
                cout << '\n';
 
                cout<<("Matrica A:\n")<<endl;
    
        for (int i = 0; i<zr; i++)
        {
            for(int j = 0;j<zr; j++)
                cout<<zla[i][j]<<"\t";
            
            cout<<"\n"<<endl;
        }
        
        printf("Matrica B:\n\n");
           for ( i = 0; i<zr; i++){
        cout<<zb[i]<<"\n";
           }
 
        
            cout << '\n'<<endl; 
            
 
     for ( zi=0;zi<zr;zi++)
     {
      zx[zi]=0;
     }
 
     if (Zeydel (zla,zb,zx,zeps))
     {
         cout << '\n'; cout << '\n'; 
         std::cout<<"Rezultat:";
         cout << '\n'; cout << '\n'; 
      for (int zi=0;zi<zr;zi++)
      {
       std::cout<<zx[zi]<<endl;
      
      }
     cout<<endl;}   
 
    
            
   system ("Pause");
    
 
        }
        
 
    double Zeydel (double zla[zn][zn], double zb[zn], double zx[zn], double ze)
    {
     int zi,zj;
     double zs1, zs2, zs, zv, zm;
    double zii=0;
     do
     {
      zm = 0;
      for (zi=0; zi< zr; zi++)
      {
       zs1=0;
       zs2=0;
 
       for (zj=0;zj<zi;zj++)  zs1+= zla[zi][zj]*zx[zj];
       for (zj=zi;zj<zr;zj++)  zs2+= zla[zi][zj]*zx[zj];
        zv=zx[zi];
        zx[zi]-=(zs1+zs2-zb[zi])/zla[zi][zi];
 
        if (fabs(zv-zx[zi])>zm)
         zm=fabs(zv-zx[zi]);
       }
        zii++;
       }
      while (zm>=ze);
 
        std::cout<<"Chislo Iteracii: "<<zii<<endl;
 
        
    
        for (int i=0; i<4; i++)
{
for (int j=0; j<4; j++)
{
    re[q1]=zla[i][j]*zx[zi]-zb[i];}}
 
            for ( i=0; i<4; i++)
{
for (int j=0; j<4; j++)
{
    cout<<re[q1]<<endl;}}
 
       return 1;
 
       
    
     }

Новые блоги и статьи Все статьи Все блоги /
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Access VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .