Построение графиков с похожими кодами - PullRequest
0 голосов
/ 07 августа 2020

Я построил свой график с помощью этого кода:

from numpy import*
from matplotlib.pyplot import*
h=6.626*10**(-34)
c=3*10**8
k=1.38*10**(-23)
t=6000
l=[]
s=arange(100,2000,1)
def fun(x,t):
    e=(2*pi*h*c**2)/(x**5*(exp((h*c)/(x*k*t))-1))
    return e
for x in arange(100*10**(-9),2000*10**(-9),1*10**(-9)):
    r=fun(x,t)
    l.append(r)
plot(s,l)
show()

Любезность RandomGuy, он дал очень компактный код:

from numpy import*
from matplotlib.pyplot import*
h=6.626*10**(-34)
c=3*10**8
T=6000
k=1.38*10**(-23)
l=linspace(100*10**(-9),2001*10**(-9),100)
E=(2*pi*h*c**2)/((l**5)*exp(h*c/(l*k*T)-1))
plot(l,E)
show

Но проблема в том, что я получаю разные максимумы для пиков. В первом он около 1, а во втором больше 2,5. Чем отличаются эти коды?

graph of first code

график второго кода

1 Ответ

3 голосов
/ 07 августа 2020

Разница происходит из аргумента экспоненты.

Во второй функции -1 находится в аргументе exp, что не относится к первой.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...