n=6%定义轮数
r=(0:n)'/n;%定义轮的半径
theta=pi*(-n:n)/n;%定义轮的扇区角
X=r*cos(theta);
Y=r*sin(theta);%定义网格顶点
C=r*cos(2*theta);%定义色图
pcolor(X,Y,C)%绘制伪彩图axis equal tight

[x y z v]=flow;%打开水流数据
h=contourslice(x,y,z,v,[1:9],[],[0],linspace(-8,2,10));%切片轮廓线
view([-12 30])