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” »
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” »