%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
from pymacula import MaculaModel
MaculaModel
generates a model with 3 random spots (by default), using default stellar and spot parameters.
model = MaculaModel() # default
ts = np.arange(0,500,0.05)
plt.plot(ts, model(ts))
[<matplotlib.lines.Line2D at 0x7fe8ef675a50>]
You can also change the default parameters; e.g. longer rotation period, different numbers of spots, different maximum spot size, etc.:
model = MaculaModel(nspots=5, star={'Peq':50}, spots={'alpha_max':10})
plt.plot(ts, model(ts))
[<matplotlib.lines.Line2D at 0x7fe8ef675b90>]