%
%            Function baseball_eval.m
%

  function F = baseball_eval(t,y,empty,k_D,k_L,g,omega)

  y1 = y(1);
  y2 = y(2);
  y3 = y(3);
  y4 = y(4);
  y5 = y(5);
  y6 = y(6);

  wx = omega(1);
  wy = omega(2);
  wz = omega(3);

  v = sqrt(y2^2 + y4^2 + y6^2);

  F(1,1) = y2;
  F(2,1) = -k_D*v*y2 + k_L*(wy*y6 - wz*y4);
  F(3,1) = y4;
  F(4,1) = -k_D*v*y4 + k_L*(wz*y2 - wx*y6);
  F(5,1) = y6;
  F(6,1) = -k_D*v*y6 + k_L*(wx*y4 - wy*y2) - g;