diffai.m
 程序代碼

function y3=diffai(x);
    y3=-0.75*sin(1.75*x)*exp(-x.^2/2)-cos(1.75*x)*exp(-x.^2/2)*x; 
fai.m
 程序代碼
function yl=fai(x)
    yl=cos(0.75.*x)+exp(-x.^2/2);


fnn.m
 程序代碼
function y2=fnn(x)
    y2=1/(1+exp(-x));