#Computes machine epsilon.  If you try to find macheps by successively
# dividing by 2 and comparing to zero, the algorithm simply returns zero
macheps=function(){
me=1.0
while(1+me>1){
me=me/2.0
}
me
}