请问高手下面的程序怎么运行啊?

yamalu

偶是刚接触matlab,弱的不知道为什么不能运行啊!!!请达人指点!!!
k=1:3000
D=0.005;
A=3.1415926/4*D*D;
%h=1;                                                           
%t=1e-5;
rouN=1.293;
miu=(4.8e-8)*20+1.71e-5;
G(1)=8.3336;
Deta_G(1)=0;
u(1)=2.5745;
rou(1)=3.553;
Re(1)=2532.7;
lamta(1)=0.044603;
  
  for j=2:3000
      %if  (lamta(j-1)==0 | u(j-1)==0)   
          %Deta_G(j)=0;
      %else
          Deta_G(j)=-A*D/(lamta(j-1)*u(j-1))*DP_t(j)*1e-3*(exp(-lamta(j-1)*u(j-1)/D*(1+j*1e-3))-1)/rouN*60000;
      
      G(j)=G(j-1)+Deta_G(j-1);
      u(j)=G(j)*rouN/(rou(j)*A*60000);
      Re(j)=rou(j)*u(j)*D/miu;
      
      %if  Re(j)<=0  
          %lamta(j)=0;
       %else if    (Re(j)>0 & Re(j)<=2300)   
              %lamta(j)=64/Re(j);
       %else  
              lamta(j)=0.3164*Re(j)^(-0.25);
       %end
           
   end
   plot(1+k*1e-3,G(k));

eight

原帖由 yamalu 于 2007-5-22 15:37 发表

                               
登录/注册后可看大图

偶是刚接触matlab,弱的不知道为什么不能运行啊!!!请达人指点!!!
k=1:3000
D=0.005;
A=3.1415926/4*D*D;
%h=1;                                                           
%t=1e-5;
rouN=1.293;
miu ...

请先看看 置顶帖：聚宝盆，然后重新整理你的问题，特别是给出错误信息，不要指望别人帮你调试