二分法で三乗根の近似値を求めるコード

Question

x＜0の場合も含めた、任意の実数値xに対して、三乗根の近似値を求めるコードを完成させるという問題なのですが、自分の書いたプログラムだと負の数の三乗根が求められません。
恐らく、探索区間が違うと思うのですが、どこが間違っているのでしょうか。

def cube_root_bisection_full_range(cube):
     pass

     epsilon = 0.000000001
     count_guesses = 0

#探索区間の上限と下限(初期値 (0、cube))

     low =0.0
     if cube >= 1:
        high = cube
     else:
        high = 1

     low < 0.0
     if cube >= 1:
        high = cube
     else:
        high < 1

     guess = (low + high)/2
     while abs(guess**3 - cube) >= epsilon  :  
        
        if guess**3 < cube:
            low = guess  
        else:
            high = guess 
        
        guess = (low + high)/2
        count_guesses += 1
        
      
        if count_guesses >1e6:
            print(f'FAIL: cannot reach the cube root {cube}')
            break
    
    　print(' count_guess =', count_guesses )
   　 if abs(guess**3 - cube) < epsilon:
        　print(guess, 'is close to the cube root of ',cube)
        　print(f'verify: {guess}**3 =',guess**3 )
   　 print('')

cube_root_bisection_full_range(0.1)    
cube_root_bisection_full_range(-0.5)   

cube_root_bisection_full_range(-2)  
cube_root_bisection_full_range(5)  
    

【結果】
count_guess = 28
0.464158883318305 is close to the cube root of 0.1
verify: 0.464158883318305**3 = 0.09999999997222532

FAIL: cannot reach the cube root -0.5
count_guess = 1000001

FAIL: cannot reach the cube root -2
count_guess = 1000001

count_guess = 32
1.7099759465781972 is close to the cube root of 5
verify: 1.7099759465781972**3 = 4.999999999135954

Comments

[Zuishin]

ループが無いようですが。

[Zuishin]

Python ってこんないい加減なインデントで動きましたっけ？

[TakaiY]

すでに指摘されているように、必要な部分が無いし、また、関数名がおかしかったり、インデントがおかしかったり、全角文字が混っていたり、処理内容も的と、これは実際に動いているソースではないようなので、これを見て何かアドバイスしても無駄になりそうです。

実際に動かしているソースを提示いただけますか。

Answer 1

def cube_root_bisection_full_range(cube):
    epsilon = 1e-9
    count_guesses = 0

    if abs(cube) < epsilon:
        print('0.0 is close to the cube root of ', cube)
        return

    #探索区間の上限と下限(初期値 (0、cube))
    if cube >= 1.0:
        low = 1.0
        high = cube
    elif 0.0 <= cube < 1.0:
        low = cube
        high = 1.0
    elif -1.0 < cube < 0.0:
        low = -1.0
        high = cube
    else:
        low = cube
        high = -1.0

    guess = (low + high)/2
    while abs(guess**3 - cube) >= epsilon  :  
        if guess**3 < cube:
            low = guess  
        else:
            high = guess 
        
        guess = (low + high)/2
        count_guesses += 1
      
        if count_guesses > 1e6:
            print(f'FAIL: cannot reach the cube root {cube}')
            break
    
    print(' count_guess =', count_guesses )
    if abs(guess**3 - cube) < epsilon:
        print(guess, 'is close to the cube root of ', cube)
        print(f'verify: {guess}**3 =', guess**3 )
    print('')

if __name__ == '__main__':
    cube_root_bisection_full_range(0.1)    
    cube_root_bisection_full_range(-0.5)   

    cube_root_bisection_full_range(-2)  
    cube_root_bisection_full_range(5)  
    cube_root_bisection_full_range(0.0000000001)  

 count_guess = 26
0.4641588844358921 is close to the cube root of  0.1
verify: 0.4641588844358921**3 = 0.10000000069455584

 count_guess = 27
-0.7937005255371332 is close to the cube root of  -0.5
verify: -0.7937005255371332**3 = -0.4999999991552862

 count_guess = 28
-1.2599210496991873 is close to the cube root of  -2
verify: -1.2599210496991873**3 = -1.999999999068104

 count_guess = 29
1.7099759466946125 is close to the cube root of  5
verify: 1.7099759466946125**3 = 5.000000000157156

0.0 is close to the cube root of  1e-10

Comments

[hagi]

elifを使えばよかったんですね！
ありがとうございました！