TREEROOT.cpp

div(J) + dp/dt = 0