title: LU decomposition in c++
body: <p>I am currently having some probs with my code to decompose an array into an upper(u) and lower(l) array.</p>
  
  <p>I am using the doolittle method</p>
  
  <p>My code:</p>
  
  <pre><code>#include &lt;iostream&gt;
  
  using namespace std;
  
  int main(){
  
      double a[10][10];
      double l[10][10];
      double u[10][10];
      double som=0;
      int RANG;
      cin &gt;&gt; RANG;
      for(int i=0; i&lt;RANG; i++){
              for(int j=0; j&lt;RANG;j++){
                      cin &gt;&gt; a[i][j];
              }
      }
      for(int i=0; i&lt;RANG; i++){
              for(int j=0; j&lt;RANG;j++){
                      u[i][j]=0;
                      l[i][j]=0;
              }
      }
  
      for(int i=0;i&lt;RANG;i++) {
      l[i][i]=1;
  
          for(int j=i;j&lt;RANG;j++) {
              for(int s=0;s&lt;i-1;s++) {
                  som+= l[i][s]*u[s][j];
              }
              u[i][j]=a[i][j]-som;
          }
  
          for(int k=i+1;k&lt;RANG;k++) {
              double som=0;
              for(int s=0;s&lt;i-1;s++) {
                  som+=l[k][s]*u[s][i];
              }
              l[k][i]=(a[k][i]-som)/u[i][i];
          }
       }
       cout &lt;&lt; "l:" &lt;&lt; endl;
      for(int i=0; i&lt;RANG; i++){
              for(int j=0; j&lt;RANG;j++){
                      cout &lt;&lt; l[i][j] &lt;&lt; "\t";
              }
              cout &lt;&lt; endl &lt;&lt; endl;
      }
      cout &lt;&lt; "u: " &lt;&lt; endl;
          for(int i=0; i&lt;RANG; i++){
              for(int j=0; j&lt;RANG;j++){
                      cout &lt;&lt; u[i][j] &lt;&lt; "\t";
              }
              cout &lt;&lt; endl &lt;&lt; endl;
      }
      return 0;
  }
  </code></pre>
  
  <p>plz help if you can...</p>
  
  <p>PS: not sure if it belongs here, might be better on the mathematics site</p>