def mittelpunkt(A,B,C,D): ''' schnitt der Diagonalen AC und BD ''' M=np.array([A-C,D-B]) R=D-C x=npl.solve(M.T,R) (np.dot(M,x)-R) S=x[0]*A+(1-x[0])*C return S