Решение СЛАУ С++

[Foralberg]

Программа запускается, но после того, как я ввожу данные она зависает. Помогите плиз исправить это.

Задание:

Решить СЛАУ A^2*X=Y^T (транспонированная), где
A = ∣2 1 5 2|
|5 2 2 6|
|2 2 1 2|
|1 3 3 1∣ ,

Y = [ 3 1 2 1 ]

Вот моя программа:

#include "stdafx.h"
#include <iostream.h>
#include <math.h>

void tr_matr(double **mtr1, int n)
{
	for(int i=0; i<n; i++)
	{
		for(int j=i+1; j<n; j++)
		{
			double p=mtr1[i][j];
			mtr1[i][j]=mtr1[j][i];
			mtr1[j][i]=p;
		}
	}
}

void mult_matr(int n, double **mtr1, double **mtr2, double **res)
{
	int i, j, k;
	for(i=0; i<n; j++)
		for(j=0; j<n; j++)
		{
			res[i][j]=0;
			for(k=0; k<n; k++)
				res[i][j]+=mtr1[i][k]*mtr2[k][j];
		}
}

void mult_matr_vect(int n, double **mtr1, double *v2, double *res)
{
	int i, j;
	for(j=0; j<n; j++)
	{
		res[j]=0;
		for(i=0; i<n; i++)
			res[j]+=mtr1[j][i]*v2[i];
	}
}

int SLAU(double **matrica_a, int n, double *massiv_b, double *x)
{
int i, j, k, r;
double c, M, max, s;
double **a, *b;
a=new double *[n];
for(i=0; i<n; i++)
a[i]=new double[n];
b=new double[n];
for(i=0; i<n; i++)
for(j=0; j<n; j++)
a[i][j]=matrica_a[i][j];
for(i=0; i<n; i++)
b[i]=massiv_b[i];
for(k=0; k<n; k++)
{
	max=fabs(a[k][k]);
	r=k;
	for(i=k+1; i<n; i++)
		if (fabs(a[i][k]>max))
		{
			max=fabs(a[i][k]);
			r=i;
		}
		for(j=0; j<n; j++)
		{
			c=a[k][j];
			a[k][j]=a[r][j];
			a[r][j]=c;
		}
		c=b[k];
		b[k]=b[r];
		b[r]=c;
		for(i=k+1; i<n; i++)
		{
			for(M=a[i][k]/a[k][k], j=k; j<n; j++)
				a[i][j]-=M*a[k][j];
			b[i]-=M*b[k];
		}
}
if(a[n-1][n-1]==0)
if(b[n-1]==0)
return -1;
else return -2;
else
{
	for(i=n-1; i>=0; i--)
	{
		for(s=0, j=i+1; j<n; j++)
			s+=a[i][j]*x[j];
		x[i]=(b[i]-s)/a[i][i];
	}
	return 0;
}
for(i=0; i<n; i++)
delete [] a[i];
delete [] a;
delete [] b;
}

int main(int argc, char* argv[])
{
	int n;
double **mas_a;
double *mas_b;
double *mas_x;
double **matr;
double v[1]={2};
cout << "Ukazhite razmernost matrici:";
cin >> n;
mas_a = new double *[n];
	for( int i=0; i<n; i++)
      mas_a[i] = new double[n];
mas_b = new double [n];
mas_x = new double [n];
matr = new double *[n];
	for(i=0; i< n; i++)
      matr[i] = new double[n];
for (i = 0; i<n; i++)
  for (int j = 0; j < n; j++)
  {
	  cout << "Vvedite element matricy: ";
	  cin >> mas_a[i][j];
  }
for (i=0; i<n; i++)
{
	cout << "Ukazhite razmernost massiva:";
    cin>> n;
	for (i=0; i<n; i++)
	{
	cout << "Vvedite element massiva: ";
    cin >> mas_b[i];
	}
}
tr_matr(mas_a, n);
mult_matr(n, matr, matr, mas_a);
mult_matr_vect(n, mas_a, v, mas_b);
SLAU(matr, n, mas_b, mas_x);
for (i=0; i<n; i++)
	cout << mas_x[i];
return 0;
}

Заранее благодарю

[VIK_aka_TOR]

зависание чаще всего в подобных кодах из-за зацикливания возникает...

[Foralberg]

Вот у меня и не получается разобраться, где зацикливание

[l1pton17]

for (i=0; i<n; i++)
{
cout << "Ukazhite razmernost massiva:";
cin>> n; // скорее всего тут ошибка
for (i=0; i<n; i++)
....
}

изменяешь переменную n, в которой хранится размерность матрицы, и вообще нехорошо в цикле менять переменную, отвечающую за выход из цикла)