double add_function(double u1, double u2)
{ 
   return u1 + u2;
}