for i=1:5
net=newff(p,t,3,{'logsig'},'trainrp');
net.trainParam.epochs =1000;
net.trainParam.max_fail=10;
net.trainParam.show=50;
net.trainParam.Ir=0.05;
net.trainParam.goal=10^-6;
net=train(net,p,t);
y=sim(net,pp);
e=tt-y;
performance(i ,: )=mse(e);
end
performance;
perf=min(performance)
dif_performance=max(performance)-min(performance)
plot(performance);