%pylab inline
from obliquity.distributions import Cosi_Distribution

cosi_dist = Cosi_Distribution((1.3,0.1),(15,0.3),(3.5,0.5))
cosi_dist.summary_plot()

cosi_dist_simple = Cosi_Distribution((1.3,0.1),(15,0.3),(3.5,0.5),veq_simple=True) #veq_simple gives faster calculation.

cosi_dist.veq_dist.plot(label='full')
cosi_dist_simple.veq_dist.plot(fig=0,label='simple (no diff. rot)',color='g')
legend()

cosi_dist.plot(label='full v_eq')
cosi_dist_simple.plot(fig=0,color='g',label='simple v_eq')
legend()