Question: Contour plot in polar coordinates

Hi, I'm trying to plot equipotential lines and streamlines in Maple, and both equations are in polar coordinates. I can plot a contourplot in cartesian coordinates, but adding "coords=polar", it throws an error. I tried to write a do loop but it wasn't having it either. any help would be appreciated

 

restart; with(LinearAlgebra); with(plots);
K := 1;
psi := K*r^(4/3)*sin(4/3*theta);
phi := K*r^(4/3)*cos(4/3*theta);

 

p1 := contourplot(psi, r = 0 .. 1, theta = 0 .. 3*Pi*(1/8));

p2 := contourplot(phi, r = 0 .. 1, theta = 0 .. 3*Pi*(1/8));
display(p1, p2);

 

Please Wait...