Код по алгоритму

[@Manya@]

Вот, собственно алгоритм,

вот код, вроде должно работать
Обрабатывается массив p_Q[deg][deg]

int * c = new int [deg];
	int i, j, r, k, flag, p;
	flag=0;
	r=0;
	for (i=0; i<deg; i++)
		c[i]=-1;
	for (k=0; k<deg; k++){
			for (j=0; j<deg; j++){
				if (this->p_Q[k][j]!=0 && c[j]<0){
					flag++;
					if (flag)
						break;
				}
				if (flag){
					for (i=0; i< deg; i++)
						this->p_Q[i][j]*=(-1/(this->p_Q[k][j]));
					for (i=0; i<deg; i++)
                                                   for(p=0; p<deg; p++)
							if (i!=j)
								this->p_Q[p][i]+=this->p_Q[p][j]*this->p_Q[k][i];
					c[j]=k;
				}
			}
			flag=0;
		}

[@Manya@]

есть вариант еще , но это не лучше

int * c = new int [deg];
	int i, j, r, k, flag, p;
	flag=0;
	r=0;
	for (i=0; i<deg; i++)
		c[i]=-1;
	for (k=0; k<deg; k++){
			for (j=0; j<deg; j++){
				if (this->p_Q[k][j]!=0 && c[j]<0){
					flag++;
					if (flag)
						break;
				}
				if (flag){
					for (i=0; i< deg; i++)
						this->p_Q[i][j]*=(-1/(this->p_Q[k][j]));
				for (i=0; i<deg; i++)
					if (i!=j)
						this->p_Q[i][j]*=this->p_Q[k][i];
				for (p=0; p<deg; p++)
					for (i=0; i<deg; i++)
						if (i!=j)
							this->p_Q[p][i]+=this->p_Q[p][j];
					c[j]=k;
				}
			}
			flag=0;
		}