如何理解三体问题中的语法错误?

时间:2019-05-08 08:35:00

标签: python python-3.x

我正在编写一个表示三个身体问题的Python代码:当重力公式输入时,问题就出现了:它说“语法无效”。这是代码段...

from pylab import *
"""
orbita vV variable tau
"""
Tmax=2.2; GM=4*pi**2
GMg=4*pi**2*9.54e-4; GMt=4*pi**2*3.015e-6          


x=1; y=0.; xg=5,26; yg=0; vx=0.; vy=pi
r=sqrt(x**2+y**2)   
rgt=sqrt((xg-x)**2+(yg-y)**2     
Fx=-GM*x/r**3; Fy = -GM*y/r**3                              
Fgtx=-GMg*(xg-x)/rgt**3; Fgty =-GMg*(yg-y)/rgt**3  

aa=0.1
v=abs(sqrt(vx**2 + vy**2))
a=abs(sqrt(Fx**2 + Fy**2)+abs(sqrt(Fgtx**2 + Fgty**2))
tau=(-v + sqrt(v**2 + 2*a*aa*abs(r)))/a 
File "<ipython-input-1-54f667811fbf>", line 21
tau=(-v + sqrt(v**2 + 2*a*aa*abs(r)))/a   
SyntaxError: invalid syntax

1 个答案:

答案 0 :(得分:1)

问题是前一行缺少)

a=abs(sqrt(Fx**2 + Fy**2)+abs(sqrt(Fgtx**2 + Fgty**2))

应该是:

a=abs(sqrt(Fx**2 + Fy**2))+abs(sqrt(Fgtx**2 + Fgty**2))

a=abs(sqrt(Fx**2 + Fy**2)+abs(sqrt(Fgtx**2 + Fgty**2)))