double modulo( double X, double Y ){
       double res;
       res = Y;
       for (;;){
           if (res == Y) break;
           res = res + Y;
       }
       return X - res;
}