10 REM
20 REM
30 REM             EXPERIMENT 3.A. HEAT CAPACITY RATIO OF GASES
40 REM                             ADIABATIC EXPANSION METHOD
50 REM
60 REM                        (LAB.14.WBASIC)
70 REM
80 REM
90 REM
100 REM
110 OPEN "e:\INPUT14.txt" FOR INPUT AS #3
120 OPEN "e:\OUTPUT" FOR OUTPUT AS #2
130 ON EOF(3) GOTO 1310
140 REM
150 REM                       "INPUT YOUR NAME"
160 INPUT #3, E$
170 REM                "INPUT THE NAME OF THE EXPERIMENT "
180 INPUT #3, a$
190 REM                  "INPUT THE NAME OF YOUR PARTNER "
200 INPUT #3, b$
210 REM                       "INPUT THE DATE "
220 LINE INPUT #3, C$
230 PRINT #2, "NAME: ", E$
240 PRINT #2, " "
250 PRINT #2, "EXPERIMENT: ", a$
260 PRINT #2, " "
270 PRINT #2, "PARTNER: ", b$
280 PRINT #2, " "
290 PRINT #2, "DATE: ", C$
300 FOR i = 1 TO 5
310 PRINT #2, " "
320 NEXT i
350 INPUT #3, bap, ebap
360 PRINT #2, "The barometric pressure during your experiment"
370 PRINT #2, "    ="; bap; "+/-"; ebap; "mm Hg"
380 PRINT #2, " "
390 INPUT #3, ddbp, eddbp
400 PRINT #2, "The density of dibutyl phthalate"
410 PRINT #2, "     ="; ddbp; "+/-"; eddbp; "g/cm3"
420 PRINT #2, " "
430 INPUT #3, dhg, edhg
440 PRINT #2, "The density of mercury"
450 PRINT #2, "    ="; dhg; "+/-"; edhg; "g/cm3"
460 PRINT #2, " "
470 INPUT #3, ep
480 PRINT #2, "The uncertainty in mm for the U-tube manometer"
490 PRINT #2, "containing the dibutyl phthalate is +/-"; ep; "mm Hg"
500 INPUT #3, ng
510 PRINT #2, " "
520 PRINT #2, ng; "gases were studied."
530 PRINT #2, " "
540 FOR i = 1 TO ng
550 INPUT #3, Q$(i)
560 PRINT #2, Q$(i)
570 PRINT #2, "========="
580 INPUT #3, m
590 PRINT #2, " ", "HEIGHTS OF DIBUTYL PHTHALATE"
600 PRINT #2, " ", "============================"
610 PRINT #2, TAB(9); "INITIAL EQUILIBRIUM"; TAB(41); "FINAL EQUILIBRIUM"
620 PRINT #2, TAB(10); "External  Internal         ";
630 PRINT #2, TAB(40); "External  Internal"
640 PRINT #2, TAB(10); " (open)   (closed)";
650 PRINT #2, TAB(40); " (open)   (closed)"
660 PRINT #2, " "
670 FOR j = 1 TO m
680 INPUT #3, poi(i, j), pci(i, j), pof(i, j), pcf(i, j)
690 ok(i, j) = poi(i, j) - pci(i, j)
700 ol(i, j) = pof(i, j) - pcf(i, j)
710 PRINT #2, "trial"; j; TAB(10); poi(i, j); TAB(20); pci(i, j); "=";
720 PRINT #2, ok(i, j); "mm DBP"; TAB(40); pof(i, j); TAB(50); pcf(i, j);
730 PRINT #2, "="; ol(i, j); "mm DBP"
740 NEXT j
750 PRINT #2, " "
760 NEXT i
770 PRINT #2, " "
780 PRINT #2, " "
790 lnp2 = LOG(bap)
800 elnp2 = ebap / bap
810 b = ddbp / dhg
820 eb = b * ((eddbp / ddbp) ^ 2 + (edhg / dhg) ^ 2) ^ .5
830 ec = ep * 2 ^ .5
840 FOR i = 1 TO ng
842 sum1(i) = 0!
844 esum1(i) = 0!
850 FOR j = 1 TO m
860 Ci(i, j) = poi(i, j) - pci(i, j)
870 ai(i, j) = b * Ci(i, j)
880 eai(i, j) = a * ((eb / b) ^ 2 + (ec / Ci(i, j)) ^ 2) ^ .5
890 p1(i, j) = bap + ai(i, j)
900 ep1(i, j) = (ebap ^ 2 + eai(i, j) ^ 2) ^ .5
910 lnp1(i, j) = LOG(p1(i, j))
920 elnp1(i, j) = ep1(i, j) / p1(i, j)
930 Cf(i, j) = pof(i, j) - pcf(i, j)
940 af(i, j) = b * Cf(i, j)
950 eaf(i, j) = a * ((eb / b) ^ 2 + (ec / Cf(i, j)) ^ 2) ^ .5
960 p3(i, j) = bap + af(i, j)
970 ep3(i, j) = (ebap ^ 2 + eaf(i, j) ^ 2) ^ .5
980 lnp3(i, j) = LOG(p3(i, j))
990 elnp3(i, j) = ep3(i, j) / p3(i, j)
1000 n(i, j) = lnp1(i, j) - lnp2
1010 en(i, j) = (elnp1(i, j) ^ 2 + elnp2 ^ 2) ^ .5
1020 d(i, j) = lnp1(i, j) - lnp3(i, j)
1030 ed(i, j) = (elnp1(i, j) ^ 2 + elnp3(i, j) ^ 2) ^ .5
1040 R(i, j) = n(i, j) / d(i, j)
1050 eR(i, j) = R(i, j) * ((en(i, j) / n(i, j)) ^ 2 + (ed(i, j) / d(i, j)) ^ 2) ^ .5
1060 sum1(i) = sum1(i) + R(i, j)
1070 esum1(i) = esum1(i) + eR(i, j) ^ 2
1080 NEXT j
1090 MR(i) = sum1(i) / m
1100 emR(i) = esum1(i) ^ .5 / m
1110 NEXT i
1120 FOR i = 1 TO ng
1130 PRINT #2, "Gas:"; Q$(i)
1140 PRINT #2, " ", "    Pressure, mm DBP", " Pressure, mm Hg"
1150 PRINT #2, " ", "    ================", " ==============="
1160 PRINT #2, " "
1170 FOR j = 1 TO m
1180 PRINT #2, "Trial ="; j
1190 PRINT #2, "Initial Equilibrium"; TAB(27); "=";
1200 PRINT #2, ok(i, j), " ", "    ="; p1(i, j)
1210 PRINT #2, "Final Equilibrium"; TAB(27); "=";
1220 PRINT #2, ol(i, j), " ", "    ="; p3(i, j)
1230 PRINT #2, " "
1240 PRINT #2, " ", " ", "Cp/Cv ="; R(i, j); "+/-"; eR(i, j)
1250 PRINT #2, " "
1260 NEXT j
1270 PRINT #2, TAB(29); "Average Cp/Cv ="; MR(i); "+/-"; emR(i)
1280 PRINT #2, " "
1290 PRINT #2, " "
1300 NEXT i
1310 CLOSE #2
1320 CLOSE #3
1330 END