[C] Algoritmo risoluzione sistemi lineari

[Neptune]

Originariamente inviato da linoma
Per il tuo metodo di risoluzione qui, nn potresti dire dov'e' l'errore?

Inanzi tutto nel penultimo ciclo, questo:

for(index=0;index<n;index++)
              {
                 var_appo=leggi_elem(M,i,index);
                 var_appo=var_appo/leggi_elem(M,i,j);
                 scrivi_elem(M,i,index,var_appo);
              }

Mi sono accorto che c'era un errore, almeno da quel che dice il flowchart, perchè devo moltiplicare per il reciproco mentre li fa una semplice divisione, e quindi l'ho sostituito con questo:

for(index=0;index<n;index++)
              {
                 var_appo=leggi_elem(M,i,index);
                 var_appo=var_appo * (1/leggi_elem(M,i,j));
                 scrivi_elem(M,i,index,var_appo);
              }

Il problema però persiste e non riesco a trovare altri errori