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

格式: DiffFourthOrder(f,x)

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

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

例子:
f = @sin(x);
t = 0:0.2:1;
df = DiffFourthOrder(f,t)//回车得到如下结果,可以对于 df与理论解sin(t)
df =
[ 0.00000000000000    0.19866928313325    0.38941843614239    0.56464270630116    0.71735577734295    0.84147085149141 ]