double area(Point a,Point b,double y)
{
	
	Point u=Point(a.x,b.y);
	double ret=fabs(cross(a,b,a,u))*y*y/((a.y-b.y)*(a.y-b.y));
	return ret/2;
}