Теоретическое поле линий тока и линий равного потенциала при бесциркуляционном обтекании цилиндра можно построить в пакете Maple с помощью команды
plots[display](
seq(plot([r*phi/(r+1/r),r*sqrt(1-(phi/(r+1/r))^2),r=1..6],numpoints=1000,colour=black),phi=0..5,0.25),
seq(plot([r*sqrt(1-(psi/(r-1/r))^2),r*psi/(r-1/r),r=1..6],numpoints=1000,colour=black),psi=0..5,0.25),
plot([cos(t),sin(t),t=0..Pi/2],thickness=2,colour=black),
view=[0..4,0..4],labels=["x/R","y/R"],scaling=constrained
);
В результате получается ортогональная сетка следующего вида.
Изображение в djvu-формате доступно здесь.