у активностей міді та нікелю
20 LET r = 8.3144: REM універсальна газова константа
30 LET t = 298: REM температура 25 С = 298 К
40 LET q1 = -16369.02: REM Q12 (1) при 25 С
50 LET q2 = 24136.67: REM Q12 (2) при 25 С
60 OPEN "Activ.txt" FOR OUTPUT AS # 1
70 PRINT # 1, "X1", "x2", "a1", "a2"
80 FOR x = .01 TO .99 STEP .01
90 REM x = x1 - мольна частка міді
100 LET y = 1 - x: REM y = X2 - мольна частка нікелю
110 LET z1 = r * T * LOG (x) + (y ^ 2) * (2 * x * q1 + (1 - 2 * x) * q2)
120 LET z2 = r * T * LOG (y) + (x ^ 2) * ((1 - 2 * y) * q1 + 2 * y * q2)
130 REM z1 = RTln a1 і z2 = RTln a2
140 LET a1 = EXP (z1/(r * T)): REM активність міді
150 LET a2 = EXP (z2/(r * T)): REM активність нікелю
160 PRINT # 1, x, y, a1, a2
170 NEXT x
180 CLOSE # 1
190 END
Додаток В
Текст програми для розрахунку активностей в потрійній системі наведено нижче:
10 CLS
20 LET r = 8.3144: REM r - Універсальна газова постійна
30 LET t = 298: REM t - температура, 25 С = 298 К
40 LET rt = r * T: REM RT
50 REM 1 - Cu, 2 - Ni, 3 - Be
60 LET x1 = .868: LET x2 = .005: LET x3 = .127: REM молярний частки в бронзі БрБ2
70 DIM q1 (3, 3): DIM q2 (3, 3)
80 REM елемент масиву q1 (i, j) - це Qij (1) при температурі t
90 REM елемент масиву q2 (i, j) - це Qij (2) при температурі t
100 LET q1 (1, 2) = 76.812 * t - 39259
110 LET q2 (1, 2) = -51.085 * t + 39360
120 LET q1 (1, 3) = -10.55 * t + 24080
130 LET q2 (1, 3) = q1 (1, 3)
140 LET q1 (2, 3) = 18.25 * t - 34260
150 LET q2 (2, 3) = q1 (2, 3)
160 REM обчислення z1 = RTln a1
170 LET z1 = rt * LOG (x1) + 2 * x1 * x2 * (1 - x1) * q1 (1, 2) + (x2 ^ 2) * (1 - 2 * x1) * q2 (1, 2)
180 LET z1 = z1 + 2 * x1 * x3 * (1 - x1) * q1 (1, 3) + (x3 ^ 2) * (1 - 2 * x1) * q2 (1, 3)
190 LET z1 = z1 - 2 * x2 * x2 * x3 * q1 (2, 3) - 2 * x2 * x3 * x3 * q2 (2, 3)
200 REM обчислення z2 = RTln a2
210 LET z2 = rt * LOG (x2) + 2 * x1 * x2 * (1 - x2) * q2 (1, 2) + (x1 ^ 2) * (1 - 2 * x2) * q1 (1, 2)
220 LET z2 = z2 + 2 * x2 * x3 * (1 - x2) * q1 (2, 3) + (x3 ^ 2) * (1 - 2 * x2) * q2 (2, 3)
230 LET z2 = z2 - 2 * x1 * x1 * x3 * q1 (1, 3) - 2 * x1 * x3 * x3 * q2 (1, 3)
240 REM обчислення z3 = RTln a3
250 LET z3 = rt * LOG (x3) + 2 * x1 * x3 * (1 - x3) * q2 (1, 3) + (x1 ^ 2) * (1 - 2 * x3) * q1 (1, 3)
260 LET z3 = z3 + 2 * x2 * x3 * (1 - x3) * q2 (2, 3) + (x2 ^ 2) * (1 - 2 * x3) * q1 (2, 3)
270 LET z3 = z3 - 2 * x1 * x1 * x2 * q1 (1, 2) - 2 * x1 * x2 * x2 * q2 (1, 2)
280 REM обчислення активностей
290 LET a1 = EXP (z1/rt)
300 LET a2 = EXP (z2/rt)
310 LET a3 = EXP (z3/rt)
320 PRINT a1, a2, a3
330 END...
