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