求解下面的編程 求助紅色區(qū)域 是

sun7854

d1=0.001,d2=1.5;

while(1)

{

d3=(d1+d2)/2.0;

Re=1.534*Q*del*1000/(d3*u);

R1=0.001;

R2=1.0;

while(1)

{

R3=(R1+R2)/2.0;if(fr(k,d3,Re,R3)>1e-6)R1=R3; if(fr(k,d3,Re,R3)<-1e-6)R2=R3;if(fabs(fr(k,d3,Re,R3))<1e-6)break;

}

R=R3; if(fd(R,z,d3)>1e-6)d1=d3;

if(fd(R,z,d3)<-1e-6) d2=d3;if(fabs(d3-d1)<1e-6&&fabs(d3-d2)<1e-6)break;

}

d=d3; printf("摩阻系數(shù)為R=%e,管徑為d=%em\n\n\n",R,d); printf("雷諾數(shù)為Re=%e\n\n\n",Re);

}

sun7854

截取了自己不懂得地方  ，發(fā)布出來全部的編程  點我QQ   謝謝前輩大神們

hqbtc

兩個while循環(huán)是死循環(huán)
因為條件為1永遠(yuǎn)為真
if語句的break條件滿足時，跳出循環(huán)
具體結(jié)果自己分析

sun7854

hqbtc 發(fā)表于 2017-3-13 17:16
兩個while循環(huán)是死循環(huán)
因為條件為1永遠(yuǎn)為真
if語句的break條件滿足時，跳出循環(huán)

程序沒問題，調(diào)用子函數(shù)算的fr和fd。就是上面的管徑d1和d2的賦值0.001和1.5，這是輸氣管道管徑的編程。發(fā)不了全部的只發(fā)了我不懂的這個