Как написать интеграл для функции sp в pysynphot - PullRequest
Новая
       59

Как написать интеграл для функции sp в pysynphot

0 голосов
/ 13 февраля 2020

В настоящее время я создаю спектры с помощью pysynphot и знаком с командой sp (900), которая возвращает значение потока фотонов при 900 нм.

Я не писал никаких математических функций в python раньше, и мне любопытно, включил ли кто-нибудь интегральную функцию для команды потока? Я ищу, чтобы получить значение для интеграла между 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)
...