read "GeneralizedShepp.txt":
F:=f->subs(x=1,diff(f,x)):

#The values of M and P may be set to whatever the user desires.
#Be warned that M+P>2000 will likely take a long time to run.

M:=700:
P:=700:
OriginalShepp:=op(PGFTableGen(M,P,F,x)):
save OriginalShepp, "OriginalData.m":