SOR方法matlab实现问题

具体问题是http://www.bsuc.cn/jpkc/dcx/resources/method/2603.htm#1上的例子下面的是网页上的SOR代码，可是我用matlab运行时总是提示“missing variable or function” ，哪位高手能帮帮忙呢！
hx=11;hy=11; %设置网格节点数
v1=ones(hy,hx); %设置行列二维数组
m=10;n=10; %横纵向网格数
%上下两行的Dirichlet条件边界值：
v1(1,:)=ones(1,hx)*50;
v1(hy,:)=ones(1,hx)*100;
%左右两列的Dirichlet条件边界值：
for i=1:hy
v1(i,1)=0;
v1(i,hx)=100;
end
%计算松弛因子
t1=(cos(pi/m)+cos(pi/n))/2;
w=2/(1+sqrt(1-t1*t1));
v2=v1;maxt=1;t=0; %初始化
k=0
while(maxt>1e-6) %由v1迭代，算出v2，迭代精度为0.000001
k=k+1 %计算迭代次数
maxt=0;
for i=2:hy-1 %从2到hy-1行循环
for j=2:hx-1 %从2到hx-1列循环
v2(i,j)=v1(i,j)+(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1)-4*v1(i,j))*w/4; %拉普拉斯方程差分式
t=abs(v2(i,j)-v1(i,j));
if(t>maxt) maxt=t;end
end
end
v1

冰火sor的问题 用MATLAB编程实现如下问题 关于sor技能加成问题 ? 用MATLAB实现克隆选择算法的布局问题 用matlab实现语音通信仿真的问题！ 如何用MATLAB实现？？？？？？？ matlab实现FT Matlab实现一个函数 暗黑SOR加点的问题．． matlab问题