In [1]:
%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
from pymacula import MaculaModel
In [5]:
model = MaculaModel() # default
ts = np.arange(0,500,0.05)
plt.plot(ts, model(ts))
Out[5]:
You can also change the default parameters; e.g. longer rotation period, different numbers of spots, different maximum spot size, etc.:
In [6]:
model = MaculaModel(nspots=5, star={'Peq':50}, spots={'alpha_max':10})
plt.plot(ts, model(ts))
Out[6]:
In [ ]: