Años que no hacía una integral doble :-)
Primero, ubico mi sistema de coordenadas en un vértice en lugar de el
baricentro y expreso cada punto de la forma:
uA+vB ; con u,v reales y A,B vectores correspondientes a los otros
vértices.
La distancia^2 al vértice es entonces:
r^2 = u^2 a^2 + v^2 b^2 + 2 uv A·B
Así que el momento de inercia respecto a ese vértice es:
Jver = AxB int(int(r^2,v,0,1-u),v,0,1) ; AxB corresponde al
jacobiano. 'x': prod. vectorial.
Jver = AxB (a^2 + b^2 + A·B)/12
Y por el teorema se Steiner:
Jbar = Jver - Rbar^2 S ; Rbar: distancia al baricentro. S:
superficie del triángulo (que se obtiene con la fórmula de Herón)
Como AxB=2S , A·B=(a^2+b^2-c^2)/2 y Rbar=|A+B|/3
la expresión final resulta:
Jbar = S (a^2+b^2+c^2)/36