如何简化这个程序, 谢谢

[slmyzcc]

%%%%% PROFILS DE TEMPERATURE DANS LES CANAUX %%%%%%
Tairv=zeros(NP0,Nbseg);
for i=1:NP0
Tdis=(T1(i)+T2(i))/2;
Tairc=Te(i);
for k=1:Nbseg
Tairv(i,k)=Tairc-273.15;
xx = Llc * k / Nbseg;
cpa=Cp_air(Tairc);
lama=lam_air(Tairc);
Rho=Rho_air(Tairc);
Nyu=Viscine_air(Tairc);
Myu=Nyu*Rho;
Myum=Viscine_air(Tdis)*Rho_air(Tdis);
Secx=PenteSec*xx+CsteSec;
Perx=PentePer*xx+CstePer;
Surx=0.5*PentePer* (2 * k - 1) * (Llc / Nbseg)^2+CstePer*Llc/Nbseg;
Dhx=4*Secx/Perx;
Ux=Qm(i)/(Rho*Secx);
Rex=Ux*Dhx/Nyu;
Nux=Va*Rex^Vb*(Myu/Myum)^Vc*(1+(xx/Llc)^Vd);
hx=lama*Nux/Dhx;
phix=hx*Surx*(Tdis-Tairc);
Tairc=Tairc+phix/(Qm(i)*cpa);
end
end