求两条直线交点

procedure P_zhixianjiaodian(const x1, y1, x2, y2, x3, y3, x4, y4: double; var Solx1, Soly1: double);
var
  a1, b1, c1, a2, b2, c2: double;
begin
  a1 := y2 - y1;
  b1 := x1 - x2;
  c1 := x2 * y1 - x1 * y2;
  a2 := y4 - y3;
  b2 := x3 - x4;
  c2 := x4 * y3 - x3 * y4;
  Solx1 := (b2 * c1 - b1 * c2) / (a1 * b2 - a2 * b1);
  Soly1 := (a1 * c2 - a2 * c1) / (a1 * b2 - a2 * b1);

end;