已知用de45求解变刚度阻尼双足机器人动力学方程，用syms求某个变量对所有变量求导MATLA

B代码如下：

function f = example(x)
syms a b c d e
f = diff(a*x^4 + b*x^3 + c*x^2 + d*x + e, x);
end

其中，`syms a b c d e` 定义了变量a、b、c、d、e为符号变量；`diff()` 函数用于求导数。

在命令行窗口中输入 `x = sym('x')`，定义变量x为符号变量。然后输入 `example(x)`，即可得到 `4*a*x^3+3*b*x^2+2*c*x+d`，即对 `a*x^4 + b*x^3 + c*x^2 + d*x + e` 对x求导的结果。