Clase del día 29 de Junio de 2021 parte 2

import scipy.integrate as integrate import matplotlib.pyplot as plt import numpy as np Las condiciones iniciales que se puede presentar en el sistema serian, t = 0 , x = x 0 , y x = 0 . Por otra parte la otra combinacion seria t = 0 , x = 0 y x ˙ = v 0 Aplicamos la primera combinacion de condiciones iniciales en la ecuacion de x teniendo: x 0 = A ∗ ( 0 ) + B ∗ ( 1 ) B = x 0 Ahora sustituimos en x ˙ 0 = A ∗ w n ∗ ( 1 ) A = 0 Definimos entonces que: x 1 = x 0 c o s ( w n t ) Continuamos sustituyendo las condiciones iniciales restantes 0 = A ∗ ( 0 ) + B ∗ ( 1 ) B = 0 Sustituyendo en x ˙ v 0 = A ∗ w n ∗ ( 0 ) A = v 0 w n De ahi tenemos que x 2 = v 0 w n s i n ( w n t ) Entonces x = x 1 + x 2 Desarrollando x = v 0 w n s i n ( w n t ) + x 0 c o s ( w n t ) pi=np.pi sqrt=np.sqrt cos=np.cos sin=np.sin def deriv_x ( d2x , phi ): x, xdot=d2x return [xdot...