- printf("adding %f %f %f\n",force.x,force.y,force.z);
- printf("total i: %f %f %f\n",ai->f.x,ai->f.y,ai->f.z);
+ printf(" adding %f %f %f\n",force.x,force.y,force.z);
+ printf(" total i: %f %f %f\n",ai->f.x,ai->f.y,ai->f.z);
+ printf(" angle: %f\n",acos(cos_theta)*360.0/(2*M_PI));
+ printf(" d ij ik = %f %f\n",d_ij,d_ik);