Estou tendo alguns problemas com a função de integração em R. Estou tentando plotar o vo integral, mas parece que não estou fazendo corretamente.

t <- seq(0, 0.04, 0.0001)
vi <- function(x) {5 * sin(2 * pi * 50 * x)}
vo <- function(x) {integrate(vi, lower=0, upper=x)$value}

test_vect = Vectorize(vo, vectorize.args='x')
plot(t, vo(t))  # should be a cosine wave
plot(t, vi(t))  # sine wave

vo deve ser uma onda senoidal, mas usandotest_vect me dá enredo errado e usandovo dá diretamente erroos comprimentos 'x' e 'y' diferem. Alguém pode, por favor, me ajudar nesse assunto?