포트란으로 아래 식에서 ir 값을 계속 바꿔서 데이타 얻어내야 하는데 답이 안 나옵니다
원래는 플랑크 식에서 나온 건데 이렇게 수정해서 풀라고 한거
상수값이나 이런거는 아래 코드식에서 쓸 value로 바꿔놨어염
제가 짠 코드:
program hw1
implicit none
integer::i
real::x,y,z,r,result1, bm, dd1, dd2, ir
CHARACTER(10) :: filename="hw1.dat"
OPEN(UNIT=2,FILE=filename,STATUS="NEW",ACTION="WRITE")
x=3.742E+08
y=1.4388E+04
z=27
r=10
PRINT*, x,y,z,r
Do i = 1, 200, 1
ir = i*0.01
bm = y/(ir*z)
dd1 = ir**5
dd2 = exp(bm)
result1 = x/(dd1*(dd2-1))
WRITE(2,*), i, bm, dd1, dd2, result1
end do
CLOSE(2)
End program hw1
이렇게 해봤는데 답이
1 53288.89 9.9999987E-11 Inf 0.000000
2 26644.45 3.1999996E-09 Inf 0.000000
3 17762.96 2.4299998E-08 Inf 0.000000
4 13322.22 1.0239999E-07 Inf 0.000000
5 10657.78 3.1249991E-07 Inf 0.000000
6 8881.481 7.7759995E-07 Inf 0.000000
7 7612.698 1.6807001E-06 Inf 0.000000
8 6661.112 3.2767996E-06 Inf 0.000000
9 5920.988 5.9048984E-06 Inf 0.000000
10 5328.889 9.9999970E-06 Inf 0.000000
11 4844.444 1.6105099E-05 Inf 0.000000
12 4440.741 2.4883198E-05 Inf 0.000000
13 4099.146 3.7129295E-05 Inf 0.000000
14 3806.349 5.3782402E-05 Inf 0.000000
15 3552.593 7.5937482E-05 Inf 0.000000
16 3330.556 1.0485759E-04 Inf 0.000000
17 3134.640 1.4198571E-04 Inf 0.000000
18 2960.494 1.8895675E-04 Inf 0.000000
19 2804.678 2.4760992E-04 Inf 0.000000
20 2664.445 3.1999990E-04 Inf 0.000000
이것도 다 나온게 아니라 38번째 밖에 또 안 나왔어여
혹시 뭐가 잘못됐는지 고수님들께 물어봅니다...
굽신
댓글 0