4、求[100,999]之间能被21整除余1的数，并存储到向量X，对应的下标存入向量Y。（用逻辑关系判断编写，使用find函数，用MATLAB语言书写）

X = [];
Y = [];
for i = 100:999
    if mod(i, 21) == 1
        X = [X i];
        Y = [Y find(X==i)];
    end
end