Постройте линии тока и стагнации с помощью MATLAB

Question

Я хочу построить линии тока и точки застоя с учетом двух линейных вихрей противоположного направления, но равной величины.

Изображение того, как это выглядит, значения циркуляции гаммы:

Как выглядит мой код:

syms gamma r theta L

r1 = sqrt(r^2+L^2+2*r*L*cos(theta));
r2 = sqrt(r^2+L^2-2*r*L*cos(theta));
psi1 = (gamma/(2*pi))*log(r1);
psi2 = (-gamma/(2*pi))*log(r2);
psi = psi1+psi2;

vel_r = (1/r)*diff(psi,theta);
vel_theta = -diff(psi,r);
disp('Velocity Component r')
pretty(vel_r)
disp('Velocity Component theta')
pretty(vel_theta)