alpha18O <- equilibrium(10)$alpha18O
fR      <- seq(0.9,1,0.001)
R0      <- 2005.2E-6
R       <-  R0*fR^(alpha18O-1)
plot(fR,R, type="l")