midpoint.m

来自「Matlab语言常用算法程序集.rar 里面有很多有用的Matlab数值运」· M 代码 · 共 12 行

function df=MidPoint(func,x0,h)
if nargin == 2 
    h = 0.1;
else if (nargin == 3 && h == 0.0)
        disp('h不能为0！');
        return;
    end
end

y1 = subs(sym(func), findsym(sym(func)),x0+h);
y2 = subs(sym(func), findsym(sym(func)),x0-h);
df = (y1-y2)/(2*h);