本文整理匯總了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])