матрица Си++

[xVeteRx]

Слушайте, такая задачка: есть треугольная матрица, даны её коеффициенты и правая часть(результаты) найти корни
#include <stdlib.h>
#include <conio.h>
#include <iostream.h>
#include <math.h>


void main()
{
clrscr();
randomize();

int a[19][20],b[20],s=0;

for (int i=0;i<19;i++)
{
for(int j=0;j<20;j++)
{
a[i][j]=rand()%3;
if (i>j) a[i][j]=0;
if (i==j && a[i][j]==0) a[i][j]+=1;
cout<<a[i][j]<<" ";
}
cout<<endl;
}

for(int n=19;n>=0;n--)
{
b[n]=(a[n][20]+s)/a[n][n];
s=s+a[n][n]*b[n];
}
cout<<"korni:";
for (int i = 0; i < 20; i++)
{
cout<<" "<<b[i];
}
getch();

}
подкинте идейку как решить,плиз

[Sazary]

Как я понимаю, речь идет о системе линейных уравнений?
Почитайте про метод Гаусса-Жордана.
А вот тут -> http://ru.wikipedia.org/wiki/Метод_Гаусса_—_Жордана очень понятно описан алгоритм.
Сам когда-то делал эту задачу.

[xVeteRx]

#include <stdlib.h>
#include <conio.h>
#include <iostream.h>
#include <math.h>


void main()
{
clrscr();
randomize();

int a[19][20],b[20],s=0;

for (int i=0;i<19;i++)
{
for(int j=0;j<20;j++)
{
a[i][j]=rand()%3;
if (i>j) a[i][j]=0;
if (i==j && a[i][j]==0) while(a[i][j]==0)a[i][j]=rand()%3;
cout<<a[i][j]<<" ";
}
cout<<endl;
}

for (int i=0;i<20;i++)
{
a[0][i]+=a[20][i]*(-a[20][20]);
}
cout<<endl<<endl;
for (int i=0;i<19;i++)
{
for(int j=0;j<20;j++)
{
cout<<a[i][j]<<" ";
}
cout<<endl;
}


getch();
}
пока так но чёт не то(

[Sazary]

Вам нужно преобразовать матрицу к единичной.
Делайте как написано в алгоритме.

[xVeteRx]

Всё вродь добил РАБОТАЕТ))))
я кароч с размером матрицы чуть ошибся и у меня числа левые из памяти вызывались вот и запутался.
#include <stdlib.h>
#include <conio.h>
#include <iostream.h>
#include <math.h>


void main()
{
clrscr();
randomize();

float a[20][21];

for (int i=0;i<20;i++)
{
for(int j=0;j<21;j++)
{
a[i][j]=rand()%100;
if (i>j) a[i][j]=0;
if (i==j && a[i][j]==0) a[i][j]+=1;
cout<<a[i][j]<<" ";
}
cout<<endl;
}


for (int i=0;i<20;i++)
{
int b=a[i][i];
for(int j=0;j<21;j++)
{
a[i][j]=a[i][j]/b;
}
}
for (int z=0; z < 19; z++)
{

for (int k = 18; k >= 0; k--)
{

float b=a[k-z][19-z];
for(int n=0;n<21;n++)
{
a[k-z][n]=a[k-z][n]-a[19-z][n]*b;


}
}
}
cout<<"korni ravni: "<<endl;
for (int i = 0; i < 20; i++)
{
cout<<"x"<<i+1<<"="<<a[i][20]<<endl;
}
getch();

}
сори за оформление только сделал рад безмерно))

[xVeteRx]

Спасиб что помог,кстат)) алгоритм пригодился)

[Sazary]

xVeteRx, Рад за вас )
Только оформляйте код в теги CODE.
Так будет намного читабельнее. )