当rb时,体密度和面密度=0.
aP=D-εE.
D={4/3*pi*（r^3-a^3)*p}/4*pi*r^2,注意要有方向的,D,E,P都是向量.
内表面p=0；
外表面p=..就是把b代入
体密度,p=p的负散度.书上有公式的