t=0:0.01:10;
y=randn(size(t));
yn=sin(y);

p=y(1:30);
t=yn(1:30);
[w1,b1,w2,b2]=initff(p,4,'tansig',1,'Purelin');
tp(1)=10;
tp(2)=1000;
tp(3)=0.01;
tp(4)=0.01;
tp=[tp(1) tp(2) tp(3) tp(4)];
[w1,b1,w2,b2,te,tr]=trainbp(w1,b1,'tansig',w2,b2,'Purelin',p,t,tp);
figure(1);
plot([0:1:te],tr);
a=simuff(p,w1,b1,'tansig',w2,b2,'Purelin');
figure(2);
plot([0:0.1:2.9],t,'+',[0:0.1:2.9],a,'-')