- printf("potential energy: %f\n",u);
- printf("total energy (1): %f\n",e+u);
- printf("total energy (2): %f\n",get_total_energy(&md));
+ printf("potential energy: %.40f [J]\n",u);
+ printf("total energy (1): %.40f [J]\n",e+u);
+ printf("total energy (2): %.40f [J]\n",get_total_energy(&md));