remove redundent U in du
first try hmc in fortran; but did not work;
correct some errors , but still did not work