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

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

[x,y]=meshgrid(10:10:1000);s=(x.^2+y.^2).^0.5;jiaodu=atan(y./x);  %atanfor 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;