X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fcomputational_physics.git;a=blobdiff_plain;f=zentral.c;fp=zentral.c;h=e7a9430014f8934e984a7203aa641ce4aec7e5a9;hp=cf071cd699b56a59a62f71240c71a2cbe70eaa78;hb=48ea798c247dfff5f6bcb8c6b51b115e829571f8;hpb=53bf1c9b8fa24799a5bb284c6e93e877dac813d0

diff --git a/zentral.c b/zentral.c
index cf071cd..e7a9430 100644
--- a/zentral.c
+++ b/zentral.c
@@ -46,7 +46,7 @@ int main(int argc,char **argv) {
 
  /* loop */
  for(i=0;i<steps;i++) {
-  r_3=x_p*x_p+y_p*y_p;
+  r_3=sqrt(x_p*x_p+y_p*y_p)*(x_p*x_p+y_p*y_p);
   f_x=-alpha*x_p/r_3; f_y=-alpha*y_p/r_3;
   x=x_p+vx_p*tau; y=y_p+vy_p*tau;
   vx=vx_p+f_x*tau; vy=vy_p+f_y*tau;