Как написать интеграл для функции sp в pysynphot
В настоящее время я создаю спектры с помощью pysynphot, и я знаком с командой "sp(900)", которая возвращает значение потока фотонов на 900 нм.
Я раньше не писал никаких математических функций на Python, и мне любопытно, есть ли у кого-нибудь интегральная функция для команды flux? Я ищу значение интеграла между 900-1100 нм. Мой текущий код сценария:
import pysynphot as S
import matplotlib.pyplot as plt
from synphot import SourceSpectrum
from synphot.models import BlackBodyNorm1D
sp3 = S.Icat('ck04models', 13000, 0.0, 4.5)
plt.semilogy(sp3.wave, sp3.flux, 'r', label='B8')
plt.xlim(800, 2000)
plt.ylim(10e3, 10e11)
plt.xlabel(sp3.waveunits)
plt.ylabel(sp3.fluxunits)
plt.legend(loc='upper right')
sp(900)