请教一个程序问题，主要是循环部分!(急)

huahua2005

初步帮你看了一下，好像是这两句有问题：
sub=cos(-(Phi(k4,k3)-pi*r**2/lambda/lpp))*jj(k3)
subn=cos(-(Phi(k4,k3+1)-pi*rn**2/lambda/lpp))*jj(k3+1)
这两个值输不出来，你自己先改改看，明天再帮你看吧

sbn

caospace wrote:
这两个应该没有问题啊
sub=cos(-(Phi(k4,k3)-pi*r**2/lambda/lpp))*jj(k3)
subn=cos(-(Phi(k4,k3+1)-pi*rn**2/lambda/lpp))*jj(k3+1)
我能保证

好好检查前面把，计算到这里Phi已经是NaN了，自然要出错了

huahua2005

终于帮你找出来了，是你的粗心哦，下面这句中rz没赋值
rsegment=(radius0/nr*rz/r0)
导致下面几句中有分母为零，所以在运行中所有与下面三个变量有关的值，均为NaN，以后要细心哦
do 130 jr3=2, (nr)
d2qd2r=(Qn(jl3,(jr3+1))+Qn(jl3,(jr3-1))-2*Qn(jl3,jr3))/(rsegment**2)
dqdr=(Qn(jl3,(jr3+1))-Qn(jl3,(jr3-1)))/(2*(jr3-1)*rsegment**2)
Qn(jl3,jr3)=Qn(jl3,jr3)+kd*(t1/nt1)*(d2qd2r+dqdr)/cp
130 continue

sbn

huahua2005 wrote:
终于帮你找出来了，是你的粗心哦，下面这句中rz没赋值
rsegment=(radius0/nr*rz/r0)
导致下面几句中有分母为零，所以在运行中所有与下面三个变量有关的值，均为NaN，以后要细心哦
do 130 jr3=2, (nr)
d2qd2r=(Qn(jl3,(jr3+1))+Qn(jl3,(jr3-1))-2*Qn(jl3,jr3))/(rsegment**2)
dqdr=(Qn(jl3,(jr3+1))-Qn(jl3,(jr3-1)))/(2*(jr3-1)*rsegment**2)
Qn(jl3,jr3)=Qn(jl3,jr3)+kd*(t1/nt1)*(d2qd2r+dqdr)/cp
130 continue

呵呵，够细心的，说实话，他原来程序书写可读性比较不好

iwtry

看到这里相去以前自己写程序时,也是这样的错误不断,我想这些问题都不是什么太难的,自己学着静下心来慢慢看,自己的程序比较好看的,学用debug一步一步来,一定可以找出来一些端倪的.之所以找不出错误,是因为自己心浮,想一下子看出错误所在,或者闲费时间.hoho,