Wondering about Domingo's Rational Mean book
Kirby Urner
urner at alumni.princeton.edu
Tue May 30 16:44:05 EDT 2000
Have folks on sci.math already discussed
http://www.etheron.net/usuarios/dgomez/Roots.htm
to anyone's recollection? I'd be interested in
comments and/or pointers to previous postings.
I don't have the necessary background to evaluate the
significance of this work, but the algorithms made
interesting programming exercises. Using them, I
was able to get:
3rd root of 2 is about:
65379522
---------
51891761
3rd root of 10 is about:
91969780593702397138462508494860
---------------------------------
42688590663356403236303435376201
You can push these things further. E.g.
10^(1/3) ~= 2196555918195778106276024585870817879185576760
----------------------------------------------
1019550942230358826790302872023722041913882251
I also liked his concise presentation of Halley's method, which
gives me the floating point nth rooth of an integer. Here's
a recursive version:
def halley(P,n,d=10,x=1):
# Source: http://www.etheron.net/usuarios/dgomez/Roots.htm
# P -- find nth root of this number
# n -- whole number root
# d -- level of depth for recursion (default 10)
# x -- initial value of x (default 1)
if d>1:
newx = ((n+1.0)*P*x + (n-1)*(x**(n+1)))/((n-1)*P + (n+1)*x**n)
return halley(P,n,d-1,newx)
else:
return x
Usage:
>>> halley(10,3) # 3rd root of 10 to depth 10
2.1544346900318838
>>> halley(2,3) # 3rd root of 10 to depth 10
1.2599210498948732
>>> halley(12931021,7) # 7th root of 12931021 to depth 10
10.336898680701067
>>> halley(12931021,7,40) # same to depth 40 (more accurate)
10.374031092075427
Sorry if this is old hat on sci.math -- thought it was pretty
cool myself.
Kirby
Cc: people at fluidiom.com (from whence I first learned of
Domingo's work):
=============
>From: "Robert Coulter" <rcoulter at mvrpc.org>
>To: <people at fluidiom.com>
>Subject: fluidiom: math stuff
>Date: Tue, 9 May 2000 16:55:14 -0400
>
>Anybody out there follow what this guy (Domingo Gomez Morin) is up to...?
>Kirby? Allen?
>Would this be useful in fluidiom ... *IS* this "New and Improved"?
>He doesn't seem bashful even as he bashes cartesian thought...
>so he oughta fit in OK around here.
>
>http://www.etheron.net/usuarios/dgomez/default.htm
>
>Core9
More information about the Python-list
mailing list