﻿功能: 求函数三阶偏导值

格式: DiffThirdOrder(f,x)

f  : 符号函数对象
x  : 数值或者矩阵变量

说明:
1. 执行成功, 返回函数在x处对应的三阶偏导数值
2. 函数f里的变量个数必须为1
3. 本函数采用数值方法计算f在x处的导数值
4. 本函数主要用于无法符号求导的表达式,进行数值求导
5. 本函数采用4阶方法求解三阶偏导
$$
f^{(3)}(x)\approx \dfrac{-f(x+3h)+8f(x+2h)-13f(x+h)+13f(x-h)-8f(x-2h)+f(x-3h)}{8h^3}
$$

例子:
f = @sin(x);
t = 0:0.2:1;
df = DiffThirdOrder(f,t)//回车得到如下结果,可以对于 df与理论解-cos(t)
df =
[-0.99999999991183   -0.98006657404376   -0.92106099135714   -0.82533562270713   -0.69670670781888   -0.54030230846172 ]