用mathematica求复合函数高阶导数,需要怎么写程序啊?求求我的论文!

例如求y=x^4+x^3+x^2+a[x]的三阶导数
D[x^4+x^3+x^2+a[x],]
模块部分可以这样：
Dn[xx_]:=Module[,hanshu=Input[Inputhanshu];
f[x_]:=hanshu;
n=Input[please Input n的值];
x0=Input[please Input x0的值];
dn=D[f[xx],];
dn0=dn/.xx->x0;
}]
执行后输入下面命令
Dn[x]
在三次出现的对话框中依次输入 x^5 3 1
(分别表示函数y=x^5 阶数3 自变量的值1)
得到
Out[15]=}
把它放入某个软件包吧

已经很长时间没有弄这个了.我看暂时是不是这样：
Check[Print[...],]
用这个代替原来的Print[...]
看f[x_] = Input["input f[x]="]
少一个冒号