How do I do an integration of a function of a dep-variable(with 2 arguments) which I have solved numerically

I have the numerical solution of Phi(r,t) which is a numerical solution of a pde..

I want to find the following energy integral at fixed time say t=10 ,

Energy = int(a*r^2*Phi(r,10)^2 + 2*b*r*Phi(r,10)^3,r=0..infinity)... I could not find help in maple.. I guess I had to do the following

1. select the numerical solution dependent only on r at t=10 slice..

2.Use numerical methods to evaluate the integral..

I don't know how to write code for this .. I would be grateful for any suggestions.

Thanks