用matlab绘制单个点电荷电场分布，写出代码

以下是MATLAB代码实现单个点电荷电场分布：

% 定义常量
epsilon = 8.854187817e-12;    % 真空介电常数
q = 1;                       % 单个点电荷电量
k = 1 / (4 * pi * epsilon);  % 库仑常数

% 定义网格范围及步长
x = linspace(-1, 1, 101);
y = linspace(-1, 1, 101);
[X, Y] = meshgrid(x, y);

% 计算电场强度
r = sqrt(X.^2 + Y.^2);
Ex = k * q .* X ./ r.^3;
Ey = k * q .* Y ./ r.^3;

% 绘制电场场线及电势等高线
figure;
contour(X, Y, sqrt(Ex.^2 + Ey.^2), 25);
hold on;
h = streamslice(X, Y, Ex, Ey, 2);
set(h, 'Color', 'r', 'LineWidth', 1.5);
title('电场分布');
xlabel('x');
ylabel('y');

以上代码会绘制出单个点电荷的电场分布图。当然，具体的图形效果还可以根据需要进行调整。