%
X=load('xfield.txt');
Y=load('yfield.txt');
VX=load('vxfield.txt');
VY=load('vyfield.txt');
%
sx=[];
sy=[];
%
k=1;
i=-2.48;
for j=-3.25:0.5:3.25
	sx(k)=i;
	sy(k)=j;
	k=k+1;
end
sx(k)=0.9;
sy(k)=0.0;
%
hold on;
streamline(X,Y,VX,VY,sx,sy);
t=0.0:0.001:2.0*pi;
plot(0.5*cos(t),0.5*sin(t),'-');
axis([-2.5,11.5,-3.5,3.5]);