function [x,y,u,v]=vtheta()

[x,y]=meshgrid([-750:1:750],[-750:1:750]);
X=complex(x,y);
        rhop=abs(X);
        ang=angle(X);
        Vtheta=(1-exp(-1.26*((rhop/5).^2)))./(2*pi*rhop);
        u=-Vtheta.*sin(ang);
        v=Vtheta.*cos(ang);
return
end
        