I am trying to plot following u_x and u_y functions, which are the x and y components of a velocity field.
The results should look like a vortex, as can be seen in the attached picture which I took from "wolframalpha.com" when I tried to compare my results.
However, instead such a simle 2D vortex I get complex velocity field which is also attached here.
I suspect it is due to wrong initialiazing the domain or wrong usage domain values in the "for loop", which is shown in the end of this message.
Could you please suggest me any possible solution.
I appreciate your time.