MATLAB数值积分求立方体电势与可视化

空间中有一均匀带电的立方体，电荷密度 ρ = 1，边长为 a = 1。 使用数值积分求空间中各点的电场和电势

主程序（根据输入坐标值求该点电势，场强）：

% 开始计时
tic;    

% 接受用户输入的坐标
x_input = input('请输入 x 坐标：');
y_input = input('请输入 y 坐标：');
z_input = input('请输入 z 坐标：');

% 调用函数计算电场和电势
[E_field, V_potential] = caculate_E(x_input, y_input, z_input);


 % 停止计时并获取经过的时间
 time_elapsed = toc;

% 输出结果
disp(['Electric Field at (', num2str(x_input), ', ', num2str(y_input), ', ', num2str(z_input), '): ', num2str(E_field)]);
disp(['Electric Potential at (', num2str(x_input), ', ', num2str(y_input), ', ', num2str(z_input), '): ', num2str(V_potential)]);
disp(['Time Elapsed: ', num2str(time_elapsed), ' sec