%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第六章例6－4－2程序
% 電位分布及電場計算
% 西安電子科技大學出版社出版  陳懷琛編著 2007年5月
%
fprintf('輸入電位分布方程 V(x,y) \n');
fprintf('例如:  log(x.^2 + y.^2) \n');
V = input(': ','s');   % 讀入字符串 V(x,y)
NGrid = 20;            % 繪圖的網格線數
xMax = 5;              % 繪圖區從 x= -xMax 到 x= xMax
xPlot = linspace( -xMax, xMax, NGrid); % 繪圖取的x值
%yPlot = xPlot; %linspace( -yMax, yMax, NGrid); % 
[x,y]=meshgrid(xPlot);
VPlot=eval(V);
[ExPlot, EyPlot] = gradient(-VPlot);		%電場是電位的負梯度
clf; subplot(1,2,1),meshc(VPlot);  % 畫含等高線的三維曲面
xlabel('x');  ylabel('y'); zlabel('電位');
% 規定等高線圖的范圍及比例
subplot(1,2,2), axis([-xMax xMax -xMax xMax]);    % 建立第二子圖
cs = contour(x,y,VPlot); % 畫等高線
clabel(cs); hold on;          % 在等高線圖上加上編號
%@ 在等高線圖上加上電場方向
quiver(x,y,ExPlot,EyPlot);  % 畫電場 E 的箭頭圖
xlabel('x');  ylabel('y');hold off;