mathematica 函数求值问题

你问的不够具体。这样问的话我只能说方法有很多，这里就随便介绍几种了。
关于多元函数的每个变量的偏微分，说白了就是求梯度，Mathematica做偏微分的函数D有5种语法，其中第4种就是对付你这情况的。代数值嘛，可以定义函数再代：
f[x, y, z, w] = Log@x + y^2 + Sin[z] + Exp@w
df[x_, y_, z_, w_] = D[f[x, y, z, w], {{x, y, z, w}}]
(*第一种*)
df[1, 2, 3, 4]
(*第二种*)
df@@{1, 2, 3, 4}
(* 第三种 *)
df[{x_, y_, z_, w_}] = D[f[x, y, z, w], {{x, y, z, w}}]
df[{1, 2, 3, 4}]
df@{1, 2, 3, 4}
(*第四种，对于这种，定不定义函数都无所谓了*)
df[x, y, z, w] /. {x -> 1, y -> 2, z -> 3, w -> 4}
(*第五种？其实和第四种没啥区别*)
df[x, y, z, w] /. Thread[{x, y, z, w} -> {1, 2, 3, 4}]
就答到这吧。有什么不懂的就仔细查查自带帮助。