%《MATLAB及其在理工課程中的應用指南》第六章例6－4－2程序
% 線電荷的電場計算
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
clear all;
q = input('線電荷密度(庫侖/m)= (書上取1) ');
L = input('線電荷半長度 (m)= (書上取或50) ');
N = input('線電荷分段數 N= (書上取50或500)  '); 
Nr = input('離原點距離分段數Nr= (書上取50) ');
E0 = 8.85e-12;           % 真空電介質常數ε0
C0 = 1/4/pi/E0	;			% 歸并常數
L0 = linspace(-L,L,N+1);	% 將線電荷分N段
L1 = L0(1:N); L2=L0(2:N+1);	% 確定每個線段的起點和終點
Lm = (L1+L2)/2;dL= 2*L/N;	% 確定每個線段的中點坐標和長度,Lm是數組
R = linspace(0,10,Nr+1);	% 將R分N+1點
for k = 1:Nr+1				% 對R的N+1點循環計算
	Rk = sqrt(Lm.^2+R(k)^2);	% 測量點到電荷段的向徑
	Vk = C0*dL*q./Rk;		% 第k個電荷段在測量點處產生的電位
	V(k) = sum(Vk);			% 對各電荷段在測量點處產生的電位求和
end
[max(V),min(V)]				% 顯示最大最小電位
plot(R,V),grid				% 繪圖
set(gcf,'color','w')        % 置圖形背景色為白色