R语言求导数，并计算最后函数值

不得不说这个函数真的太管用了，如果参数过多，需要求导，需要验证可以用这些代码

#1、调包
library(symengine)

#2、定义函数参数为符号
x <- Symbol(“x”)
y <- Symbol(“y”)
theta <- Symbol(“theta”)

#3、写出函数定义
f<-x^2+2*(theta2)xy+3*y

#4、对函数求1阶导数
f1<- D(f, “theta”)

#5、令theta=5，计算出带给定theta的函数值
df_dtheta_value <- subs(f1, theta, 5)
df_dtheta_value

注：
如果想求二阶导数，可以用
f1<- D(f, “theta”)
df_dtheta_2nd <- D(f1, “theta”)