Массив

    #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;
     double zb[zn];
     double zla[zn][zn];
   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<<"Вы выбрали метод Зейделя";
                cout << '\n'; 
 
                std::cout << "Введите NUM уравнений:";
                std::cin >>  zr;
            
         for (int zi = 0; zi < zr; zi++)
            for (int zj = 0; zj < zr; zj++)
            {{
                 std::cout << "Заполните матрицу A [" << zi+1 << "][" << zj+1 << "] = ";
                std::cin >> zla[zi][zj];
            
            }}
            for (int zi = 0; zi < zr; zi++) {
                std::cout << "Вектор B [" << zi << "] = ";
                std::cin >> zb[zi];
            }
 
            // Выводим матрица и вектор 
 
                cout<<("Матрица А\n");
    
        for (int i = 0; i<zr; i++)
        {
            for(int j = 0;j<zr; j++)
                cout<<zla[i][j]<<"\t";
            
            cout<<"\n";
        }
        
        printf("Вектор B\n");
           for (int i = 0; i<zr; i++){
        cout<<zb[i]<<"\n";
           }
 
        
            cout << '\n'; 
            
 
     for (int zi=0;zi<zr;zi++)
     {
      zx[zi]=0;
     }
 
     if (Zeydel (zla,zb,zx,zeps))
     {
         cout << '\n'; cout << '\n'; 
         std::cout<<"Результат";
         cout << '\n'; cout << '\n'; 
      for (int zi=0;zi<zr;zi++)
      {
       std::cout<<zx[zi]<< ' ';
      }}   
 
    
            
   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<<"Число интерации "<<zii;
       return 1;
     }

double *xz=new double[zn];
for(int i=0;i<zn;i++)
double **zla=new int*[zn];

double **zla
for(int i=0;i<zn;i++)
zla[i]=new int*[zn];

double**zla=new double *[zn];
for(int i=0;i<zn;i++)
zla[i]=new double[zn];