В данном множестве точек узнать, какие из них создадут треугольник наименьшей площади

#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
 
 
int _tmain(int argc, _TCHAR* argv[])
{
    float x[255];
    float y[255];
    float a;
    float b;
    float c;
    int n = 0;
    float asd = 18;
    float s;
    float st;
    float p;
    int i; 
    int j;
    int k;
        printf("vvedite koordinatu tojek\n");
        scanf_s("%d, &n");
        if (n>=3&&n<=100) {
        for(i=1;i<=n;i++) {
        scanf_s("%f, %f\n", &x[i], &y[i]);
        }
 
        s = exp(asd);
    
        for(j=1;j<=n;j++)
            {
                for(k=1;k<=n;k++)
                {
                    if(i!=j||j!=k||j!=k) {
                    
                    a=sqrt((x[i]-x[j])*(x[i]+x[j])-(y[i]-y[j])*(y[i]+y[j]));
                    b=sqrt((x[i]-x[k])*(x[i]+x[k])-(y[i]-y[k])*(y[i]+y[k]));
                    c=sqrt((x[j]-x[k])*(x[j]+x[k])-(y[j]-y[k])*(y[j]+y[k]));
                    p=(a+b+c)/2;
                    st=sqrt(p*(p-a)*(p-b)*(p-c));
                    if(st<s) {
                        s=st;
                    }
                    }
                }
 
        }
    }
        
        printf("%.2f\n", s);
    return 0;
}

scanf_s("%d", &n);

for(i = 0; i < n; i++) 
{
    scanf_s("%f%f", &x[i], &y[i]);
}

#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
 
int _tmain(int argc, _TCHAR* argv[])
{
    float x[100];
    float y[100];
    float a;
    float b;
    float c;
    int n = 0;
    float asd = 18;
    float s;
    float st;
    float p;
    int i; 
    int j;
    int k;
        printf("vvedite koordinatu tojek\n");
        scanf_s("%d", &n);
        if (n>=3&&n<=100) {
        for(i=1;i<n;i++) {
        scanf_s("%f, %f\n", &x[i], &y[i]);
        }
 
        s = exp(asd);
    
        for(j=0;j<=n;j++)
            {
                for(k=0;k<=n;k++)
                {
                    if(i!=j||j!=k||j!=k) {
                    
                    a=sqrt((x[i]-x[j])*(x[i]+x[j])-(y[i]-y[j])*(y[i]+y[j]));
                    b=sqrt((x[i]-x[k])*(x[i]+x[k])-(y[i]-y[k])*(y[i]+y[k]));
                    c=sqrt((x[j]-x[k])*(x[j]+x[k])-(y[j]-y[k])*(y[j]+y[k]));
                    p=(a+b+c)/2;
                    st=sqrt(p*(p-a)*(p-b)*(p-c));
                    if(st<s) {
                        s=st;
                    }
                    }
                }
 
        }
    }
        
        printf("%.2f\n", s);
    
return 0;
system("PAUSE"); 
getch();
}

for(i=1;i<n;i++) {
 scanf_s("%f, %f\n", &x[i], &y[i]);
 }