这是通过话梯度和等高线画电场线和等势线的程序 
[x,y]=meshgrid([-1:.05:1]); 
clf 
U=1./sqrt(x.^2+y.^2); 
[dx,dy]=gradient(U,.05,.05); 
contour(x,y,U) 
hold on 
 
axis([-1 1 -1 1]); 
quiver(x,y,dx,dy,50,'y') 
colormap hsv 
hold off 
但是画不出电场线, 
把U=1./sqrt(x.^2+y.^2);(line 3)的方程改成U=x.^2+y.^2后梯度的向量就画出来了 
为什么哦,是不是梯度场的大小原因? 
等待高手......
		 
		
		
		
		
		
		
		
	 |