大佬们，关于微分方程拟合问题求教

xlist = Table[i, {i, 0, 20}];
ylist = Table[3*i^3 + 2*i^2, {i, 0, 20}]
data = Table[{xlist[[i]], ylist[[i]]}, {i, 1, 21}]；
eq = y'[x] == a*x^m + b*x^n
pfun = ParametricNDSolveValue[{eq, y[0] == 0}, y, {x, 0, 20}, {a, b, m, n}]
sol = NonlinearModelFit[data, pfun[a, b, m, n][x], {a, b, m, n}, x]
为什么这么简单的a=9,b=4,m=2,n=1，程序就是拟合不出来呢？大佬们看看我的代码有没有问题？