五点法求函数的导数

function df=FivePoint（funcx0typeh）
if nargin == 3 
    h = 0.1;
else if （nargin == 4 && h == 0.0）
        disp（‘h不能为0！‘）;
        return;
    end
end

y0 = subs（sym（func） findsym（sym（func））x0）;
y1 = subs（sym（func） findsym（sym（func））x0+h）;
y2 = subs（sym（func） findsym（sym（func））x0+2*h）;
y3 = subs（sym（func） findsym（sym（func））x0+3*h）;
y4 = subs（sym（func） findsym（sym（func））x0+4*h）;
y_1 = subs（sym（func） findsym（sym（func））x0-h）;
y_2 = subs（sym（func） findsym（sym（func））x0-2*h）;
y_3 = subs（sym（func） findsym