\\ Problem 1 by the Ooura-Mori method read("ooura_mori.gp") read("lambertw.gp"); {func(x)=local(w,term); w=lambertw2(x); w/(1+w)/x } {p1(M)=q_ossinfM("func",M,1,1/2)} gettime; for (j=1,11, M=2^j; default(realprecision,M/2); print("M= ",M," p1(M)=",p1(M)); print("milliseconds=",gettime)); quit