斑竹和高手能帮菜鸟小弟吗，万分感激！

利郎

我这个程序是自己编的，因为是刚开始学习就自己编了这个BP网络。但是开始运行就不可以，为什么啊？
clear
p=0:0.001:1;
T=sin(10*p);
[w1,b1,w2,b2,w3,b3]=initff(p,2,'tansig',3,'purelin',2,'logsig');
err_goal=0.001;
max_epoch=500;
lr=0.05;
for epoch=1:max_epoch;
A1=tansig(w1*p,b1);
A2=purelin(w2*A1,b2);
A3=logsig(w3*A2,b3);
E=T-A3;
D3=deltalin(A3,E);
D2=deltalin(A2,D3,w3);
D1=deltatan(A1,D2,w2);
[dw1,db1]=learnbp(p,D1,lr);
[dw2,db2]=learnbp(A1,D2,lr);
[dw3,db3]=learnbp(A2,D3,lr);
w1=w1+dw1;
b1=b1+db1;
w2=w2+dw2;
b2=b2+db2;
w3=w3+dw3;
b3=b3+db3;
SSE=sumsqr(T-purelin(w2*tansig(w1*p,b1),b2)-logsig(w3*purelin(w2*A2,b2),b3));
if SSE>err_goal
   continue
else
    break
end
end
w1
b1
w2
b2
w3
b3