jestem początkująca w c++, chodzę do 1 klasy lo, ale bardzo zaintrygowała mnie ta metoda, poniżej wklejam kod
(Rozwiązywanie układów równań metodą eliminacji Gaussa)
int main()
{
double a1,a2,b1,b2,c1,c2,wx,wy,wz,w,x,y,t,z,b,j,i,n;
cout<<"Podaj a1: ";
cin>>a1;
cout<<"Podaj b1: ";
cin>>b1;
cout<<"Podaj c1: ";
cin>>c1;
cout<<"Podaj a2: ";
cin>>a2;
cout<<"Podaj b2: ";
cin>>b2;
cout<<"Podaj c2: ";
cin>>c2;
{
void gauss(double **t, double *b)
for (int i=1; i<=n-1; i++);
for (int j=i+1; j<=n; j++);
for (int k=i+1; k<=n; k++)
t[j-1][k-1] -= t[i-1][k-1]*t[j-1][i-1]/t[i-1][i-1];
b[j-1] -= b[i-1]*t[j-1][i-1]/t[i-1][i-1];
}
b[n-1] /= t[n-1][n-1];
s=0;
for (int i=1; i<=n-1; i++){
for (int j=n+1-i; j<=n; j++)
s += t[n-i-1][j-1]*b[j-1];
b[n-i-1] = (b[n-i-1] - s)/t[n-i-1][n-i-1];
if(w!=0)
{
cout<<"x = "<<wx/w<<endl;
cout<<"y = "<<wy/w<<endl;
cout<<"z = "<<wz/w<<endl;
}
else
if(wx==0&&wy==0&&wz==0)
cout<<"Uklad ma nieskonczenie wiele rozwiazan"<<endl;
else
cout<<"Uklad sprzeczny"<<endl;
system("PAUSE");
return 0;
}
Chciałabym poprosić o wskazanie błędów i o wskazówki tego co dodać do kodu aby dobrze wyliczał. Proszę o jak najszybsze odpowiedzi.