vector3{
x,y,z
abs=√(x^2+y^2+z^2)
}

F,a,v,p　vector直交座標
p1,p2　質点

p=p2-p1
F=(p/pabs)(-GmM/pabs^2 )
a=F/m
v=v+at
p0=p0+vt