matlab循环给数组赋值

hang=size(I,1); %I的行数
lie=size(I,2); %Id列数
for (ii=1:1:hang)
for (jj=1:1:lie)
if ((ii-ox)*(ii-ox)+(jj-oy)*(jj-oy) P(ii,jj)=1
else
P(ii,jj)=0
end
end
end

这么个循环怎么总是不完
P数组出来一次接着一直出 这是为什么？

我是想把I这个图中的一个圆内赋值1，其余部分赋值0，ox，oy是圆心坐标R是半径

N=4000，q=0.98
for n=1:N;
for i=1:n;
j=i-1;

if
j<=0
a(i,n+1)=n^(q+1)-(n-q)*(n+1)^q;
else
a(i,n+1)=(n-j+2)^(q+1)+(n-j)^(q+1)-2*(n-j+1)^(q+1);
end

end
end
希望对你能有所帮助。