for i in range(max_degree): poly_features[:, i] /= math.gamma(i + 1)

这段代码是将多项式特征矩阵中每列的每个元素除以该列的阶乘。具体来说，对于矩阵中的每个元素 $x$，它所在的列的阶数为 $i$，那么 $x$ 除以 $i$ 的阶乘（$i!$）就是该元素的新值。这样可以将原始特征矩阵转换为多项式特征矩阵，以便更好地拟合非线性模型。其中，math.gamma(i + 1) 是 Python 标准库中的 gamma 函数，计算 $i+1$ 的阶乘。