Thank you so much for all of your help! I went with: f <- function(x) body(f) <- D(expression(x^2),"x") This is SO COOL! Thanks again, Erin mailto: hodgess at gator.uhd.edu