ua and u are global variables
ua and u are global variables

They are declared in constants and defined in the main, so we don't need
to pass it as an argument of the subroutine.

cphyc/n-body,cphyc/n-body