主题：我已经定义了实数的范围，可是实数的平方和却为零 求大大帮忙，谢谢
dingsheng1206
[专家分：0] 发布于 2011-06-06 16:17:00
integer k
real(8) b1,b0,e,w0,w,b,c,m,a
b0=234
if (k.eq.1) then
w0=6.0
w=22.2
b1=40.63
e=0.9047
else if (k.eq.2) then
w0=5.4
w=29.26
b1=44.17
e=0.9017
else if (k.eq.3) then
w0=4.7
w=36.33
b1=47.7
e=1.1522
end if
b=2*b1*e*w0*w*w+b0*(w0*w0-w*w)
c=b1*w*(w0*w0-w*w)-2*b0*e*w0*w
m=b**2+c**2
write(*,*) "input k please"
read(*,*)k
write(*,*) b,c,m
write(*,*) "input k please"
read(*,*) k
write(*,*) b,c,m
write(*,*) "input k please"
read(*,*) k
write(*,*) b,c,m
end
上面是我的程序，我已经定义了实数的范围，可是算出来结果却是平和为零，谢谢大大帮忙
回复列表 （共5个回复）
沙发
cgl_lgs [专家分：21040] 发布于 2011-06-06 21:40:00
程序是自上而下运行的，您这个逻辑明显不对啊：）
板凳
fortran2008 [专家分：750] 发布于 2011-06-07 11:02:00
integer k
real(8) b1,b0,e,w0,w,b,c,m,a
b0=234
10 write(*,*) "input k please"
read(*,*) k
if (k.eq.1) then
w0=6.0
w=22.2
b1=40.63
e=0.9047
else if (k.eq.2) then
w0=5.4
w=29.26
b1=44.17
e=0.9017
else if (k.eq.3) then
w0=4.7
w=36.33
b1=47.7
e=1.1522
end if
b=2*b1*e*w0*w*w+b0*(w0*w0-w*w)
c=b1*w*(w0*w0-w*w)-2*b0*e*w0*w
m=b**2+c**2
write(*,*) b,c,m
goto 10
end
4 楼
dingsheng1206 [专家分：0] 发布于 2011-06-08 08:16:00
[quote]程序是自上而下运行的，您这个逻辑明显不对啊：）[/quote]
谢谢你啊，
5 楼
dingsheng1206 [专家分：0] 发布于 2011-06-08 08:16:00
[quote]integer k
real(8) b1,b0,e,w0,w,b,c,m,a
b0=234
10 write(*,*) "input k please"
read(*,*) k
if (k.eq.1) then
w0=6.0
w=22.2
b1=40.63
e=0.9047
else if (k.eq.2) then
w0=5.4
w=29.26
b1=44.17
e=0.9017
else if (k.eq.3) then
w0=4.7
w=36.33
b1=47.7
e=1.1522
end if
b=2*b1*e*w0*w*w+b0*(w0*w0-w*w)
c=b1*w*(w0*w0-w*w)-2*b0*e*w0*w
m=b**2+c**2
write(*,*) b,c,m
goto 10
end[/quote]
谢谢你啊，
我来回复
您所在位置：