10  REM ****  HP TIME-SHARED BASIC PROGRAM LIBRARY  *******************
20  REM
30  REM             KINET:  KINETIC REACTION
40  REM
50  REM             36309  (A833)  REV A -- 10/71
60  REM
70  REM ****  CONTRIBUTED PROGRAM  ************************************
80  REM KINET - COPYRIGHT 1971, POLYTECHNIC INSTITUTE OF BROOKLYN
90  REM  KINET  HOWARD SHANNON  HARBORFIELDS HS  8/15/68 
100  REM      REVISED 7/28/69  PIB J. MARCHISOTTO 
110  REM  REVISED BY C.LOSIK  7-28-70 
120  REM  F,K ARE DEFINED BELOW; L IS THE TIME INCREMENT; 
130  REM  G IS THE DECAY CONSTANT; D IS THE FRACTION OF 'P' 
140  REM  IF L IS TOO SMALL, ADJUST LINE 530. 
150  REM  IF YOU WANT A VARIABLE TIME INTERVAL, MAKE 530 LET L=.69/G
160  REM     THIS PROGRAM STUDIES THE KINETICS OF A SINGLE SPECIES
170  REM (A) GOING TO A SINGLE SPECIES (P),AND APPROACHES EQUILIBRIUM 
180  REM WITH (P).
190  REM          IE. ISOMER EQUILIBRIUM
200  REM      IF THE EQUILIBRIUM CONSTANT IS VERY LARGE (K>10,000), 
210  REM IT CAN BE ASSUMED THAT ALL OF THE REACTANT GOES TO PRODUCT.  
220  REM THIS PROGRAM CAN THEN BE USED FOR RADIOACTIVE DECAY. 
230  REM      AN INPUT OF THE FORWARD RATE CONSTANT AND THE EQUILIBRIUM 
240  REM CONSTANT WILL GIVE A PRINTOUT OF THE CONCENTRATION OF (A)
250  REM AND (P) COMPARED TO THE INITIAL CONCENTRATION OF (A) AT 10 
260  REM EQUAL TIME INTERVALS AS IT APPROACHES 99.9PERCENT TO EQUILIBRIUM.
270  PRINT "  FOR THE EQUILIBRIUM PROBLEMS YOU ARE ABOUT TO DO, THE "
280  PRINT "  DATA MAY BE PRESENTED IN THE FOLLOWING MANNER:  "
290  PRINT "     (INDICATE YOUR CHOICE BY NUMBER)"
300  PRINT 
310  PRINT "            CHOICE 1 = TABLE OF DATA"
320  PRINT "            CHOICE 2 = GRAPH OF DATA"
330  PRINT "            CHOICE 3 = TABLE AND GRAPH OF DATA"
340  PRINT "            CHOICE 4 = END PROGRAM"
350  PRINT 
360  PRINT " WHAT IS YOUR CHOICE";
370  INPUT Q
380  IF Q>4 THEN 870
390  IF Q<1 THEN 870
400  IF Q=4 THEN 890
410  PRINT 
420  PRINT " LET F = THE FORWARD RATE CONSTANT"
430  PRINT " LET K = THE EQUILIBRIUM CONSTANT FOR THE REACTION A =P"
440  PRINT " TYPE IN THE CONSTANTS F AND K IN THAT ORDER."
450  INPUT F,K
460  PRINT 
470  PRINT 
480  PRINT "               ******************************"
490  PRINT 
500  LET H=K/(K+1)
510  LET G=F/H
520  LET L=.05
530  LET T=-L
540  IF Q=2 THEN 680
550  PRINT " LET A1 = ORIGINAL CONCENTRATION OF A"
560  PRINT " LET A = PERCENT CONCENTRATION OF A  (A/A1*100)"
570  PRINT " LET P = PERCENT CONCENTRATION OF P  (P/A1*100)"
580  PRINT 
590  PRINT " TIME ","  A  ","  P  "
600  PRINT " ---- "," --- "," --- "
610  FOR I=0 TO 10
620  LET T=T+L
630  LET D=H*(1-EXP(-G*T))
640  PRINT T,(1-D)*100,D*100
650  NEXT I
660  LET T=-L
670  IF Q=1 THEN 350
680  PRINT 
690  PRINT "                 PERCENT CONCENTRATION OF A(*) AND P(+)"
700  PRINT 
710  PRINT " ","0        25        50        75        100"
720  PRINT "  TIME","I---------I---------I---------I---------I"
730  FOR I=0 TO 10
740  LET T=T+L
750  LET D=H*(1-EXP(-G*T))
760  LET X=1-D
770  PRINT I*L,"I";
780  IF INT(40*D+.5)=INT(40*X+.5) THEN 820
790  IF INT(40*D+.5)>INT(40*X+.5) THEN 840
800  PRINT TAB(40*D+14.5);"+";TAB(40*X+14.5);"*"
810  GOTO 850
820  PRINT TAB(40*D+14.5);"."
830  GOTO 850
840  PRINT TAB(40*X+14.5);"*";TAB(40*D+14.5);"+"
850  NEXT I
860  GOTO 350
870  PRINT " YOUR CHOICE MUST BE A NUMBER BETWEEN 1 AND 4, TRY AGAIN."
880  GOTO 350
890  END 
