Блог

Привет, я делаю свой проект, и я хочу помочь мне, пожалуйста, с построением функции излучения безмассовой частицы из круглого огненного шара без скорости. И сравните его с функцией излучения от огненного шара со скоростью vx = 0,5. если функция излучения S (x, k) = 1 / e ^ (Egamma — P ^ yU_y) / T, предполагая, что частицы представляют собой фотон массы и E = P ^ y = 1 ГэВ, T = 0,2GeV? Это то, с чего я начинаю, но я не знаю, как я могу задействовать функцию эмиссии,

 import matplotlib.pyplot as plt
import numpy as np
from numpy import ma
from matplotlib import ticker, cm

N = 100
x = np.linspace(-3.0, 3.0, N)
y = np.linspace(-2.0, 2.0, N)

X, Y = np.meshgrid(x, y)

# A low hump with a spike coming out.
# Needs to have z/colour axis on a log scale so we see both hump and spike.
# linear scale only shows the spike.
Z1 = np.exp(-X**2 - Y**2)
Z2 = np.exp(-(X * 10)**2 - (Y * 10)**2)
z = Z1   50 * Z2

# Put in some negative values (lower left corner) to cause trouble with logs:
z[:5, :5] = -1

# The following is not strictly essential, but it will eliminate
# a warning.  Comment it out to see the warning.
z = ma.masked_where(z <= 0, z)


# Automatic selection of levels works; setting the
# log locator tells contourf to use a log scale:
fig, ax = plt.subplots()
cs = ax.contourf(X, Y, z, locator=ticker.LogLocator(), cmap=cm.PuBu_r)


cbar = fig.colorbar(cs)

plt.show()