def METHOD_NAME(self, points):
    """Fit plane to points
    Parameters
    ----------
    points : list of ndarray
        Each array is a vector containing the x,y,z coordinates of 1 point
    Returns
    -------
    np.ndarray (float64)
        (vector of length 4) Coefficients of plane
    """
    axyz = np.ones((len(points), 4))
    axyz[:, :3] = points
    m = np.linalg.svd(axyz)[-1][-1, :]
    if m[0] < 0:
        m = m * (-1)
    self.from_coeffs(m)
    return self.coeffs