function ax = fx(vx,vy)
  m = 1;
  b = 1;
  g = 10;
  vter = sqrt(m*g/b);
  V = sqrt(vx.^2+vy.^2);
  ax = -b/m*V.*vx;
endfunction