[x,y] = meshgrid(-5:0.1:5,-5:0.1:5);
r = sqrt(x.^2+y.^2);
theta = atan2(y, x);
gamma = 2;
B = 2;
u = 1;
psi = u.*r.*sin(theta)-B./r.*sin(theta)-gamma/(2*pi).*log(r);
contour(x,y,psi)
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!