% 定义常数
epsilon_0 = 8.854e-12; % 真空介电常数
% 定义参数
L = 2; % 线电荷的长度(单位:米)
lambda = 1e-9; % 线电荷的电荷密度(单位:库仑/米)
% 定义计算区域
x = linspace(-5, 5, 100); % x坐标范围
y = linspace(-5, 5, 100); % y坐标范围
[XX, YY] = meshgrid(x, y); % 创建坐标网格
% 计算电场和电势
E = zeros(size(XX)); % 存储电场值
V = zeros(size(XX)); % 存储电势值
for i = 1:numel(x)
for j = 1:numel(y)
% 计算点(x(i), y(j))处的电场和电势
r = sqrt(XX(i, j)^2 + YY(i, j)^2); % 点(x(i), y(j))到原点的距离
theta = atan2(YY(i, j), XX(i, j)); % 点(x(i), y(j))与x轴正方向的夹角
% 计算电势
V(i, j) = (lambda / (4 * pi * epsilon_0)) * log(r / sqrt(r^2 + L^2));
% 计