如果仅考虑xk,应该比较简单，前两句即可。
function [xk,avg] = sequenceOfNumber(k)
xk = (9812*(1:k).^6 +(1:k).^2)./(2*(1:k).^5-(1:k).^4+1); %求序列
xk(xk<793903)=0; %判断，<793903的值改成0
for n=1:k
avg=mean(xk(1:n)) %求xk(1:n)的均值
if avg>793903 %如果均值大于793903，则结束
break;
end
end
xk=xk(1:n); %返回的是均值>793903的最小序列
end
