麻烦高手改下这个MATLAB,

　　fcu0=26.8;
　　c=35;
　　d=12;
　　ti=53.0;
　　deltacr=0.012*c/d+0.00084*fcu0+0.018;
　　Kcl=1.6;
　　lambda0=7.53*Kcl*4.2*(0.75+0.0125*17)*(0.75-0.45)^(2/3)*c.^(-0.675).*fcu0.^(-1.8);
　　tcr=ti+deltacr./lambda0;
　　lambda1=(4.5-340*lambda0).*lambda0;
　　rhos=zeros(1,100);%%定义rhos的大小
　　delta0=0;%%初始化delta0
　　for t=1:100%%从1开始循环,默认步长为1
　　if tti&t0.25&q00.3&q0>0.44)
　　kb=0.92-0.63*q0;
　　elseif find(o10.25&q0(i)