clear; clc;
Ux=0.25;
Qp=1;a=1;Q=Qp*a/pi();
[x,y]=meshgrid(-2:0.1:2,-2:0.1:2);
u=Ux+Q*(y.^2-x.^2)./(x.^2+y.^2).^2;
v=-2*Q*x.*y./(x.^2+y.^2).^2;
uplot=u./(u.^2+v.^2).^0.5;
vplot=v./(u.^2+v.^2).^0.5;
quiver(x,y,uplot,vplot,0.5)

ys=-2:.2:2;
xs=-2*ones(size(ys));
streamline(x,y,u,v,xs,ys)

ys=0.1;
xs=-((Q/Ux)-ys^2)^0.5;
streamline(x,y,u,v,xs,ys)
ys=-0.1;
streamline(x,y,u,v,xs,ys)

