Question: Integrating numerical solution

sol := dsolve([ode2, theta(0) = (1/180)*Pi*10, (D(theta))(0) = 0], numeric)

sol(0)

[t = 0., theta(t) = .174532925199433, diff(theta(t), t) = 0.]

given a numerical solution in the above form, how can I integrate a function

f(theta(t)) numerically?

I tried to isolate theta(t) first:

x:=t->rhs(op(2,sol(t)))

But I can't integrate

phi:=t->int(x(t),t=0..t)

Please Wait...