Is there an R package that fits multiple functions to find the best fit between two variables?

Question

I have two variables, x and y, and I would like to find the function that best describes how y varies as a function of x (y~x).

I know I can visually look at the x,y