d
r
d
t
=
p
m
d
p
d
t
=
−
∇
[
V
(
r
)
+
W
(
r
,
E
)
]
∇
⋅
(
R
2
p
)
=
0