使用numpy.linalg.lstsq进行4阶2D数据拟合时输出系数顺序的疑问

我在用NumPy的numpy.linalg.lstsq处理4阶2D数据的最小二乘拟合时，对返回的系数顺序有点困惑。

我已经拿到了拟合后返回的系数（具体数值在下方），能看出来第一个数值≈2936397对应的是x⁰y⁰项，但剩下的系数对应的(x,y)阶数到底是按什么规则排列的呢？

我原本猜测顺序是这样的：

Coefficient,x-order,y-order
(0,0), (1,0), (2,0), (3,0), (4,0), (0,1), (1,1), (2,1), (3,1), ... (3,4), (4,4)

但实际测试后发现这个猜想不对，有没有懂的朋友帮忙解答一下？提前感谢！

2936397.6942270664
-15121.957463652841
9.03102972423578
-0.012272870202833481
0.00010249506244899587
-114608.90940722093
468.86962365232756
-0.17328878542700762
-0.00025389866132422223
-4.259024111721344e-08
2460.5890909300783
-7.124697975197747
0.001634958550849398
4.363571901584845e-07
3.193596492720019e-08
-23.96706910671856
0.03889662608251326
8.892484402167611e-07
-3.9098540582926944e-07
-7.002745405770838e-08
0.08027271303336647
-3.2947609724942595e-06
-1.3261478670756333e-07
-1.8548462321632542e-08
2.465621946612373e-07

备注：内容来源于stack exchange，提问作者abc def