求线与线的夹角

typedef pair<double,double> CPoint;typedef pair<CPoint,CPoint> CLine;double angle(CLine l,CLine m){    return acos(fabs(project(l.second-l.first,m.second-m.first)/length(l.second-l.first)));}