civ3's

数值微分.m
function g=dfs(f,x,h,s)
switch s
case 1,
g=(feval(f,x+h)-feval(f,x))/h;
case 2,
g=(feval(f,x+2*h)+feval(f,x)-2*feval(f,x+h))/(h*h);
otherwise,
end