x <- roct()
exp(x+x) - exp(x)*exp(x) # zero to numerical precision
jj <- exp(log(x)/2) # use sqrt() here
jj*jj-x # also small
y <- roct()
exp(x+y) - exp(x)*exp(y) # some rules do not operate for onions
max(Mod(c(sin(asin(x))-x,asin(sin(x))-x))) # zero to numerical precision
Run the code above in your browser using DataLab