contrôle

# contrôle global
print M * fp
matrix([[ 1616411.9125,1617482.5175,1617503.8225,1616433.2175],
        [ 9173672.32999999,9173641.84,9174389.85,9174420.34]])
 
# contrôle par point
xfp = np.matrix([[1],[8810],[1]])
xtp  = M * x0
print xtp
matrix([[ 1616433.2175],
        [ 9174420.34  ]])
 
# nouveau point
xfp = np.matrix([[0],[8810],[1]])
xtp = M * xfp
print xtp
matrix([[ 1616433.13258526],
        [ 9174420.34241831]])