[求助]已解决
程序如下:function q=Ex0001()
% Ex0001 收缩喷管的流量压力关系图
k=1.3;
zeta1=0:0.01:1;
if (zeta1>0.546)&(zeta1<1)
q=sqrt((2/(k-1))*((k+1)/2).^((k+1)/(k-1))*(zeta1.^(2/k)-zeta1.^((k+1)/k)));
elseif (zeta1>0)&(zeta1<0.546)
q=1;
else zeta1==1
q=0;
end
axis([0 2 0 2]);
plot(zeta1,q)
[此贴子已经被作者于2006-5-26 21:02:17编辑过]