（未完成）matlab--矩阵运算技巧

//2014年7月28日入“未完成”

[x,y]=meshgrid(10:10:1000);
s=(x.^2+y.^2).^0.5;
jiaodu=atan(y./x);  %atan
for i =-n:n  %-n:n
    jiaodu2=i/(n+1)*pi/4
    newy=s*sin(jiaodu2);
    newx=s*cos(jiaodu2);
    m=airflow(newx,newy,0);
    if(jiaodu+jiaodu2 >= 3/4*pi)%t=t+0.1087.*(jiaodu+jiaodu2 >= 3/4*pi).*m;
        t=t+0.1087*airflow(newx,newy,0);
    elseif(jiaodu+jiaodu2 > 1/4*pi)%t=t+0.5435.*(jiaodu+jiaodu2 < 3/4*pi && jiaodu+jiaodu2 > 1/4*pi).*m;
        t=t+0.5435*airflow(newx,newy,0);
    elseif(jiaodu+jiaodu2 <= 1/4*pi)%t=t+0.3478.*(jiaodu+jiaodu2 <= 1/4*pi).*m;