matlab符号函数怎么求阶乘

syms n;
b=sym('(n+1)!');
c=symsum((-1)^nb/n^(n+1),n,1,inf)%显然求不出来
%但可以求数值
vpa(c)

a=input('请输入一个数：');
while a<0 || fix(a)~=a
disp('您输入的数不是正整数')
a=input('请重新输入：');
end
a
disp('a!=')
disp(factorial(a))
（运行后：）
请输入一个数：-20
您输入的数不是正整数
请重新输入：-12
您输入的数不是正整数
请重新输入：23
您输入的数不是正整数
请重新输入：-43
您输入的数不是正整数
请重新输入：12
a =
12
a!=
479001600

prod()这个命令对所有数字相乘，当然，两个阶乘相乘能用一个prod表示，如：
>> prod([1:5])
ans =
120
>> prod([prod(1:3) prod(1:5)])
ans =
720

---