docs: public void EBEsMat3DL_iRig_jRig(int e) throws JMetalException{
      // The element element 3D than form rigid matrix in the local coordinates
      // i: rigid
      // j: rigid
      // esfuerzos en nudo j por reacciÃƒÂ³n de desplazamientos en i + esfuerzo en i
      // double [][]Kii = new double [5][5];
      // Elements
      // the Element long
      double Lij=Element_[e][L_];
      // Secction
      int idx = (int)Element_[e][INDEX_];
      double S=Groups_[idx][AREA];
      // inertia in axis local z
      double Iz=Groups_[idx][Iz_];
      // inertia in axis local y
      double Iy=Groups_[idx][Iy_];
      // inertia torsion
      double Ip=Groups_[idx][It_];
      // elastic modulus (Young)
      double E=Groups_[idx][E_];
      //  elastic transversal modulus
      double G=Groups_[idx][G_];
      // esfuerzos en nudo j por reacciÃƒÂ³n de desplazamientos en i + esfuerzo en i
      Kii[0][0] = E * S / Lij;
      Kii[0][1] = 0;
      Kii[0][2] = 0;
      Kii[0][3] = 0;
      Kii[0][4] = 0;
      Kii[0][5] = 0;
      Kii[1][0] = 0;
      Kii[1][1] = 12 * E * Iz / Math.pow(Lij, 3);
      Kii[1][2] = 0;
      Kii[1][3] = 0;
      Kii[1][4] = 0;
      Kii[1][5] = -6 * E * Iz / Math.pow(Lij, 2);
      Kii[2][0] = 0;
      Kii[2][1] = 0;
      Kii[2][2] = 12 * E * Iy / Math.pow(Lij, 3);
      Kii[2][3] = 0;
      Kii[2][4] = 6 * E * Iy / Math.pow(Lij, 2);
      Kii[2][5] = 0;
      Kii[3][0] = 0;
      Kii[3][1] = 0;
      Kii[3][2] = 0;
      Kii[3][3] = G * Ip / Lij;
      Kii[3][4] = 0;
      Kii[3][5] = 0;
      Kii[4][0] = 0;
      Kii[4][1] = 0;
      Kii[4][2] = 6 * E * Iy / Math.pow(Lij, 2);
      Kii[4][3] = 0;
      Kii[4][4] = 4 * E * Iy / Lij;
      Kii[4][5] = 0;
      Kii[5][0] = 0;
      Kii[5][1] = -6 * E * Iz / Math.pow(Lij, 2);
      Kii[5][2] = 0;
      Kii[5][3] = 0;
      Kii[5][4] = 0;
      Kii[5][5] = 4 * E * Iz / Lij;
  
      // esfuerzos en nudo i por reacciÃƒÂ³n de desplazamientos en j
      Kij[0][0] = E * S / Lij;
      Kij[0][1] = 0;
      Kij[0][2] = 0;
      Kij[0][3] = 0;
      Kij[0][4] = 0;
      Kij[0][5] = 0;
      Kij[1][0] = 0;
      Kij[1][1] = 12 * E * Iz / Math.pow(Lij, 3);
      Kij[1][2] = 0;
      Kij[1][3] = 0;
      Kij[1][4] = 0;
      Kij[1][5] = -6 * E * Iz / Math.pow(Lij, 2);
      Kij[2][0] = 0;
      Kij[2][1] = 0;
      Kij[2][2] = -12 * E * Iy / Math.pow(Lij, 3);
      Kij[2][3] = 0;
      Kij[2][4] = -6 * E * Iy / Math.pow(Lij, 2);
      Kij[2][5] = 0;
      Kij[3][0] = 0;
      Kij[3][1] = 0;
      Kij[3][2] = 0;
      Kij[3][3] = G * Ip / Lij;
      Kij[3][4] = 0;
      Kij[3][5] = 0;
      Kij[4][0] = 0;
      Kij[4][1] = 0;
      Kij[4][2] = -6 * E * Iy / Math.pow(Lij, 2);
      Kij[4][3] = 0;
      Kij[4][4] = -2 * E * Iy / Lij;
      Kij[4][5] = 0;
      Kij[5][0] = 0;
      Kij[5][1] = -6 * E * Iz / Math.pow(Lij, 2);
      Kij[5][2] = 0;
      Kij[5][3] = 0;
      Kij[5][4] = 0;
      Kij[5][5] = 2 * E * Iz / Lij;
  
      // esfuerzos en nudo j por reacciÃƒÂ³n de desplazamientos en i
      Kji[0][0] = E * S / Lij;
      Kji[0][1] = 0;
      Kji[0][2] = 0;
      Kji[0][3] = 0;
      Kji[0][4] = 0;
      Kji[0][5] = 0;
      Kji[1][0] = 0;
      Kji[1][1] = 12 * E * Iz / Math.pow(Lij, 3);
      Kji[1][2] = 0;
      Kji[1][3] = 0;
      Kji[1][4] = 0;
      Kji[1][5] = -6 * E * Iz / Math.pow(Lij, 2);
      Kji[2][0] = 0;
      Kji[2][1] = 0;
      Kji[2][2] = -12 * E * Iy / Math.pow(Lij, 3);
      Kji[2][3] = 0;
      Kji[2][4] = -6 * E * Iy / Math.pow(Lij, 2);
      Kji[2][5] = 0;
      Kji[3][0] = 0;
      Kji[3][1] = 0;
      Kji[3][2] = 0;
      Kji[3][3] = G * Ip / Lij;
      Kji[3][4] = 0;
      Kji[3][5] = 0;
      Kji[4][0] = 0;
      Kji[4][1] = 0;
      Kji[4][2] = -6 * E * Iy / Math.pow(Lij, 2);
      Kji[4][3] = 0;
      Kji[4][4] = -2 * E * Iy / Lij;
      Kji[4][5] = 0;
      Kji[5][0] = 0;
      Kji[5][1] = -6 * E * Iz / Math.pow(Lij, 2);
      Kji[5][2] = 0;
      Kji[5][3] = 0;
      Kji[5][4] = 0;
      Kji[5][5] = 2 * E * Iz / Lij;
  
      // esfuerzos en nudo i por reacciÃƒÂ³n de desplazamientos en j + esfuerzo en j
      Kjj[0][0] = E * S / Lij;
      Kjj[0][1] = 0;
      Kjj[0][2] = 0;
      Kjj[0][3] = 0;
      Kjj[0][4] = 0;
      Kjj[0][5] = 0;
      Kjj[1][0] = 0;
      Kjj[1][1] = 12 * E * Iz / Math.pow(Lij, 3);
      Kjj[1][2] = 0;
      Kjj[1][3] = 0;
      Kjj[1][4] = 0;
      Kjj[1][5] = -6 * E * Iz / Math.pow(Lij, 2);
      Kjj[2][0] = 0;
      Kjj[2][1] = 0;
      Kjj[2][2] = 12 * E * Iy / Math.pow(Lij, 3);
      Kjj[2][3] = 0;
      Kjj[2][4] = 6 * E * Iy / Math.pow(Lij, 2);
      Kjj[2][5] = 0;
      Kjj[3][0] = 0;
      Kjj[3][1] = 0;
      Kjj[3][2] = 0;
      Kjj[3][3] = G * Ip / Lij;
      Kjj[3][4] = 0;
      Kjj[3][5] = 0;
      Kjj[4][0] = 0;
      Kjj[4][1] = 0;
      Kjj[4][2] = 6 * E * Iy / Math.pow(Lij, 2);
      Kjj[4][3] = 0;
      Kjj[4][4] = 4 * E * Iy / Lij;
      Kjj[4][5] = 0;
      Kjj[5][0] = 0;
      Kjj[5][1] = -6 * E * Iz / Math.pow(Lij, 2);
      Kjj[5][2] = 0;
      Kjj[5][3] = 0;
      Kjj[5][4] = 0;
      Kjj[5][5] = 4 * E * Iz / Lij;
  
      // PrintArchTxtMKLB(e);
  
    }