#!/usr/bin/env python
# coding: utf-8

# # **Tarea No. 1 - Física II**
# 
# 
# ***

# # Ecuación de Kepler

# In[3]:


from IPython.display import YouTubeVideo
YouTubeVideo('KmzFDEu2RoA')


# $$ M = E - e\sin(E) $$

# # Solución por el método de Newton - Raphson

# In[12]:


from IPython.display import HTML
HTML('<iframe src=https://es.wikipedia.org/wiki/M%C3%A9todo_de_Newton width=800px height=400px >')


# In[11]:


from mpl_toolkits.mplot3d import Axes3D
from matplotlib.mlab import griddata
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter
import matplotlib.pyplot as plt
import numpy as np
import sys
import math

M = 1.829166028
e = 0.093
E = M

for i in range(10):
    E = E - ((E-e*math.sin(E)- M)/(1-e*math.cos(E)))
E = E*180/3.14159
print("la anomalia excentrica es")
print(E)


# In[ ]: