Dans P1::intgrad2

Cette fonction calcule l'integrale sur le domaine du carre du gradient de la fonction. Il faut maintenant que cela soit multiplié par nu et qu'il lui soit ajouté l'integrale sur le domaine de a fois la fonction au carre. Les constantes nu et a doivent donc etre des parametres de la fonction.

P1::intgrad2(const float& nu, const float& a)

Le principe de calcul de la 2eme integrale I est le suivant

I=0;
for(int k=0; k<g->nt; k++)
{  
  int i0 = g->no(g->t[k].v[0]);
  int i1 = g->no(g->t[k].v[1]);
  int i2 = g->no(g->t[k].v[2]);
  I+= g->t[k].area * a *( (cc[i0]+cc[i1])*(cc[i0]+cc[i1])
         + (cc[i1]+cc[i2])*(cc[i1]+cc[i2]) + (cc[i2]+cc[i0])*(cc[i2]+cc[i0]) 
         )/12;
 }