function sdot=f(stato,t)
  g = 9.8;
  x  = stato(1);
  y  = stato(2);
  vx = stato(3);
  vy = stato(4);
  
  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); 
  
  sdot(1) = vx;
  sdot(2) = vy;
  sdot(3) = ax;
  sdot(4) = ay;
endfunction
  
  g = 9.8;
  t = 0:0.01:10;
  dt=t(2)-t(1);

  posizione_x = [];
  posizione_y = [];
  s(1) = 5;
  s(2) = 0;
  s(3) = 0;
  s(4) = 0;


  for now=t
    posizione_x=[posizione_x s(1)];
    posizione_y=[posizione_y s(2)];
    s1 = f(s,t);
    s(1) = s(1) + s1(1) * dt;
    s(2) = s(2) + s1(2) * dt;
    s(3) = s(3) + s1(3) * dt;
    s(4) = s(4) + s1(4) * dt;
  end;
