Colony

  float rayHits  (double Or_x,double  Or_y,double  Oc_xdouble ,double  Oc_y,double  Vr_x,double  Vr_y)

  /* Or_x = Origin of ray, x-coordinate

         Or_y = Origin of ray, y-coordinate

         Oc_x = Origin of circle, x-coordinate

         Oc_y = Origin of circle, y-coordinate

         Vr_x = Direction of ray, x

         Vr_y = Direction of ray, y

  */

  {

        if( ((2*Vr_x*(Or_x - Oc_x))^2 - 4*Vr_x*(Or_x - Oc_x)^2) < 0) // Discriminant < 0 

        {

                return -1;

        }

        if( ((2*Vr_x*(Or_x - Oc_x))^2 - 4*Vr_x*(Or_x - Oc_x)^2) = 0) // Discriminant = 0

        {

                return -((2*Vr_x*(Or_x - Oc_x))^2 - 4*Vr_x*(Or_x - Oc_x)^2) / (2*Vr_x));

        }

        \

  }