{

return 0.5 * ((get<0>(t[1]) - get<0>(t[0])) * (get<1>(t[2]) - get<1>(t[0]))

- (get<0>(t[2]) - get<0>(t[0])) * (get<1>(t[1]) - get<1>(t[0])));

{

//triangle<boost::geometry::point_xy<double> > t;

triangle<boost::tuple<double, double> > t;

t[0] = boost::make_tuple(0, 0);

t[1] = boost::make_tuple(5, 0);

t[2] = boost::make_tuple(2.5, 2.5);

std::cout << "Triangle: " << boost::geometry::dsv(t) << std::endl;

std::cout << "Area: " << boost::geometry::area(t) << std::endl;

//boost::geometry::point_xy<double> c;

boost::tuple<double, double> c;

boost::geometry::centroid(t, c);

std::cout << "Centroid: " << boost::geometry::dsv(c) << std::endl;

return 0;