T=[0.6 1];
net=newlin(minmax(P),1,0,0.01);
net=init(net);
net.trainParam.epochs=500;
net=train(net,P,T);
R=sim(net,P);
plot(P,R);
hold on
plot(P,T,'r+');
E=T-R;
SSE=sumsqr(E);
w_range=-1:0.1:1;
b_range=-1:0.1:1;
ES=errsurf(P,T,w_range,b_range,'purelin');