for(int i=0; i<=1; ++i){
     for(int l=0; l<=1 ; ++l){
      for(int p=0; p<NB_CORRESP; ++p){

       point1(0) = list1(p,0); point1(1) = list1(p,1); point1(2) = list1(p,2);
       point2(0) = list2(p,0); point2(1) = list2(p,1); point2(2) = list2(p,2);
       point3(0) = list3(p,0); point3(1) = list3(p,1); point3(2) = list3(p,2);

       for(int k=0; k<=2; ++k){
        A(2,l*p) = point1(k)*point2(i)*point3(2);
        A(i,l*p) = - (point2(2)*point3(2));
        A(2,2*p) = -(point2(i)*point3(l));
        A(i,2*p) = point2(2)*point3(l);

        // EQ += point1(k)*(point2(i)*point3(2)*T(2,l,k) - point2(2)*point3(2)*T(i,l,k) - point2(i)*point3(l)*T(2,2,k) + point2(2)*point3(l)*T(i,2,k));
       }
      }
     }
    }