`vecF=(d)/(dt)((m_0vecv)/(sqrt(1-v^2/c^2)))=m_0(vecv)/(sqrt(1-v^2/c^2))+m_0vecv/c^2vecv*vecv(1)/((1-v^2/c^2)^(3//2))`
Thus `vecF_(_|_)=m_0(vecw)/(sqrt(1-beta^2)), vecw=vecv, vecw_|_vecv`
`vecF_(||)=m_0(vecw)/((1-beta^2)^(3//2)), vecw=vecv, vecw_(||)vecv`