Boa tarde caros membros

Alguém poderia auxiliar no seguinte código...

Tenho um modelo linear conforme expresso  no objeto"y". Neste código gostaria de obter uma integral "genérica" semelhante à derivada génerica obtida em D(y, "x1"). A idéia é integrar (dydx) em função de x1. x1 é o intervalo de no instante 1 e 2

Segue o código

###

y <- expression( a + b/x1 + c*log(x2) + d*((1/x1)*log(x2)) )
D(y, "x1")
dydx1 <- function(x, a=a, b=b, c=c, d=d)
{
 -(d * (1/x1^2 * log(x2)) + b/x1^2)
}
integrate(dydx, lower = x1 <- 1, upper = x1 <- 2) # Aqui que está o interesse em generalizar a integral independente dos valores dos parametros (a, b, c, d). O que se deseja é uma expressão.
### 

Obrigado

Samuel