本文整理汇总了Python中planeta.Planeta.verlet_2方法的典型用法代码示例。如果您正苦于以下问题:Python Planeta.verlet_2方法的具体用法?Python Planeta.verlet_2怎么用?Python Planeta.verlet_2使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类planeta.Planeta
的用法示例。
在下文中一共展示了Planeta.verlet_2方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: rk4
# 需要导入模块: from planeta import Planeta [as 别名]
# 或者: from planeta.Planeta import verlet_2 [as 别名]
vy = sp.zeros(pasos)
energia=sp.zeros(pasos)
[x[0],y[0],vx[0],vy[0]] = condicion_inicial
energia[0]=p.energia_total()
''' Se define el segundo valor del vector yactual mediante rk4 (en este caso) puesto que verlet2 lo necesita para la iteracion'''
p.avanza_rk4(dt)
resultados = p.y_actual
x[1] = resultados[0]
y[1] = resultados[1]
vx[1] = resultados[2]
vy[1] = resultados[3]
''' Se obtienen los arreglos con todos sus valores mediante el metodo de verlet 2'''
for i in range (2,pasos):
p.verlet_2(dt,x[i-2],y[i-2])
resultados = p.y_actual
x[i] = resultados[0]
y[i] = resultados[1]
vx[i] = resultados[2]
vy[i] = resultados[3]
energia[i] = p.energia_total()
''' Se grafican los resultados, en este caso energia en funcion del tiempo y la trayectoria para el metodo de verlet 2 con alpha=0'''
plt.figure(2)
plt.subplot(2, 1, 1)
plt.subplots_adjust(hspace=.5)
plt.scatter(x,y,6,c=u'darkred',label = "trayectoria")
plt.title(" Solucion verlet 2 alfa = 0 ")
plt.xlabel("x [m]")
plt.ylabel("y [m]")