cheng2.asv

高 斯 迭 代

q1=1;
q2=-3;
a=0.05;
xx0=0.1;
yy0=0.08;
k=9e9;
ezplot('-3*9/sqrt((xx-1)^2+yy^2)+1*9/sqrt((xx+1)^2+yy^2)-0',[-1,1,-0.8,0.8]);
hold on;
ezplot('-3*9/sqrt((xx-1)^2+yy^2)+1*9/sqrt((xx+1)^2+yy^2)-0',[-1,1,-0.8,0.8]);
hold on;
%y=linspace(-yy0,yy0,20);
%[X,Y]=meshgrid(x,y);
%r2=sqrt((X + a).^2 + Y.^2);
%r1=sqrt((X - a).^2 + Y.^2);
%U=k*q1./r1+k*q2./r2;
%u0=8e11;
%u=linspace(u0,-u0,11);
%contour(X,Y,U,u,'r-');
hold on;
%[Ex,Ey]=gradient(-U);
%E=sqrt(Ex.^2+Ey.^2);
%Ex=Ex./E;
%Ey=Ey./E;
%quiver(X,Y,Ex,Ey);
C2=0;
for i=1:30   
ezplot('-3*(xx-1)/sqrt((xx-1)^2+yy^2)+1*(xx+1)/sqrt((xx+1)^2+yy^2)-C2',[-1,1,-0.8,0.8]);
hold on;
C=C+0.1;
end
hold off;