파이썬 뉴턴 방법은

Question

내가 도움을 커뮤니티 많이 요청 된 발행하고 나는 그것을 모두 감사

그래서 파이썬에서하지만 누군가가 그 위에 볼 수 있었다 작동하지 않는 어떤 이유로 뉴턴 방법을 해결하는 프로그램에서 일하고 필자 부디? 당신에게 감사 =)

import sympy 
from collections import defaultdict 

def main(): 
    dir(sympy) 
    print ("NEWTONS METHOD") 
    print ("Write your expression in terms of 'x' ") 
    e = sympy.sympify(raw_input("input expression here: ")) 
    f = sympy.Symbol('x') 
    func1 = e 
    func1d = sympy.diff(e,f) #takes the dirivative of the function 
    print ("the dir of your function = "), func1d 
    x = input("number to substitute for x: ") 
    a = input("how many digits would you like to round to [recomended at least 4]") 
    func1sub = func1.subs({'x':x}) #substitutes the user given value of x into the equation 
    func1dsub = func1d.subs({'x':x}) #substitutes the user given value of x into the equation 
    func1sub = float(func1sub) 
    func1dsub = float(func1dsub) 
    func1sub = round(func1sub) 
    func1dsub = round(func1dsub) 
    round(func1sub,a) 
    round(func1dsub,a) 
    n = x - (func1sub/func1dsub) 
    x1 = 0 
    x2 = 0 
    n = x - (func1sub/func1dsub) 
    x1 = n 
    x1 = round(x1) 
    n = x2 - (func1sub/func1dsub) 
    x2 = n 
    x2 = round(x2) 
    while 0 == 0: 
     if abs(x1-x2) < .0001: 
     print x1 
     break 
     else: 
     n = x2 - (func1sub/func1dsub) 
     x2 = n 
     if abs(x - n) < .03: 
     print x 
    if func1dsub == 0: 
     print ("ERROR CAN NOT DIVIDE BY 0") 
main() 

Answer 1

당신은 여기에 무한 루프가 있어요 :

n = x2 - (func1sub/func1dsub) 
x2 = n 

그리고 당신의 루프 조건은 다음과 같습니다

while 0 == 0: 

    if abs(x1-x2) < .0001: 
     print x1 
     break 

    else: 
     n = x2 - (func1sub/func1dsub) 
     x2 = n 

    if abs(x - n) < .03: 
     print x 

중요한이 루프의 일부가 될 것으로 보인다 abs(x1-x2) < .0001이므로 다음을 다시 작성하겠습니다.

while abs(x1 - x2) >= .0001: 
    x2 -= (func1sub/func1dsub) 
print x1 

그래서 x2 -= (func1sub/func1dsub)은 x2을 잘못 입력하고 있습니다. 나는이 같은 인쇄 문을 추가하고 값이 실제로 수렴하고 있는지 확인 거라고 :

또한

while abs(x1 - x2) >= .0001: 
    x2 -= (func1sub/func1dsub) 
    print (x1, x2) 

, 나는 뉴턴의 방법이 익숙하지 해요,하지만 코드에서 절대 변경 func1sub/func1dsub하지만 안 모든 반복마다 변경 하시겠습니까?