示例#1
0
void EliminaGaussiana(Matriz& m) {
    for (int i = 0; i < m.size(); ++i) {
        int fila_mayor = i;
        for (int j = i + 1; j < m.size(); ++j)
            if (fabs(m[fila_mayor][i]) <
                fabs(m[j][i])) fila_mayor = j;
        swap(m[i], m[fila_mayor]);
        if (Igual(m[i][i], 0)) continue;
        for (int j = m[i].size() - 1; j >= i; --j)
            m[i][j] /= m[i][i];
        for (int j = 0; j < m.size(); ++j) {
            if (i == j || Igual(m[j][i], 0)) continue;
            double ratio = m[j][i] / m[i][i];
            for (int k = i; k < m[j].size(); ++k)
                m[j][k] -= m[i][k] * ratio;
        }
    }
}
void EliminacionGaussiana(Matriz& m) {
  for (int i = 0; i < m.size(); ++i) {
    // <comment>
    int fila_mayor = i;
    for (int j = i + 1; j < m.size(); ++j)
      if (fabs(m[fila_mayor][i]) <
        fabs(m[j][i])) fila_mayor = j;
    swap(m[i], m[fila_mayor]);
    // </comment>

    if (EsCero(m[i][i])) continue;
    for (int j = m[i].size() - 1; j >= i; --j)
      m[i][j] = m[i][j] / m[i][i];
    for (int j = 0; j < m.size(); ++j) {
      if (i == j || EsCero(m[j][i])) continue;
      for (int k = m[j].size() - 1; k >= i; --k)
        m[j][k] = m[j][k] - m[i][k] * m[j][i];
    }
  }
}
示例#3
0
void printVector(Matriz &matriz){
int fila = 0,columna = 0;
     for(int i=0; i<matriz.size(); ++i){
          for(int j=0; j<matriz[0].size(); ++j){
//               cout << "[" << matriz[i][j] << "]";
               cout << "[" << matriz.at(i).at(j) << "]";
               fila += matriz.at(i).at(j);
          }
          columna = sumaColumna(matriz, i);
          cout << " suma por fila= " << fila << ", suma por columna= " << columna;
          columna = 0;
          fila = 0;
          cout <<"\n";
     }
}