Расчет зависимости энергетической светимости абсолютно черного тела от температуры R(T). Python

[Utvhf]

Нужно написать программу на Pythone без использования numpy. Задание с приложением в файле.
from matplotlib import pyplot as plt
import math
c=0
m=500
v=3*10**8
h=6.63*10**-34
k=1.38*10**-23
T=500
listT=[]
listR=[]
while T<2000:
T=T+5
listT.append(T)
#print(listT)
for T in listT:
d=100*10**-6
h=5*10**-9
for i in range(0, 2*m+1):
l=c+h*i
y=(2*math.pi*(v**2)*h)/(l**5*(math.e**((h*v)/(l*k*T))-1))
integral=h/3*(4*np.sum(y[1::2])+2*np.sum(y[2:2*m-1:2])+y[0]+y[-1])
listR.append(integral)
print(listR)
plt.plot(listT, listR)
Вот мой код. Я не понимаю как его написать, так чтобы он работал