Ableitung von Funktionen mehrerer Veränderlicher

restart:
f:=(x,y,z)->             ;#(* Funktion eingeben *)
n:=        ; # (* n-te Ableitung *)
fn:=D[1$n,2$n,3$n](f);

#(* Beispiel *)
# Berechnen der 3.Ableitung an den Werten x=y=Pi/4 und z=2
f:=(x,y,z)-> sin(x)*cos(2*y)*z^4;
n:=3;
x1:=Pi/4:
y1:=Pi/4:
z1:=2:
fn:=D[1$n,2$n,3$n](f)(x1,y1,z1);