請各位大大 幫我看為啥這一段程式會出錯 謝謝

clear

N=[ 1 2 3 4 5
1 2 3 4 5
2 5 4 6 9
9 4 6 9 7 ]

[m,n]=size(N)
for i=1:m
for j=1:n
d1(i,j)=abs(((N(i-1,j) N(i,j-1))/2)-N(i,j) );
if d1(i,j)>1
n1(i,j)=ceil(log2(d1(i,j)));
d2(i,j)=2^n1(i,j) 1;
N(i,j)=(N(i-1,j) N(i,j-1))/2-d2(i,j);
else
N(i,j)=N(i,j)
end
end
end

會出現 ??? Subscript indices must either be real positive integers or logical
請問哪裡有出錯呢!!??

像d1(i,j) 和n1(i,j) d2(i,j) 這種有需要加(i,j)嗎
還是只要打d1就好呢?

還有
else N(i,j)=N(i,j) 這段 需要嗎 (這邊是說 d1(i,j)沒有>1 原本地方的數值不變)

謝謝大家^^

矩陣元素 N(0,0) 不存在.

只存在 N(1,1) ~ N(4,5)