Функции для вычисления периметра и площади треугольника

#include <stdio.h>
#include <math.h>
float dlina(float q1, float w1, float q2, float w2)
{
             return sqrt((q1-q2)*(q1-q2)+(w1-w2)*(w1-w2));
      }
      
      float pl(float a,float b,float c)
      {
            float p;
            p=(a+b+c)/2;
            return sqrt(p*(p-a)*(p-b)*(p-c));
            }
            
            float pl1(float x1,float y1,float x2,float y2,float x3,float y3)
            {
                  float s;
                  s=pl(dlina(x1,y1,x2,y2),dlina(x1,y1,x3,y3),dlina(x2,y2,x3,y3)); 
                  return s;
                  }     
            
            int main()
            {
                int i,l,j;
                float a[2],b[2],x[2],y[2],s1;
                printf("vv 1 treugolnika: ");
                for (i=0;i<3;i++)
                scanf("%f %f",&a[i],&b[i]);
                printf("\n vv tr: ");
                for(j=0;j<3;j++)
                printf("\n vv 2 treugolnika: ");
                scanf("%f %f",&x[j],&y[j]);
                l=0;
                for (i=0;i<3;i++)
                for(j=0;j<3;j++)
                if((fabs(pl1(x[0],y[0],x[1],y[1],x[2],y[2])-(pl1(a[i],b[i],x[1],y[1],x[2],y[2]))-
                (pl1(a[i],b[i],x[0],y[0],x[1],y[1]))-
                (pl1(a[i],b[i],x[0],y[0],x[2],y[2])))<=0.0001)||(fabs(pl1(a[0],b[0],a[1],b[1],a[2],b[2])-pl1(x[j],y[j],a[1],b[1],a[2],b[2])-
                pl1(x[j],y[j],a[0],b[0],a[1],b[1])-pl1(x[j],y[j],a[0],b[0],a[2],b[2]))<=0.0001))
                {printf("yes");
                s1=fabs(pl1(a[0],b[0],a[1],b[1],a[2],b[2])-pl1(x[0],y[0],x[1],y[1],x[2],y[2]));
                printf("\n %f", s1);}
                else printf("no");
                getchar();
                getchar();
                }