请找出1到1000中满足a^2-100b=a的正整数a,b

请找出1到1000中满足a^2-100b=a的正整数a,b。编写函数依次返回行向量a，b（其中a(i），b(i)为第i组解），其中数组a中的元素是递增排列。
我写的是
a(1)=1;b(1)=1;
for k=1:1000
a(k)=k;
for m=1:1000
b(m)=m;
if a(k)^2-100*b(m)==a(k)
x=[k,m];
disp(x)
end
end
end
输出的结果是单个的行向量，怎么将所有行向量整合成一个矩阵？请给个代码参考一下，刚学matlab两天，不是很了解

x_all = [];
a(1)=1;b(1)=1;
for k=1:1000
    a(k)=k;
    for m=1:1000
        b(m)=m;
        if a(k)^2-100*b(m)==a(k)
            x = [k,m];
            x_all = [x_all;x];
        end
    end
end
disp(x_all);