Hn=10;        %隱層神經元個數(Hidden number)
On=4;        %輸出層神經元個數(Output number)
fid = fopen('inputdata.txt');                
F= (fscanf(fid,'%g',[On 1]))'

load W;
for j=1:Hn                    
    temp2=Whi(:,j);        %取出Whi第j列，即第j個神經元所對應的權值
    net_h(j)=F*temp2+thh(j);                       
    Y(j)=f1(net_h(j));      %求出隱層的輸出
end         
for k=1:On
    temp3=Woh(:,k);         %取出Woh第k列，即第k個神經元所對應的權值
    net_o(k)=Y*temp3+tho(k);            
    Z(k)=f2(net_o(k));       %求出輸出層的輸出    
end
Z