这是我的计算过程，我想给图片里的 方程加的边界条件是当x^2+y^2=0.03时，D1=1，因为自变量是x,y,z，请问这种边界条件怎么加呢？先蟹蟹大家啦
eqn6 = (-8.30996*10^6 - 1/(x^2 + y^2)) D1[x, y,
z] + (0. + 2.371*10^7 I)
\!\(\*SuperscriptBox[\(D1\),
TagBox[
RowBox[{"(",
RowBox[{"0", ",", "0", ",", "1"}], ")"}],
Derivative],
MultilineFunction->None]\)[x, y, z]
DSolve[eqn6 == 0, D1[x, y, z], {x, y, z}]

这是我后面写的，主要问题是它总是提示应该使用方程或者方程列表，是边界条件哪里有问题吗

……DSolve还不支持区域语法。换成极坐标可能有点希望。

用 NDSolveValue，支持 Element[ {x,y} ,Disk[]] 的语法。