Poryjemy..zrobaczymy

1. ok, ale w książce , Soulburner 21/08/08 21:32
   fragment dotyczący xn jest podany "luzem", nie w pętli. Zamiana //x[n]= na
   x[n]= tabA[n][n+1] / tabA[n][n];
   nie pomaga w obliczeniach - iksy wciąż wychodzą zerowe :/

   You must gather your party before
   venturing forth.

   1. No pewnie ze luzem , JOTEM 22/08/08 00:05
      i tak ma być - to rekurencyjne liczenie poza pierwszym elementem - ten jest wyliczany jawnie i poza pętlą.
      Można włączyć do pętli, ale trzeba indeksy odpowiednio rozszerzyć.
      
      Poza tym, Twoja książka ma błąd i Twój kod ma inne jeszcze błędy. Analizuję ten starszy kod:
      - brak wyliczeń pierwszego wyrazu (xn)
      - błąd w książce w pętli
      jest for (p=i+1; p<=1; p++)
      powinno być for (p=i+1; p<=n-1; p++)
      - zastanów się, czy pętle mają być od i=1 do n czy od i=0 do n-1
      bo jesteś niekonsekwentny (w książce masz np od 1 do n-1, a Ty używasz od 0 do n-1 - jeden obieg za dużo).
      
      To tyle - popraw, to zadziała.

      Poryjemy..zrobaczymy

2. ok, coś tam wymodziłem i nawet działa, tylko , Soulburner 21/08/08 23:06
   wyniki zaokrągla do liczb całkowitych, podczas gdy liczenie na piechotę daje np.
   x3 = 2,8
   zaś program wyświetla x3 = 3
   
   Tutaj zmodyfikowany kod programu: http://daath.x.pl/kodzik2.cpp
   
   Ciekawostka - mityczny //x[n]= okazał się zbędny, bo... czy linijka:
   x[n-1] = tabA[n-1][n] / tabA[n-1][n-1]; //jak widać, trochę inna niż w poście wyżej
   jest, czy jej nie ma, wyniki wszystkie są :P

   You must gather your party before
   venturing forth.

   1. Tak mniej więcej , JOTEM 22/08/08 00:30
      krytyczny fragment po poprawnym przeindeksowaniu z książki z indeksowanie pascalowym (1..n) do indeksowania c++(0..n-1):
      
      x[n-1] = (tabA[n-1][n] - w) / tabA[n-1][n-1];
      for (i=n-2; i>=0; i--)
      {
      double w=0;
      for (p=i+1; p<n; p++)
      {
      w = w + tabA[i][p] * x[p];
      }
      x[i] = (tabA[i][n] - w) / tabA[i][i];
      }

      Poryjemy..zrobaczymy

      1. nie usunałem w z x(n-1) , JOTEM 22/08/08 00:32
         x[n-1] = tabA[n-1][n] / tabA[n-1][n-1];

         Poryjemy..zrobaczymy