una solucion

float area_total=0;
 for (int k=0;k<num_tri;k++){
     float sp=0;
     int vert1,vert2,vert3;
     float a,b,c,area;
     float dx,dy;
     vert1=T[k][0];
     vert2=T[k][1];
     vert3=T[k][2];
     dx= P[vert2][0]-P[vert1][0];
     dy= P[vert2][1]-P[vert1][1];
     a=sqrt(dx*dx+dy*dy);
     dx= P[vert3][0]-P[vert2][0];
     dy= P[vert3][1]-P[vert2][1];
     b=sqrt(dx*dx+dy*dy);     
     dx= P[vert1][0]-P[vert3][0];
     dy= P[vert1][1]-P[vert3][1];
     c=sqrt(dx*dx+dy*dy);
     sp=0.5*(a+b+c)     ;
     area = sqrt(sp*(sp-a)*(sp-b)*(sp-c));
     area_total+=area;
 }
 cout<<"Area total"<<area_total<<endl;