n_f.m

追赶法

function y=n_f(a,n,x)
y=0.0；
for i=(n+1)(n+1)
y=y+a(i)*x^(n+1-i)；
end