Python中的牛顿方法

Question

我在python中编写一个程序，它将使用newtons方法解决零问题。 我完成了它的粗略版本，然后实现了几个不同的东西，并想知道我是否需要强制/改变它。 （该主题的知识可能会有所帮助）

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) 


   print ("the dirivative of your function = "), func1d

   x = input("number to substitude for x: ")


   func1sub = func1.subs({'x':x})
   func1dsub = func1d.subs({'x':x})
   n = x - float(func1sub/func1dsub)



   while n != x:
      func1sub = func1.subs({'x':x})
      func1dsub = func1d.subs({'x':x})
      n = x - float(func1sub/func1dsub)
      print n



main()

1）首先我想知道，因为n和x的值可能并不总是完全相同，我必须使用round函数进行舍入。

2）看了这个之后，我觉得我的while循环没有解决它应该解决的问题，它应该解决什么是x ，你可以在函数中将它插入x ，输出将是x 。 我可以通过将值附加到数组然后查看它们是一个数字的倍数实例来实现吗？

Answer 1

Q1：我同意。

Q2：在循环顶部附近x = n 。 你想继续前进。

Answer 2

首先， x的值在你的while循环中没有变化。

其次，终止标准应该是“当前”估计与“前一个”估计应该“足够接近”，即......你应该使用类似的东西：

while abs(current - previous) > tolerance: