%方法1：
syms x y z;du_dx=diff((x^2+y^2+z^2)^(1/2),x)%給出x的偏導數

du_dy=diff((x^2+y^2+z^2)^(1/2),y)%給出y的偏導數

du_dz=diff((x^2+y^2+z^2)^(1/2),z)%給出z的偏導數

%方法2：本例也可以使用jacobian函數求偏導數
syms x y z;
jacobian((x^2+y^2+z^2)^(1/2),[x y z])