Вычисление площади треугольника по координатам x,y

#include <iostream>
#include <fstream>
#include <cmath>
#include <stdexcept>
 
struct Point
{
    explicit Point(int const x = 0, int const y = 0) : _x(x), _y(y) {}
    int _x, _y;
};
 
std::istream& operator>> (std::istream& ist, Point& rhs)
{
    return ist >> rhs._x >> rhs._y;
}
 
struct Triangle
{
    explicit Triangle(Point const& a = Point(), Point const& b = Point(), Point const& c = Point()) : _a(a), _b(b), _c(c) {}
 
    double perimeter() const
    {
        return sqrt(pow(_b._x - _a._x, 2) + pow(_b._y - _a._y, 2)) +
                sqrt(pow(_c._x - _b._x, 2) + pow(_c._y - _b._y, 2)) +
                sqrt(pow(_a._x - _c._x, 2) + pow(_a._y - _c._y, 2));
    }
 
    double area() const
    {
        double hp = perimeter() / 2;
        return sqrt(hp *
                    (hp - sqrt(pow(_b._x - _a._x, 2) + pow(_b._y - _a._y, 2))) *
                    (hp - sqrt(pow(_c._x - _b._x, 2) + pow(_c._y - _b._y, 2))) *
                    (hp - sqrt(pow(_a._x - _c._x, 2) + pow(_a._y - _c._y, 2))));
    }
    Point _a, _b, _c;
};
 
std::istream& operator>> (std::istream& ist, Triangle& rhs)
{
    ist >> rhs._a >> rhs._b >> rhs._c;
    if(!ist)
        throw std::runtime_error("operator>> failed");
    return ist;
}
 
int main()
{
    try
    {
        std::fstream fs("in.txt", std::ios::in);
        if(fs.is_open())
        {
            Triangle tr;
            fs >> tr;
            fs.clear();
            fs.close();
            fs.open("out.txt", std::ios::out);
            double area = tr.area();
            std::cout << area << "\n\n";
            fs << area << "\n\n";
            fs.close();
        }
        else
            std::cerr << "Unable to open input file\n";
    }
    catch(std::exception const& exc)
    {
        std::cerr << exc.what() << "\n\n";
    }
    return 0;
}

#include "stdafx.h"
#include <math.h>
#include <iostream>
using namespace std;
#include <string>
int _tmain(int argc, _TCHAR* argv[])
{
 
 
{  double x1,x2,x3,y1,y2,y3, c,p,s,sg,pp,a,b;
  cout<<"x1="; cin>>x1;
  cout<<"y1="; cin>>y1;
  cout<<"x2="; cin>>x2;
  cout<<"y2="; cin>>y2;
  cout<<"x3="; cin>>x3;
  cout<<"y3="; cin>>y3;
 a = sqrt(((x1-x2)*(x1-x2))+((y1-y2)*(y1-y2))); cout<<" a=" <<a<<endl ;
 b = sqrt(((x1-x3)*(x1-x3))+((y1-y3)*(y1-y3))); cout<<" b=" <<b<<endl ;
 c = sqrt(((x2-x3)*(x2-x3))+((y2-y3)*(y2-y3))); cout<<" c=" <<c<<endl ;
 p = a+b+c;
 
pp=p/2;//считаем полупериметр
sg=sqrt(pp*(pp-a)*(pp- b )*(pp-c)); 
s=0.5*(a*b);
 cout<<"ploshad`=" <<s<<endl;
cout<<"ploshad` gerona=" <<sg<<endl;}
system("pause");
    return 0;
}