本文整理汇总了Python中planeta.Planeta.avanza_verletvelocity方法的典型用法代码示例。如果您正苦于以下问题：Python Planeta.avanza_verletvelocity方法的具体用法？Python Planeta.avanza_verletvelocity怎么用？Python Planeta.avanza_verletvelocity使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类planeta.Planeta的用法示例。
在下文中一共展示了Planeta.avanza_verletvelocity方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: range
# 需要导入模块: from planeta import Planeta [as 别名]
# 或者: from planeta.Planeta import avanza_verletvelocity [as 别名]
x = sp.zeros(pasos)
y = sp.zeros(pasos)
vx = sp.zeros(pasos)
vy = sp.zeros(pasos)
r = sp.zeros(pasos)
energia = sp.zeros(pasos)
perihelio = [[], [],[] ] # se trabajara con dt,dy y dx
''' Condiciones iniciales'''
[x[0],y[0],vx[0],vy[0]] = condicion_inicial
r[0] = sp.sqrt(x[0]**2+y[0]**2)
energia[0] = p.energia_total()
''' Creacion de arreglos de valores para [x,y,vx,vy] con metodo de verlet velocity'''
for i in range (1,pasos):
p.avanza_verletvelocity(dt)
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()
r[i] = sp.sqrt(x[i]**2+y[i]**2)
'''Tratamos de buscar el perihelio'''
b = 0.0009
valor =10
if valor-b<r[i-1] and r[i-1]<valor+b:
perihelio[0].append(p.t_actual)
perihelio[1].append(x[i-1])
perihelio[2].append(y[i-1])