R中函数D若使用自定义函数function如何实现求导？

比如有：
myfun <- function(x){
    2*x^2
}
D(expression(myfun), "x")却得不到求导的结果。如何处理？

我给的只是一个简单的自定义函数，所以可以使用
D(expression(2*x^2), "x")解决。但是实际使用中myfun自动的较为复杂，因此不适合直接在D中直接写，而需要先根据上边的方式自定义。用function编写函数后改怎么求偏导呢？

你可以试试从这个例子改

1. grad <- function(scal){
   
2.     scal <- substitute(scal)
   
3.     sapply(c('x', 'y', 'z'), function(v) D(scal, v))
   
4. }
   
5. grad(x*y*z)

复制代码