带电小圆柱的静电场

%  高等电磁理论 第一次作业%  %  设带电小圆柱几何中心为坐标原点，长轴为z轴% % 圆柱参数clear;clf;clc;cR=0.5;    % 半径cL=4;    % 长度cRho=1;  % 电荷密度cK=1;    % 系数% 坐标轴[crdX,crdY,crdZ]=meshgrid(-5:1:5);% 求电势ticphi=(cRho/cK)*arrayfun(@(cX,cY,cZ)quadl(@(x)arrayfun(...    @(zz)quad2d(@(y,x)1./sqrt((cX-x).^2.+(cY-y).^2.+(cZ-zz).^2),...    -cR, cR, @(y)-sqrt(cR^2-y.^2), @(y)sqrt(cR^2-y.^2)), x), -cL/2, cL/2), crdX, crdY, crdZ);toctic[u,v,w]=gradient(-phi); % 对电势求梯度得到电场