  g = 9.8;
  t=0:0.01:10;
  dt=t(2)-t(1);

  posizione_x = [];
  posizione_y = [];
  vx = 0;
  vy = 0;
  x  = 5;
  y  = 0;

  for now=t
    posizione_x=[posizione_x x];
    posizione_y=[posizione_y y];
    x  = x + vx * dt;
    y  = y + vy * dt;
    ay = - g * (x*x)/(x*x + y * y) - (vx * vx + vy * vy) * y / (x * x + y * y);
    ax =   g * (x*y)/(x*x + y * y) - (vx * vx + vy * vy) * x / (x * x + y * y); 
    vy = vy + ay * dt;
    vx = vx + ax * dt;
  end;
