@Carl Love Thanks for this, the issue persists unfortunately. Also, this new code doesn't seem to run for me, I get the error

Error, (in LinearAlgebra:-Norm) expects its 1st argument, A, to be of type {Matrix, Vector}, but received 6.

I'm running Maple 2019 on Windows 10.

@Carl Love Thanks for this! This is much cleaner. There was one adjustment that I forgot to metnion: if I want to add a global piece of text which is always at the top of Graphs.txt when I run the code, (say, another 5 lines of code that will not depend on n) how would I do this? Thanks again for your help.

EDIT: I also noticed that the code is not working for me for n>5. Maple gives the error

Error, (in unknown) Maple was unable to allocate enough memory to complete this computation.  Please see ?alloc

Is this as issue with Maple itself or my computer? If I have several processes (Adobe, Firefox windows etc.) running in the background can this effect the performamce?

@mmcdara Thank you, this is quite helpful! I am very new to Maple, so I appreciate your help. Is there a way I can automate this procedure for other polynomials B? I gave B(4) in my original post, but I have defined the following:

for n from 0 to 10 do
    d(n) := factor(expand(sum((-1)^k*k!*Stirling2(n, k)/binomial(k + p + 1, k), k = 0 .. n)));
end do;
for n from 0 to 10 do
    B(n) := collect(sum(binomial(n, k)*d(k)*x^(n - k), k = 0 .. n), x, factor);
end do;

How can I get I automatically produce a plot (with p<5000) for each B(n)?