如何用R语言绘制变上限积分图?

R语言下,求y随着q变化的函数图,这个代码一直提示错误图片说明

y <- function(x,q){(x^2)*q}

f <- function(q){integrate(y,0,q)}

plot(f,0,100)

https://bbs.pinggu.org/thread-3379941-1-1.html