homework_2018_1_09_1035441007.ipynb

Kernel: Python 3

from scipy import integrate
import numpy as np
import matplotlib.pyplot as plt
%pylab inline

Populating the interactive namespace from numpy and matplotlib

#Defino mis parámetros y cndiciones iniciales
g=-9.8 #[m/s**2]
y0=200 #[m]
a=0 #Límite inferior del intervalo de integración
b=3 #Límite superior del intervalo de integración

#Definimos la función a la que equivale mi EDO
def f(y, t):
    return g * t

#Creamos t's que prtenezcan al intervalo [0,3]
t=np.arange(a, b, 0.01)

#Integrando

y= integrate.odeint(f, y0, t)

plt.plot(t,y)
plt.ylabel("Posición [m]")
plt.xlabel("Tiempo [s]")
plt.title("Posición vs. Tiempo: Aplicación a las EDO")
plt.show()