这是我的代码：
import numpy as np
import matplotlib.pyplot as plt
h=6.63e-34
k=1.38e-23
#h=6.62607015e-34 # Planck constant, units of J*s
#k=1.380649e-23 # Boltzmann constant, units of J/K
c=3.0e8 # speed of light M/s
T1=300
T2=1000
T3=10000
T4=3*10000
T5=100000
T6=3*100000
T7=1000000
T8=3*1000000
x1=np.linspace(10**7, 10**21, 99999)
x2=np.linspace(10**7, 10**21, 999999)
x3=np.linspace(10**7, 10**21, 999999999)
y1=((2 * h * x1**3) / c**2) * (1 / (np.exp((h * x1) / (k * T8))-1))
y2=((2 * h * x2**3) / c**2) * (1 / (np.exp((h * x2) / (k * T7))-1))
y3=((2 * h * x3**3) / c**2) * (1 / (np.exp((h * x3) / (k * T1))-1))
plt.loglog(x1,y1)
plt.loglog(x2,y2)
plt.loglog(x3,y3)
plt.ylim(10**(-15),10**1)
plt.show()
我想要用python把这个双log log图画出来，但是问题很大，x的取点数要很大才能把T1=300k画出来，画出来的图也不平滑，求助各位大佬指点一下
