向量为什么会有未定义的问题？？

[gkyingkyin]

我编的程序，最后用到向量r,按说向量应该不用定义吧，可是不知为什么，一运行就显示

??? Undefined function or variable 'r'.

Error in ==> fengle at 248
r

请高人解答，谢谢！！急！！
以下是我的程序：

I=imread('ggg.bmp');
J1=I;
J2=I;
J3=I;
M1=rgb2gray(J1);
N1=histeq(M1);
for i=200:400
if N1(i,1)<50
Q=i;
end;
end;
I=double(I);
for i=1:480
for j=1:640
IR=I(i,j,1)/255;
IG=I(i,j,2)/255;
IB=I(i,j,3)/255;
k=0.21267/[0.21267*IR+0.71516*IG+0.07217*IB];
R=IR*k;
G=IG*k;
B=IB*k;
x=0.41245*R+0.35758*G+0.18042*B;
y=0.21267*R+0.71516*G+0.07217*B;
z=0.01933*R+0.11919*G+0.95023*B;
a=500*[myfun(x/0.95045)-myfun(y)];
b=200*[myfun(y)-myfun(z/1.08875)];
L=myfun(y)*116-16;
if (I(i,j,1)-I(i,j,3)>25)&(I(i,j,1)-I(i,j,2)>25)&(b/a>0)&(b/a<1/2)&(a^2+b^2>25)&i<Q
J1(i,j,1)=255;
J1(i,j,2)=0;
J1(i,j,3)=0;
end;
end;
end;
for i=1:480
for j=1:640
if (I(i,j,1)-I(i,j,3)>25)&(I(i,j,1)-I(i,j,2)>25)
J2(i,j,1)=255;
J2(i,j,2)=255;
J2(i,j,3)=255;
end;
end;
end;
M2=rgb2gray(J2);
P=rgb2gray(J3);
N2=histeq(M2);
for i=200:400
if N2(i,1)<50
Q=i;
end;
end;
c=0;
A=0;
B=0;
for i=Q-50:Q-10
for j=2:630
if N2(i,j)<90&(N2(i-2,j)-N2(i,j)>35)&(N2(i+2,j)-N2(i,j)>35)
c=c+1;
A=A+j;
B=B+i;
end;
end;
end;
C=A/c;
D=B/c;
E=0;
F=0;
for i=Q-50:Q-10
for j=2:630
if N2(i,j)<90&(N2(i-2,j)-N2(i,j)>35)&(N2(i+2,j)-N2(i,j)>35)
E=E+(j-C)*(i-D);
F=F+(j-C)^2;
end;
end;
end;
G=E/F;
H=D-G*C;
for i=Q-50:Q-10
for j=2:630
if abs(G*j-i+H)/((G^2+1)^(1/2))<=0.5
J1(i,j,1)=255;
J1(i,j,2)=255;
J1(i,j,3)=255;
N2(i,j)=255;
P(i,j)=255;
end;
end;
end;%得到最底下的边框
for i=1:400
if N2(i,2)==255;
R=i;
end;
end;
c=0;
A=0;
B=0;
for i=R-160:R-120
for j=2:630
if N2(i,j)<90&(N2(i-2,j)-N2(i,j)>35)&(N2(i+2,j)-N2(i,j)>35)
c=c+1;
A=A+j;
B=B+i;
end;
end;
end;
C=A/c;
D=B/c;
E=0;
F=0;
for i=R-160:R-120
for j=2:640
if N2(i,j)<90&(N2(i-2,j)-N2(i,j)>35)&(N2(i+2,j)-N2(i,j)>35)
E=E+(j-C)*(i-D);
F=F+(j-C)^2;
end;
end;
end;
G=E/F;
H=D-G*C;
for i=R-160:R-120
for j=2:630
if abs(G*j-i+H)/((G^2+1)^(1/2))<=0.5
J1(i,j,1)=255;
J1(i,j,2)=255;
J1(i,j,3)=255;
N2(i,j)=255;
P(i,j)=255;
end;
end;
end;%得到最上面的边
for i=1:200
if N2(i,2)==255
S=i;
end;
end;
T=(R+S)/2;
T=fix(T);
for j=1:640
J1(T,j,1)=255;
J1(T,j,2)=255;
J1(T,j,3)=255;
N2(T,j)=255;
P(T,j)=255;
end;%得到中间的边
c=0;
A=0;
B=0;
for i=1:480
for j=1:640
if (I(i,j,1)-I(i,j,3)>25)&(I(i,j,1)-I(i,j,2)>25)
c=c+1;
A=A+j;
B=B+i;
end;
end;
end;
C=A/c;
D=B/c;
E=0;
F=0;
for i=1:480
for j=1:640
if (I(i,j,1)-I(i,j,3)>25)&(I(i,j,1)-I(i,j,2)>25)
E=E+(j-C)*(i-D);
F=F+(j-C)^2;
end;
end;
end;
G=E/F;
H=D-G*C;
for i=1:480
for j=10:630
if abs(G*j-i+H)/((G^2+1)^(1/2))<=0.5
if M2(i-25,j-2)-M2(i-25,j)>30&M2(i-25,j+2)-M2(i-25,j)>30&M2(i-25,j)<200&M2(i-26,j)<200&M2(i-24,j)<200&M2(i-25,j-1)<200&M2(i-25,j+1)<200&M2(i-26,j-1)<200&M2(i-26,j+1)<200&M2(i-24,j-1)<200&M2(i-24,j+1)<200
for k=i-200:i+50
J1(k,j,1)=255;
J1(k,j,2)=255;
J1(k,j,3)=255;
N2(k,j)=255;
P(k,j)=255;
end;
end;
end;
end;
end;%得到拟合的红线,与边框竖线的交叉点并上下延长
T=P;
h=1;
for i=1:Q
for j=10:630
if P(i,j)==255&P(i,j-1)==255&P(i,j+1)==255
p(h)=i;
h=h+1;
end;
end;
end;
U=min(p);
V=max(p);
for i=1:480
for j=1:640
if i<U|i>V
P(i,j)=0;
J1(i,j,1)=0;
J1(i,j,2)=0;
J1(i,j,3)=0;
end;
end;
end;
h=1;
for i=5:Q
for j=1:640
if T(i,j)==255&T(i-1,j)==255&T(i+1,j)==255
q(h)=j;
h=h+1;
end;
end;
end;
X=min(q);
Y=max(q);
for i=1:480
for j=1:640
if j<X|j>Y
P(i,j)=0;
J1(i,j,1)=0;
J1(i,j,2)=0;
J1(i,j,3)=0;
end;
end;
end;
for i=1:480
for j=1:640
if J1(i,j,1)==255|(J1(i,j,1)==0&J1(i,j,2)==0&J1(i,j,2)==0)
else J1(i,j,1)=0;
J1(i,j,2)=0;
J1(i,j,3)=0;
end;
end;
end;
imwrite(J1,'vv.bmp');%得到只包含分数栏和分数的图片
k=(T+S)/2;
k=fix(k);
h=1;
for j=1:640
if J1(k,j,1)==255&J1(k,j,2)==255&J1(k,j,3)==255&J1(k,j+1,1)==0&J1(k,j+1,2)==0&J1(k,j+1,3)==0
r(h)=j;
h=h+1;
end;
end;
r