Resolución de Sistemas de Ecuaciones Lineales y Cálculo Numérico

Problema 41

if(A.dim1()==0 || A.dim1()!=A.dim2()) return 0.;

if(A.dim1()==1) return A[0][0];

real determinante=0.;
  for(int k=0;k    Array2D B(A.dim1()-1,A.dim1()-1);
    for(int i=0;i      for(int j=0;j        if(j        else  B[i][j]=A[i+1][j+1];
      }
    }
    if(k%2==0) determinante+=A[0][k]*mn_determinante_recursivo(B);
    else determinante-=A[0][k]*mn_determinante_recursivo(B);
  }
  return determinante; Seguir leyendo “Resolución de Sistemas de Ecuaciones Lineales y Cálculo Numérico” »