把代码改了一下:
# -*- coding: utf-8 -*-
import math
import sys
import time
import numpy as np
def area():
x1=1 #起始区间
x2=4.5 #结束区间
dx=0.0001 #步长,步长越小越精确
y1=0
y2=0
y3=0 #结果
for x in np.arange(x1, x2, dx):
y1=y1+dx*x*x #细小矩形取左边为高
y2=y2+ dx*(x+0.5*dx)*(x+0.5*dx) #取中间
y3=y3+dx*(x+dx)*(x+dx)
print (y1, y2, y3)
start = time.clock()
area()
end = time.clock()
print (end-start)