% Rotasjonsfelt

[x,y]=meshgrid(-30:2:30,-20:2:20);      % definerer omraadet xy 

omega = 1;                              % setter konstanten omega

u=-omega*y;                             % x-komp av vektor
v=omega*x;                              % y-komp av vektor

quiver(x,y,u,v);                    % quiver plotter vektorer
axis equal;
hold on

% konturlinjer/str鴐linjer
[x,y]=meshgrid(-30:0.1:30, -20:0.1:20); % definerer omraadet xy
rsq=((x.*x)+(y.*y));                    % regner ut r^2 = x*x+y*y
psi=0.5*omega*rsq;                      % str鴐funksjonen/feltfunksjonen

t = [10 50 100 200 300 400];            % velger hvilke konturlinjer som skal plottes
contour(x,y,psi,t);