鉛の絶対エントロピーの計算

その温度範囲で位相変化を起こさない物質を温度T1からT2まで加熱したときの
エントロピー変化は次の式で与えられる。
10 ' Program entropy
20 ' 実験データに対してsimpson則を適用する例
30 DIM X(100),Y(100)
40 DATA 0,0,5,.061,10,.28,15,.4666,20,.54,25,.564,30,.55,50,.428
45 DATA 70,.333,100,.245,150,.169,200,.129,250,.105,298,.089
50 N = 14 ' データ数
55 PRINT"入力データをそのまま表示したもの:"
60 FOR I=1 TO N
65 READ X(I), Y(I)
70 PRINT X(I), Y(I)
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 '
140 S = S+(X(I+1)-X(I))*(Y(I+1)-Y(I))*0.6667 ' この部分は放物線1,4の
150 S = S+(X(I+2)-X(I+1))*(Y(I+1)-Y(I+2))*0.6667 ' 下の面積を計算する
160 NEXT I
170 PRINT : PRINT "エントロピー (変化) は :";S
180 END

entropy.bas

問題 298Kにおける鉛のエントロピーを計算しなさい。与えられたデータは温度Tを独立変数とみなしたときのCp/T-Tの実験値からなっている。
ここで、SはエントロピーでCpは定圧熱容量である。
T Cp/T
0 0
5 .061
10 .28
15 .4666
20 .54
25 .564
30 .55
50 .428
70 .333
100 .245
150 .169
200 .129
250 .105
298 .089

2002 Copyright 木原 寛