範例EX4-1 利用梯行積分法求方程式近似解---Python語言

範例EX4-1 利用梯行積分法求方程式近似解---Python語言

# 利用梯行積分法求方程式近似解

from math import *

def F(x):
    temp=(exp(x))
    return temp

n=10
a=0.0
b=1.0
sum1=0.0
h=(b-a)/n
x=a
for i in range (1,n):
    x=x+h
    sum1=sum1+F(abs(x))
    tn=(h/2.0)*(F(fabs(a))+F(fabs(b))+2.0*sum1)

print("T{%d}={%10.6f}\n" %(n,tn))

輸出畫面
======= RESTART: H:/2018-09勤益科大數值分析/數值分析/PYTHON/EX4-1.py ==============
T{10}={  1.719713}

>>>