质心法求光斑中心中遇到的问题

a=imread('e:\01.bmp');
    b=medfilt2(a,[9,9]);
[m,n]=size(b);
for i=1:m;
    for j=1:n;
   e(i,j)=b(i,j)-67;   
        end
end
su=0;
di=0;
    for i=1:m;
    for j=1:n;
f==j*e(i,j)
z=e(i,j);
su=su+f;
di=di+z;
        end
    end
x=su/di
我编的求光斑中心的部分程序，不知为什么su和di的值都是255