Объедините geom_function и geom_ribbon в ggplot2

Я хотел бы выделить область между двумя функциональными кривыми, построенными с помощью функции geom_function.

Например, я создаю две кривые:

ggplot() + geom_function(fun = function(x){x^2}) + geom_function(fun = function(x){x^2 + 1})

Как я могу выделить область между этими кривыми? Я думаю, что geom_ribbon может быть ответом, но весь пример кода, который я видел для geom_ribbon, использовал набор данных для сопоставления. Как я могу сделать это с помощью функции geom_function, то есть кривых, основанных только на функции, а не на фактических данных?