Question: Using numeric with dsolve and producing a polarplot...

Hi,

Quick question.  If I was to solve a second order ODE using the dsolve function, and I was to use the numeric option such that I get a Procedure as the result.  Lets say I run something like sol:=dsolve(...,numeric);  If theta is the parameter to sol, then I can obviously call sol(Pi) to return the solution at Pi, I will receive the following result [theta=Pi, r(theta)=..., D(r)(theta)=...].

Now, lets say I wanted a Polarplot of [theta,r(theta)], how can I acheive this?

I have tried several ways, my thinking was something like polarplot(rhs(evalf(sol(theta))[2],theta=0..2*Pi); but I get the following "Error, invalid input: rhs received eqn_n(theta)[2], which is not valid for its 1st argument, expr".

Would anyone of you clever people have any pointers on how to polarplot a numeric solution?

Thankyou.