ym[i + 1] = ym[i] + vym[i] * dt + 0.5 * aym[i] * dt**2
A_xm, A_ym, rm = Acc_Grav_jup(xm[i + 1], ym[i + 1])
axm[i + 1] = A_xm
aym[i + 1] = A_ym
vxm[i + 1] = vxm[i] + 0.5 * dt * axm[i] + 0.5 * dt * axm[i + 1]
vym[i + 1] = vym[i] + 0.5 * dt * aym[i] + 0.5 * dt * aym[i + 1]
outf = open('prj5e_ear_mars_classtrial.txt', 'a')
outf.write("{:.6f}".format(x[i]))
outf.write("\t\t")
outf.write("{:.6f}".format(y[i]))
outf.write("\n")
outf.close()
plt.figure(1)
nam = Planet.name(p2)
plt.plot(x, y, label=nam)
print(xj, yj)
#plt.plot(xj, yj, label='jupiter')
#plt.plot(xm, ym, label='mars')
plt.xlabel("x[AU]")
plt.ylabel("y[AU]")
plt.legend(loc='upper right')
plt.title("solar sys model")
scatter([0], [0], s=20, color='yellow')
plt.savefig("solar sys model.png")
plt.show()
if __name__ == "__main__":
t2 = time.time()
print(t2 - t1)
plt.figure(1)
