Remarque sur Ikl

Dans l'expression précédente, si i0,i1,i2 sont les numéros des 3 sommets du triangle Tk

 Ikl = (t[l].area/3) * (f[i0] + f(i1))/2 *(1/2) *(!(i2==k))
                  +(f[i1] + f(i2))/2 *(1/2) *(!(i0==k))
                  +(f[i2] + f(i0))/2 *(1/2) *(!(i1==k));

car wk vaut 0 sur l'arête opposée a k et 1/2 sur les 2 autres arêtes. Cette expression se simpifie en

 Ikl = t[l].area * (f[i0] + f(i1) + f[i2])/12
                  + t[l].area * f[k] / 12;