08/29/2007, 06:19 AM

jaydfox Wrote:Using gp, I computed Andrew's solution for base e, using a 50x50 matrix. (Side question: Andrew, how much faster are other libraries at solving these large matrices?)Jay,

could you provide me your gp-code for comparision? You may put it here or send me an email.

Also I would like to check, whether your last computations can be smoothed by accelerating convergence using Euler-summation of terms. But I've only Pari/GP and maxima(the latter with extremely little experience) and I don't know how complex your computations were and how much work it would be to get it implemented.

What do you think?

Gottfried

Gottfried Helms, Kassel