非線性方程式 exp(x)-3x^2=0
在 -0.5 與 -0.4之間有一根存在 若取 ε=0.01
Python程式
import math
#define math.exp(x) - 3.0*math.pow(x,2))
#====================================
def f_cal(x) :
temp= ( math.exp(x) - 3.0*x*x )
return temp
a=-0.5
b=-0.4
x=a
i=1
print(" i \t\tx \t f(x)\n")
while x<=b:
s = "{%5d} {%8.2f} {%2.5f}" %(i,x,f_cal(x))
print(s)
i+=1
x=x+0.01
======================================================
n >= (ln (b-a) -ln (ε)) / ln(2) =3.321 所以n=4
輸出畫面
======== RESTART: H:/2018-09勤益科大數值分析/數值分析/PYTHON/Pr2-1.py ===========
i x f(x)
{ 1} { -0.50} {-0.14347}
{ 2} { -0.49} {-0.10767}
{ 3} { -0.48} {-0.07242}
{ 4} { -0.47} {-0.03770}
{ 5} { -0.46} {-0.00352}
{ 6} { -0.45} {0.03013}
{ 7} { -0.44} {0.06324}
{ 8} { -0.43} {0.09581}
{ 9} { -0.42} {0.12785}
{ 10} { -0.41} {0.15935}
>>>
沒有留言:
張貼留言