例題4-1 計算 exp(x) 在[0 , 1]的定積分

例題4-1 計算 exp(x) 在[0 , 1]的定積分

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

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: F:\2018-09勤益科大數值分析\數值分析\PYTHON\EX4-1.py =============
T{10}={  1.719713}

>>>