ส่งงาน week12

week12/Initial-Value.py

@@ -27,13 +27,17 @@ def bracket(f, **kwargs):
     xi -= h
     print('Solution xi',xi,f(xi))
 bracket(f,a=-2,b=2,h=0.01)
+#Output Solution xi 1.0000000000000022 0.0
 
 #Golden Section Search
+
 import math
+from math import log as ln
 R = (-1 + 5**0.5)/2
-def n(a=-2,b=2,ebs=10**-9) :
+eps=10**-9
+def n(a=-2,b=2,eps=10**-9) :
     from math import log as ln
-    return int(-ln(R)*ln(eps/abs(b-a)))
+    return int(ceil(-2.078087)*ln(eps/abs(b-a)))
 def search(f,a,b,tol=1.0e-9):
     nIter = n(a,b)
     0.618033989
@@ -55,3 +59,4 @@ def search(f,a,b,tol=1.0e-9):
     if f1 < f2: return x1,f1
     return x2,f2
 search(f,-2,2)
+#Output (0.9999999940181203, 0.0)