Почему вы не печатаете значения во время каждой итерации? Вы увидите, что происходит, затем вы можете попробовать отладить его и переосмыслить свое решение
Просто сделайте:
def cube_root(num):
import random
g = random.randint(1,num+1)
new_g = float(num/(g*g))
print("g", g, "new_g", new_g, "num", num)
check_ifclose(g,new_g,num)
def check_ifclose(g,new_g,num):
from math import isclose
print("g", g, "new_g", new_g, "num", num)
if isclose(g, new_g, abs_tol=0.0000001):
print("The cube root of ",num,"is ",float(new_g))
else:
g = float(new_g)
new_g = float(num/(g*g))
check_ifclose(g,new_g,num)
Тогда для cube_root(16)
вы получите вывод:
>>> cube_root(16)
g 2 new_g 4.0 num 16
g 2 new_g 4.0 num 16
g 4.0 new_g 1.0 num 16
g 1.0 new_g 16.0 num 16
g 16.0 new_g 0.0625 num 16
g 0.0625 new_g 4096.0 num 16
g 4096.0 new_g 9.5367431640625e-07 num 16
g 9.5367431640625e-07 new_g 17592186044416.0 num 16
g 17592186044416.0 new_g 5.169878828456423e-26 num 16
g 5.169878828456423e-26 new_g 5.986310706507379e+51 num 16
g 5.986310706507379e+51 new_g 4.464794497196387e-103 num 16
g 4.464794497196387e-103 new_g 8.02633041618099e+205 num 16
g 8.02633041618099e+205 new_g 0.0 num 16
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "<stdin>", line 6, in cube_root
File "<stdin>", line 9, in check_ifclose
File "<stdin>", line 9, in check_ifclose
File "<stdin>", line 9, in check_ifclose
[Previous line repeated 8 more times]
File "<stdin>", line 8, in check_ifclose
ZeroDivisionError: float division by zero
Как видите, это похоже на то, как Размик сказал о небольших числах.
Но почему?
Потому что ваше решение ошибочно. Вы уверены, что должны рассчитать new_g
таким образом? Если вы посмотрите здесь, он становится слабым - очень маленьким, затем очень большим, затем снова маленьким ... Почему бы вам не попробовать, например, установить new_g как среднее от старого и текущего нового, чтобы сбалансировать его?
РЕДАКТИРОВАТЬ: решил проверить мое последнее предложение, теперь он отлично работает.
def cube_root(num):
import random
g = random.randint(1,num+1)
new_g = float(num/(g*g))
print("g", g, "new_g", new_g, "num", num)
check_ifclose(g,new_g,num)
def check_ifclose(g,new_g,num):
from math import isclose
print("g", g, "new_g", new_g, "num", num)
if isclose(g, new_g, abs_tol=0.0000001):
print("The cube root of ",num,"is ",float(new_g))
else:
g = float(new_g)
new_g = (g+float(num/(g*g)))/2
check_ifclose(g,new_g,num)
>>> cube_root(16)
g 3 new_g 1.7777777777777777 num 16
g 3 new_g 1.7777777777777777 num 16
g 1.7777777777777777 new_g 3.420138888888889 num 16
g 3.420138888888889 new_g 2.3939850310351836 num 16
g 2.3939850310351836 new_g 2.592869433541053 num 16
g 2.592869433541053 new_g 2.4863846555355664 num 16
g 2.4863846555355664 new_g 2.537249169224923 num 16
g 2.537249169224923 new_g 2.51131729000868 num 16
g 2.51131729000868 new_g 2.524147960725603 num 16
g 2.524147960725603 new_g 2.517700180896523 num 16
g 2.517700180896523 new_g 2.52091579334918 num 16
g 2.52091579334918 new_g 2.519305938864469 num 16
g 2.519305938864469 new_g 2.520110351423893 num 16
g 2.520110351423893 new_g 2.51970801680198 num 16
g 2.51970801680198 new_g 2.5199091519863974 num 16
g 2.5199091519863974 new_g 2.5198085763676823 num 16
g 2.5198085763676823 new_g 2.519858862169772 num 16
g 2.519858862169772 new_g 2.519833718766991 num 16
g 2.519833718766991 new_g 2.5198462903429375 num 16
g 2.5198462903429375 new_g 2.519840004523604 num 16
g 2.519840004523604 new_g 2.519843147425431 num 16
g 2.519843147425431 new_g 2.519841575972557 num 16
g 2.519841575972557 new_g 2.519842361698504 num 16
g 2.519842361698504 new_g 2.519841968835408 num 16
g 2.519841968835408 new_g 2.5198421652669256 num 16
g 2.5198421652669256 new_g 2.519842067051159 num 16
The cube root of 16 is 2.519842067051159