请帮忙分析调试这段程序：
Ir=0.05;
err_goal=0.001;
max_epoch=10000;
X=[1 2;-1 1;1 3];
T=[1 1;1 1];
[M,N]=size(X);q=10;[L,N]=size(T);
Wij=rand(q,M);
Wki=rand(L,q);
b1=zeros(q,1);b2=zeros(L,1);
for epoch=1:max_epoch
Oi=tansig(Wij*X,b1);
Ok=purelin(Wki*Oi,b2);
E=T-Ok;
deltak=deltalin(Ok,E);
deltai=deltalog(Oi,deltak,Wki);
[dWki,db2]=learnbp(Oi,deltak,Ir);
Wki=Wki+dWki;b2=b2+db2;
[dWij,db1]=learnbp(X,deltak,Ir);
Wij=Wij+dWij;
b1=b1+db1;
SSE=sumsqr(T-purelin(Wki*tansig(Wij*X,b1),b2));
if(SSE<err_goal)break;end
end
epoch
X1=X;
Oi=tansig(Wij*X1,b1);
Ok=purelin(Wki*Oi,b2);
Ok

注：q=10不能更改

请各位高手赐教！
邮箱：cat-1218@163.com，qq：849624371