无法解释的错误　　帮帮忙　等指教

下面的程序错误怎么也开不出来 什么原因 #define PAI 3.14159265 #include<math.h> main() {double N,Np,Si,u1,Qi,Ai,Rei,Pri,Nui,hi0,twi,Nwi,hi,twi1,Twi,Hi; double A=1,do=1,L=1,wh=1,phm=1,Si=1,cphm=1,th1=1,th20=1,di=1,u1=1,N1=900,rd1=1,Tmh=1; int n;

N=(int)(A/(pai*do*L))+1; Np=2;

Si=pai*di*di*N/(Np*4); u1=wh/(phm*Si); Qi=wh*cphm*(th1-th20); Ai=pai*di*L*N;

Rei=(di*u1*phm)/N1; Pri=(cphm*N1)/rd1; Nui=0.023*pow(Rei,0.8)*pow(Pri,0.33); hi0=Nui*rd1/di;

twi=Tmh-Qi/(hi0*Ai); Nwi=1000; hi=hi0*pow((N1/Nwi),0.14); twi1=Tmh-Qi/(hi*Ai); for(n=0;fabs(twi-twi1)>=0.5;n++) {twi=twi1; Nwi=niandu1(twi); hi=hi0*pow((N1/Nwi),0.14); twi1=Tmh-Qi/(hi*Ai); };

Twi=twi; Hi=hi; printf("Twi=%f,Hi=%f\n",Twi,Hi); printf("n=%d",n); }

----------------解决方案--------------------------------------------------------
头都看晕了,变量太多.
----------------解决方案--------------------------------------------------------
do=1；这个是最不应该的了！！！

niandu1(twi);什么意思？？？？

其他的自己看看！
----------------解决方案--------------------------------------------------------