I’m realizing now that the solution in the example does have the jacobian in it, but for some reason my correct equation below is still not working. Are there other issues with polar coordinates I am not aware of?
0 = int2d(Th)(x^2*(a*dx(u)*dx(v) + dy(u)*dy(v) )) // bilinear form
+ int2d(Th)( (x^2*cos(y)*dx(u)*v + x*cos(y)*v*u) // linear form
+ any boundary conditions.