clc;
clear;
load P;
load T;
p1=P';
t1=T';
p11=p1(60:350);
t11=t1(60:350);
[Pn,minp,maxp,Tn,mint,maxt]=premnmx(p11,t11);
%Pn2=Pn(100:150);
%Tn2=Tn(100:150);
spread=7;%取0.2 ，0.5 ，0.75，1，1.5
net=newrbe(Pn,Tn,spread)
y=sim(net,Pn);
 %error=y-Tn;
 %res=norm(error);
    figure(1)
    plot(Pn,Tn,'-r');
    %title('函數的實際值');
    hold on
     plot(Pn,y);
    title('函數的實際值和函數的實際值對比圖');
    figure(2)
   plot(0:290,y-Tn);
    title('網絡的誤差曲線');