In [1]:
import batman
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
In [2]:
params = batman.TransitParams()
params.t0 = 0.0
params.per = 1.05
params.rp = 0.1
params.a = 15.0
params.inc = 88.0
params.ecc = 0.05
params.w = 90.0
params.u = [0.1, 0.3]
params.limb_dark = "quadratic"
In [3]:
t = np.linspace(-0.05, 0.05, 100)
m = batman.TransitModel(params, t)
flux = m.light_curve(params)
In [4]:
plt.plot(t, flux)
plt.xlabel("time from central transit")
plt.ylabel("relative flux")
plt.ylim(0.985, 1.005)
plt.show()
In [ ]: