frame1.axhline(y=0, color='k', ls='--', alpha=.3)
t_sample = np.linspace(min(x), max(x), num=10001, endpoint=True)
# Planet 1 #
Planet1 = Model(P1=P1 / 100,
tau1=tau1 / 1000,
k1=k1 / 100,
w1=w1,
e1=e1,
P2=P2 / 100,
tau2=tau2 / 1000,
k2=0,
w2=w2,
e2=e2,
offset1=offset1,
offset2=0)
y1 = Planet1.get_value(t_sample)
plt.plot(t_sample, y1, 'b-.', alpha=.3, label='Planet1')
# Planet 2 #
Planet2 = Model(P1=P1 / 100,
tau1=tau1 / 1000,
k1=0,
w1=w1,
e1=e1,
P2=P2 / 100,
tau2=tau2 / 1000,
k2=k2 / 100,
w2=w2,
e2=e2,
offset1=0,
offset2=offset2)
y2 = Planet2.get_value(t_sample)
if 0:
# Planet 1 #
Planet1 = Model(P1=P1 / 100,
tau1=tau1 / 1000,
k1=k1 / 100,
w1=w1,
e1=e1,
P2=P2 / 100,
tau2=tau2 / 1000,
k2=0,
w2=w2,
e2=e2,
offset1=offset1,
offset2=0)
y1 = Planet1.get_value(t_sample)
plt.plot(t_sample, y1, 'b-.', alpha=.3, label='Planet1')
# Planet 2 #
Planet2 = Model(P1=P1 / 100,
tau1=tau1 / 1000,
k1=0,
w1=w1,
e1=e1,
P2=P2 / 100,
tau2=tau2 / 1000,
k2=k2 / 100,
w2=w2,
e2=e2,
offset1=0,
offset2=offset2)
y2 = Planet2.get_value(t_sample)