R code for intraocular lens calculations