Please make sure that your code reproduces the error message when copy-pasting it into a new file.
You are missing several key definitions, such as a mesh (can be a built in one such as a UnitSquareMesh
) and the variables theta
and r
. Note that you can remove the line solve(a==L, w, bcs)
, as the code is only dependent on the split of w
, not the actual values.
For more information about how to easily get your questions answered, see: Read before posting: How do I get my question answered?