>> a=0;
>> for i=1:20
a=a+factorial(i);
end
>> format long
>> a
a =
2.561327494111820e+018