金属の絶対エントロピーの計算

10 'PROGRAM entropy2
20
' 実験データに対してsimpson則を適用する例
30 DIM X(100),Y(100)
40 DATA 0,0, 5,.24, 10,.64, 15,1.36, 20,2.31, 25,3.14, 30,4.48, 50, 9 .64
45 DATA 70,15.7, 100,20.2,150, 22!, 200,23.4, 250,24.3, 298,25.5
50 N = 14
55 PRINT "入力データをそのまま表示したもの:"
60 FOR I=1 TO N
65 READ T, CP
70 PRINT T, CP
75 X(I)=T : IF T<> 0 THEN Y(I)=CP/T ELSE Y(I)=0
80 NEXT I
90 FOR I=0 TO N-2 STEP 2
'この部分は矩形
100 S = S + (X(I+1)-X(I))*(Y(I))
' ブロックの面積を
110 S = S + (X(I+2)-X(I+1))*Y(I+2)
' 計算する
120 NEXT I
130 FOR I=0 TO N-2 STEP 2
' この部分は放物線1,4の
140 S = S + (X(I+1)-X(I))*(Y(I+1)-Y(I))*.6667
'下の面積を計算する
150 S = S + (X(I+2)-X(I+1))*(Y(I+1)-Y(I+2))*.6667
'演習 2-1参照
160 NEXT I
170 PRINT : PRINT "エントロピー (変化) は :";S
180 END

entropy2.bas

問題 ある光沢を有する金属の熱容量が次に示すような温度の関数であるとする。   298Kにおけるエントロピーを計算しなさい。標準エントロピーの表を調べ、この金属が何か推定しなさい。

T Cp/T
0 0
5 .24
10 .64
15 1.36
20 2.31
25 3.14
30 4.48
50 9.64
70 15.7
100 20.2
150 22
200 23.4
250 24.3
298 25.5
2002 Copyright 木原 寛