f:exp(-sqrt(x1*x1+x2*x2+x3*x3));
df1:diff(f,x1,1);
df2:diff(f,x2,1);
df3:diff(f,x3,1);
df11:diff(f,x1,2);
df12:diff(f,x1,1,x2,1);
df13:diff(f,x1,1,x3,1);
df22:diff(f,x2,2);
df23:diff(f,x2,1,x3,1);
df33:diff(f,x3,2);
df111:diff(f,x1,3);
df112:diff(f,x1,2,x2,1);
df113:diff(f,x1,2,x3,1);
df122:diff(f,x1,1,x2,2);
df123:diff(f,x1,1,x2,1,x3,1);
df133:diff(f,x1,1,x3,2);
df222:diff(f,x2,3);
df223:diff(f,x2,2,x3,1);
df233:diff(f,x2,1,x3,2);
df333:diff(f,x3,3);