I have an integral with a change parameter (y) as follows.
y <- seq(1,2,0.01) integrand_1 <- function(t){ A1 <- 30*(y*30^(1/2)-t)^2 A2 <- (