Igor Schein on Fri, 8 Aug 2003 02:14:45 -0400


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: extreme elllseries() behavior


On Sat, Aug 02, 2003 at 04:00:01PM +0200, Karim BELABAS wrote:
> On Fri, 1 Aug 2003, Igor Schein wrote:
> > ? forstep(k=23,24,.1,print1(k" ");print(elllseries(e,k)))
> > 23 0.9999999999999715783271077903
> > 23.10000000000000000000000000 0.9999999999999752574949227059
> > 23.20000000000000000000000000 0.9999999999999784603968424610
> > 23.30000000000000000000000000 0.9999999999999812486647844467
> > 23.40000000000000000000000000 0
> > 23.50000000000000000000000000
> >   ***   object too big, length can't fit in a codeword
> 
> And what is 'e' here ?

sorry, ellinit([0,0,1,0,10])

Using cut-off point !=1 allows to work around it, but the higher s is,
the further A needs to be from 1, slowing things down considerably, as
documented.  In the actual example, error, though too general, at least
prevents the disaster, while value of 0 at s=23.4 is unacceptable IMO.

Thanks

Igor