threepoint2.m

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

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

for i=1:5
     y(i) = subs(sym(func), findsym(sym(func)), x0-2*h+i*h-h);
end
hd = 1/h/h;

switch type
    case 1,
        df = (y(3) - 2*y(4) + y(5))*hd;   %用第一个公式求二阶导数
       
    case 2,
        df = (y(2) - 2*y(3) + y(4))*hd;          %用第二个公式求二阶导数
        
    case 3,
        df = (y(1) - 2*y(2) + y(3))*hd;   %用第三个公式求二阶导数
         
end