import numpy as np                       # llama a la libreria numpy (operaciones numericas)
import matplotlib.pyplot as plt          # llama a la libreria grafica

# Este codigo muestra los problenas de estabilidad de la integracion por Euler explícito.
#
# Vamos a definir una funcion maricial que realiza la siguente cuenta
#
# [0    1].[  y(t)  ]
# [-w^2 0] [dy(t)/dt]
#

tita0=0.1
omega2=2.0

n=1000
h=(2*np.pi)/n

y=np.linspace(tita0,tita0,n)
ydot=np.linspace(0,0,n)

i=1

while i<n:
	y[i]=y[i-1]+h*ydot[i-1]
	ydot[i]=-omega2*h*y[i-1]+ydot[i-1]
	i=i+1

plt.figure()                                    # inicio nuevo grafico
plt.plot(y,ydot)                                # grafico de diagrama de fases
plt.show()                                      # mostrar en pantalla