mathematica最小二乘法曲线拟合问题

mathematica最小二乘法曲线拟合问题
x = Table[100 + 10*i, {i, 0, 9}];
y = {45, 51, 54, 61, 66, 70, 74, 78, 85, 89};
xy = Table[{x[[i]], y[[i]]}, {i, 1, 10}]
ListPlot[xy, PlotStyle -> PointSize[0.015]]
q[a_, b_] := Sum[(b + ax[[i]] - y[[i]])^2, {i, 1, 10}]
Solve[{D[q[a, b], a] == 0, D[q[a, b], b] == 0}, {a, b}]
按理说应该可以解出来a,b,为什么出现了这个：
“Part::partd: 部分指定 ax[[1]] 比对象深度更长. >>
Part::partd: 部分指定 ax[[2]] 比对象深度更长. >>
Part::partd: 部分指定 ax[[3]] 比对象深度更长. >>
General::stop: 在本次计算中,Part::partd 的进一步输出将被抑制. >>
Solve::svars: 方程可能无法给出所有 "solve" 变量的解. >>
哪里出了问题