Beispiel - SAM - ETH Zürich

5
6 % plot tangent field
7 f i g u r e(’Name’,’LV field’); hold on;
8 % set up the sampling grid
9 N = 8; [X,Y] = meshgrid(0:1.5/N:1.5,0:1.5/N:1.5);
0 U = zeros( s i z e(X)); V = zeros( s i z e(Y));
1 % get velocity vectors
2 f o r i=0:N-1
3 f o r j=0:N-1
4 x = [1;fn(X(i+1,j+1),Y(i+1,j+1))];
5 x = 0.3*x/norm(x);
6 U(i+1,j+1) = x(1); V(i+1,j+1) = x(2);
7 end
8 end
9
0 % plot velocity vectors
1 quiver(X,Y,U,V,’b-’);
2 % set axes, labels, ...
3 set(gca,’fontsize’,14); axis([0 1.5 0 1.5]);
4 x l a b e l(’{\bf t}’); y l a b e l(’{\bf y}’);
5 % Create EPS output file
6 p r i n t -depsc2 ’riccatti2.eps’
Numerische
Mathemtik
R. Hiptmair
rev 35327,
25. April
2011
1.1
p. 17