Vector vector_reflect(Vector L, Vector n) {

   Vector reflected_vector;

   reflected_vector = vector_subtract(L, vector_scale(n, 2.0*vector_dot_product(L, n)));

   return reflected_vector;
}