1  REM 73.05.19, JHP
70  DIM V[6],W[29],I[6],Z[1]
80  DIM A[2,20],S[3]
85  DIM A$[72]
90  DIM T[9]
140  REM-DATA CONTAINS INITIAL VALUES FOR ALL VARIABLES AND 
150  REM-PRESET VALUES FOR OPERATION OF THE STABLIZING SYSTEM 
155  MAT  READ A
160  DATA 20,30,40,50,60,70,80,90,100,110
170  DATA 120,130,140,150,160,170,180,190,200,250
180  DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
185  MAT  READ W
190  DATA 20,100,100,100,-2,8,5000,2.2,1.3,40
195  DATA 40,40,39.2,38.2,36,29,19,11,8,6
196  DATA 5,4,4,4,4,4,4,4,4
197  PRINT "MCMAN WAS DEVELOPED AND IS CURRENTLY USED IN THE MEDICAL"
198  PRINT "CURRICULUM AT MCMASTER UNIVERSITY, HEALTH SCIENCES CENTRE,"
199  PRINT "HAMILTON, ONTARIO, CANADA."
200  PRINT 
201  PRINT "NEITHER HEWLETT-PACKARD NOR MCMASTER UNIVERSITY ASSUME"
202  PRINT "ANY RESPONSIBILITY FOR THE OPERATION AND CORRECTNESS OF THIS MODEL"
203  PRINT 
205  U1=.1
210  MAT I=ZER
220  REM-INITIALIZE ALL VALUES AND INDICES
230  J5=-15
240  J4=0
250  J7=0
260  C5=0
270  J3=-60
280  K2=0
290  K4=0
300  H1=72
310  C1=163
320  D1=16
330  T[9]=W[9]
335  FOR I1=1 TO 3
336  S[I1]=W[I1]
337  NEXT I1
360  REM******RESTORE PREVIOUS VALUES 
380  FOR I1=1 TO 9
381  T[I1]=W[I1]
382  NEXT I1
390  J1=1
400  FOR I1=10 TO 29
410  A[2,J1]=W[I1]
420  J1=J1+1
430  NEXT I1
440  REM-******TYPEWRITER INTERACTION PART 1
450  PRINT  USING 3000
455  GOSUB 3500
490  GOTO J9 OF 500,9999
500  PRINT 
510  CHAIN X5=T[9]
515  PRINT  USING 3005
520  GOSUB 3500
521  PRINT 
522  P9=J9
523  IF P9=2 THEN 550
524  PRINT  USING 3190
525  PRINT  USING 3200
530  PRINT  USING 2890
531  PRINT  USING 2900
532  PRINT  USING 2905
540  PRINT  USING 2840
550  REM-START OF MAIN PROGRAM CALCULATION VARIABLES
560  REM- EACH 2 SECS FOR 1 MINUTE.
570  FOR M2=1 TO 30
580  A5=T[9]*(40+T[1])*.017
590  X4=X5*(40+T[1])*.017
600  GOSUB 5060
610  J3=J3+2
620  J5=J5+1
630  REM- INSERT RANDOMIZATION OF SYSTEMIC ARTERIAL PRESSURE(S1)
635  A9=.97
636  B9=1.03
640  GOSUB 8070
645  A2=Z9
650  S1=S1*A2
660  IF (J5<1) THEN 680
670  GOSUB 7070
680  IF (M2 >= 30) THEN 1410
690  REM-******DYING ROUTINE
700  IF (T[9]>.1) THEN 730
710  T[9]=0
720  GOTO 850
730  IF (S1 >= 53) THEN 770
740  REM-THIS GRADUALLY AND PROPORTIONATELY REDUCES CARDIAC 
745  REM- FUNCTION EACH 2 SEC.
750  T[9]=T[9]*(.947+.001*S1)
760  GOTO 850
770  IF T[9] <= X4 THEN 810
780  T[9]=X4
790  GOTO 850
800  REM- THIS GRADUALLY BRINGS CARDIAC FUNCTION BACK IF IT IS DOWN 
810  T[9]=T[9]*(1+.0005*(S1-53))
820  J4=1
830  REM- THIS STARTS A COUNTER TO FIND HOW LONG S1 IS LESS THAN 20 
840  REM- AND IF TOO LONG GOES INTO DYING ROUTINE 
850  IF (S1 >= 20) THEN 1170
860  J4=J4+1
870  IF (J4<35) THEN 1170
880  I4=0
890  I6=0
910  PRINT  USING 3380
920  E1=C2-T[5]
930  IF (T[7] >= 3300) THEN 960
940  PRINT  USING 3420
960  IF (E1<25) THEN 990
970  PRINT  USING 3250
990  IF (T[6]>1.6) THEN 1020
1000  PRINT  USING 3430
1020  T[5]=T[5]/2
1030  T[6]=0
1040  A5=0
1050  T[4]=0
1060  S3=0
1070  H1=0
1080  GOSUB 5060
1090  IF (X5>.3) THEN 1120
1110  PRINT  USING 3260
1120  PRINT  USING 2760;I6
1121  PRINT  USING 2770;I4,S1
1122  PRINT  USING 2771
1123  PRINT  USING 2780;V2
1124  PRINT  USING 2785
1125  PRINT  USING 2790;C2,C3
1126  PRINT  USING 2795
1127  PRINT  USING 2800;S3,H1
1130  PRINT  USING 2810;D1
1131  PRINT  USING 2820;T[8],A5
1132  PRINT  USING 2830
1135  PRINT 
1140  PRINT  USING 2920
1141  PRINT  USING 2930;T[2]
1142  PRINT  USING 2935
1143  PRINT  USING 2940;T[3]
1144  PRINT  USING 2945;T[4]
1145  PRINT  USING 2950;T[5]
1146  PRINT  USING 2960;T[6],T[7]
1150  PRINT 
1160  GOTO 230
1170  IF (C3 >= 1) THEN 1190
1180  J4=J4+1
1190  S1=S1-C5
1200  GOSUB 6050
1210  S1=S1+C5
1220  REM-6050 HAS PRODUCED A PREDICTED ARTERIAL RESISTANCE, BUT TO
1230  REM-GIVE A REALISTIC TIME TO WORK,INSERT A TIME DELAY OF 8 
1240  REM-"SECONDS" BY FEEDING VALUES THROUGH A DELAY LINE, AS FOLLOWS 
1250  S[3]=S[2]
1260  S[2]=S[1]
1270  S[1]=A3*A2
1280  REM- ADJUST ARTERIAL RESISTANCE TOWARDS ITS RECALCULATED VALUE 
1290  REM- DERIVED FROM FIGURES SUPPLIED BY THE STABILIZING NETWORK.
1300  REM- BUT GRADUALLY WITH A TIME FACTOR SUPPLIED BY THE INDEX U1 
1310  T[1]=T[1]+(S[3]-T[1])*U1
1320  REM- T(2) IS THE PERCENTAGE ARTERIAL RESISTANCE ADJUSTMENT 
1330  REM- WHICH CAN BE MANUALLY SUPERIMPOSED ON THE STABLIZING
1335  REM- SYSTEM. 
1340  D1=T[2]*.01*(12.18+T[1]/5+C5*.2)
1350  REM- VENOUS CAPACITANCE IS ALSO AFFECTED BY ARTERIAL PRESSURE
1360  REM-AND IS MADE TO CHANGE OVER TWOFOLD RANGE IN RECIPROCAL 
1370  REM- ADJUSTMENT BY ARTERIAL RESISTANCE 
1380  C1=.48*(310+1500/(T[1]+4))
1385  NEXT M2
1390  REM- APPLY TESTS FOR OUTPUT OF SYMPTOMS
1400  REM- FIND EFFECTIVE PRESSURE ACTING ACROSS PULMONARY CAPILLARIES 
1410  E1=(V2+(S1-V2)*.1)-T[5]
1415  PRINT 
1420  K1=0
1430  K3=0
1440  IF E1 >= 25 THEN 880
1460  IF S1>35 THEN 1500
1470  PRINT  USING 3270
1480  K1=1
1490  GOTO 1990
1500  IF (E1<21) THEN 1540
1510  PRINT  USING 3370
1520  K1=1
1530  GOTO 1570
1540  IF (E1<17) THEN 1570
1550  PRINT  USING 3290
1560  K3=1
1570  IF (C2<21) THEN 1610
1580  PRINT  USING 3240
1590  K1=1
1600  GOTO 1640
1610  IF (C2<18) THEN 1640
1620  PRINT  USING 3410
1630  K3=1
1640  IF C3>2 THEN 1680
1650  PRINT  USING 3300
1660  K1=1
1670  GOTO 1740
1680  IF (C3>3) THEN 1710
1690  PRINT  USING 3390
1691  PRINT  USING 3400
1700  K3=1
1710  IF (C3<8) THEN 1740
1720  PRINT  USING 3310
1730  K3=1
1740  IF (S1>51) THEN 1780
1750  PRINT  USING 3350
1751  PRINT  USING 3360
1760  K1=1
1770  GOTO 1810
1780  IF (S1>75) THEN 1810
1790  PRINT  USING 3340
1800  K3=1
1810  IF (S1<160) THEN 1850
1820  PRINT  USING 3330
1830  K1=1
1840  GOTO 1880
1850  IF (S1<145) THEN 1880
1860  PRINT  USING 3320
1870  K3=1
1880  IF (K1 <= K3) THEN 1900
1890  K3=1
1900  K5=K1+K3
1910  IF (K2 <= K5) THEN 1940
1920  PRINT  USING 3440
1921  PRINT  USING 3450
1930  GOTO 1990
1940  IF (K2 <= K1) THEN 1970
1950  PRINT  USING 3460
1951  PRINT  USING 3470
1960  GOTO 1990
1970  IF (K4 <= K3) THEN 1990
1980  PRINT  USING 3480
1981  PRINT  USING 3490
1990  REM
2000  IF C3 <= 1 THEN 880
2010  U2=S1/D1
2020  IF (U2>3.3) THEN 2040
2030  PRINT  USING 3280
2040  K2=K1
2050  K4=K3
2060  REM- CALCULATE STROKE VOLUME FROM CARDIAC OUTPUT AND HEART RATE
2070  S3=C3*1000/H1
2080  REM******OUTPUT VARIABLES VALUES 
2090  PRINT  USING 2760;I6
2091  PRINT  USING 2770;I4,S1
2092  PRINT  USING 2771
2093  PRINT  USING 2780;V2
2094  PRINT  USING 2785
2095  PRINT  USING 2790;C2,C3
2096  PRINT  USING 2795
2097  PRINT  USING 2800;S3,H1
2100  PRINT  USING 2810;D1
2101  PRINT  USING 2820;T[8],A5
2102  PRINT  USING 2830
2105  PRINT 
2110  PRINT  USING 2920
2111  PRINT  USING 2930;T[2]
2112  PRINT  USING 2935
2113  PRINT  USING 2940;T[3]
2114  PRINT  USING 2945;T[4]
2115  PRINT  USING 2950;T[5]
2116  PRINT  USING 2960;T[6],T[7]
2145  PRINT 
2150  PRINT  USING 2980
2151  GOSUB 3500
2160  GOTO J9 OF 2460,2240
2180  REM******NO CHANGE TO BE MADE
2190  PRINT  USING 3000
2215  GOSUB 3500
2220  IF J9=2 THEN 9999
2225  PRINT 
2226  PRINT 
2227  GOTO 515
2230  REM******CHANGES TO STABILIZING SYSTEM 
2240  IF (J7>0) THEN 2190
2250  PRINT  USING 3040
2251  PRINT  USING 3050
2255  GOSUB 3500
2260  GOTO J9 OF 2280,2190
2280  PRINT  USING 3060
2281  PRINT  USING 3070
2285  GOSUB 3500
2290  GOTO J9 OF 2310,2370
2310  FOR I1=1 TO 20
2320  A[2,I1]=40
2325  NEXT I1
2330  J7=1
2340  C5=0
2345  PRINT 
2346  PRINT 
2350  PRINT  USING 3120
2360  GOTO 515
2370  IF (C5 >= 1) THEN 2190
2380  PRINT  USING 3080
2381  PRINT  USING 3090
2385  GOSUB 3500
2390  GOTO J9 OF 2410,2190
2410  C5=45
2415  PRINT 
2416  PRINT 
2420  PRINT  USING 3130
2421  PRINT  USING 3140
2430  GOTO 515
2450  REM******CHANGES TO BE MADE
2460  PRINT  USING 3010
2485  INPUT M1
2488  IF M1=0 THEN 2680
2500  REM******ASK FOR CHANGED VALUES
2530  IF M1<0 OR M1>6 THEN 2645
2550  PRINT  USING 3020;M1
2570  INPUT T[M1+1]
2571  GOTO M1 OF 2572,2572,2572,2580,2590,2596
2572  IF T[M1+1] >= 20 AND T[M1+1] <= 300 THEN 2600
2573  PRINT "MEDICALLY UNACCEPTABLE VALUE, BE MORE REASONABLE"
2574  GOTO 2570
2580  IF T[5]<0 AND T[5]>-6 THEN 2600
2585  GOTO 2573
2590  IF T[6] >= 4 AND T[6] <= 20 THEN 2600
2595  GOTO 2573
2596  IF T[7]>2499 AND T[7] <= 7000 THEN 2600
2597  GOTO 2573
2600  IMAGE "ERROR IN ",3D,"TH FACTOR - REENTER ALL FACTOR NUMBERS"
2610  REM-THIS DETECTS WHETHER MANUAL ADJUSTMENT OF CARDIAC
2620  REM- CONTRACTILITY HAS BEEN MADE OR NOT. IS SO THE VALUE IS
2630  REM- STORED FOR READJUSTMENT WHENEVER ARTERIAL PRESSURE IS NOT 
2635  REM-TOO HIGH OR TO LOW TO ALLOW IT.
2640  IF (M1-3)<0 THEN 2670
2641  IF (M1-3)=0 THEN 2650
2642  IF (M1-3)>0 THEN 2670
2645  PRINT  USING 2600;M1
2646  GOTO 2485
2650  X5=T[4]*.01*W[9]
2660  T[9]=X5
2670  T[8]=T[3]*W[8]*.01
2675  GOTO 2460
2680  PRINT 
2685  PRINT 
2690  PRINT  USING 3160;T[2],T[3],T[4],T[5],T[6],T[7]
2700  GOTO 515
2760  IMAGE #,"ARTERIAL (AORTIC) PRESSURE"/"SYSTOLIC=",3D.D,"   DIAST"
2770  IMAGE #,"OLIC=",3D.D3X,"MEAN=",6D.D," MM HG"/"CENTRAL VENOUS"
2771  IMAGE #," (RIGHT"
2780  IMAGE #," ATRIAL) PRESSURE=",6D.D," MM HG"/"MEAN CAPILLARY"
2785  IMAGE #," PRESSURE="
2790  IMAGE 6D.D," MM HG"/"CARDIAC OUTPUT=",5D.D," L/MIN"
2795  IMAGE #,"STROKE VOLUME="
2800  IMAGE 6D.D," ML/BEAT     HEART RATE=",3D.D," BEATS/MIN"
2810  IMAGE #,"ARTERIAL RESISTANCE=",6D.D7X,"VENOUS RESISTANCE=",5D.D
2820  IMAGE #,5D.D," MM HG/L/MIN"/"CARDIAC CONTRACTILITY=",5D.D
2830  IMAGE " L/MIN/MM HG"
2840  IMAGE "SECS ",14(".    ")
2870  IMAGE "DO YOU WANT ANOTHER SUBJECT?" 
2890  IMAGE#,"TIME 0   20   40   60   80   100  "
2900  IMAGE #,"120",2X,"140",2X,"160",2X,"180",2X,"200",2X,"220",2X
2905  IMAGE"240  260"
2920  IMAGE #,"FINAL VALUES FOR THIS RUN WERE:"/"1. ARTERIAL RE"
2930  IMAGE #,"SISTANCE=",4D," % OF NORMAL"/"2. VENOUS "
2935  IMAGE #,"RESISTANCE="
2940  IMAGE #,4D," % OF NORMAL"/"3. CARDIAC CONTRACTILITY="
2945  IMAGE #,4D," % OF NORMAL"/"4. MEAN INTRATHORACIC PRESSURE="
2950  IMAGE #,4D.D," MM HG"/"5. LIMITING CARDIAC INPUT PRESSURE="
2960  IMAGE 6D.D," MM HG"/"6. BLOOD VOLUME=",5D," ML"
2980  IMAGE #,"DO YOU WANT TO CHANGE ANY OF THESE FACTORS"
3000  IMAGE #,"DO YOU WANT TO CONTINUE"
3005  IMAGE #,"DO YOU WANT A PLOT"
3010  IMAGE #,"CHANGE FACTOR (1..6 OR 0 TO QUIT)"
3020  IMAGE #,"FACTOR ",D
3030  IMAGE "A NUMBER MUST BE ENTERED"
3040  IMAGE "DO YOU WANT TO CHANGE THE OPERATION OF THE SYSTEMIC"
3050  IMAGE #,"ARTERIAL BAROCEPTORS"
3060  IMAGE "DO YOU WANT TO CUT THE BUFFER NERVES AND PUT THE B.P."
3070  IMAGE #,"STABILIZING SYSTEM PERMANENTLY OUT OF ACTION"
3080  IMAGE "DO YOU WANT TO RESET BAROCEPTOR SENSITIVITY TO THAT"
3090  IMAGE #,"EXPECTED OF A HYPERTENSIVE PATIENT"
3100  IMAGE " "
3110  IMAGE  "NO CHANGE"
3120  IMAGE "YOU HAVE CUT THE BUFFER NERVES"
3130  IMAGE "YOU HAVE ALTERED THE STABILIZING SYSTEM TO THAT OF A"
3140  IMAGE "HYPERTENSIVE PATIENT"
3150  IMAGE "   FACTORS ",6(10D)
3160  IMAGE "1=",4D,", 2=",4D,", 3=",4D,", 4=",4D,", 5=",4D,", 6=",4D
3190  IMAGE #,"(XXX) SYSTOLIC/DIASTOLIC PRESSURES - MM HG   "
3200  IMAGE "(.) HEART RATE - BEATS/MIN"
3210  IMAGE #," ALTERED TO ",6(10D.2DXE)
3230  IMAGE " "
3235  IMAGE #,"*"
3240  IMAGE ">>> MY LEGS ARE LIKE TREE TRUNKS"
3250  IMAGE "NECROPSY SHOWED GROSS PULMONARY EDEMA"
3260  IMAGE "NECROPSY SHOWED SEVERE HEART DISEASE."
3270  IMAGE "YOUR PATIENT HAS (MERCIFULLY) LOST CONSCIOUSNESS."
3280  IMAGE "YOUR PATIENT IS PASSING HARDLY ANY URINE."
3290  IMAGE ">>> I AM RATHER SHORT OF BREATH"
3300  IMAGE ">>> I FEEL REALLY TERRIBLE"
3310  IMAGE ">>> MY HEART IS POUNDING"
3320  IMAGE ">>> I HAVE A HEADACHE"
3330  IMAGE ">>> I HAVE A SPLITTING HEADACHE"
3340  IMAGE ">>> I FEEL FAINT WHEN I TRY TO STAND UP"
3350  IMAGE #,">>> I FEEL AS IF I AM GOING TO DIE. MY EYES ARE GOING "
3360  IMAGE "DARK."/">>> WHAT IS HAPPENING TO ME?"
3370  IMAGE ">>> I CAN HARDLY BREATHE AT ALL"
3380  IMAGE "YOUR PATIENT HAS DIED,"
3390  IMAGE #,">>> I AM VERY TIRED AND EXHAUSTED. I DON'T FEEL"
3400  IMAGE " LIKE DOING ANYTHING."
3410  IMAGE ">>> MY LEGS ARE SWELLING UP"
3420  IMAGE"APPARENTLY FROM HEMORRAGIC SHOCK."
3430  IMAGE "NECROPSY SHOWED PERICARDIAL TAMPONADE."
3440  IMAGE ">>> GOD BLESS YOU, DOCTOR. I FEEL REALLY WELL AGAIN."
3450  IMAGE ">>> IT'S LIKE A MIRACLE!"
3460  IMAGE ">>> I FEEL BETTER, BUT NOT RIGHT YET."
3470  IMAGE ">>> CAN'T YOU DO SOMETHING ELSE FOR ME?"
3480  IMAGE ">>> THAT'S BETTER DOCTOR, BUT ARE YOU GOING TO DO ANY"
3490  IMAGE ">>> OTHER NASTY THINGS TO ME?"
3500  REM-INPUT 'YES' OR 'NO' RETURNS J9=1 FOR 'YES',J9=2 FOR 'NO'
3502  J9=1
3504  INPUT A$
3505  IF A$[1,1]="n" THEN 3514
3506  IF A$[1,1]="N" THEN 3514
3507  IF A$[1,1]="y" THEN 3516
3508  IF A$[1,1]="Y" THEN 3516
3510  PRINT "TYPE 'YES' OR 'NO', PLEASE"
3512  GOTO 3502
3514  J9=2
3516  RETURN 
5023  REM
5060  REM- THIS TRIES VARIOUS VALUES OF CARDIAC INPUT PRESSURE (V2)
5070  REM- AND BY TRIAL AND ERROR FINDS THE UNIQUE VALUE TO SATISFY
5080  REM- ALL THE EQUATIONS 
5120  B2=T[5]
5130  REM- IF THE HEART IS ALMOST NON-FUNCTIONAL, THIS 
5135  REM- STOPS IT ALTOGETHER.
5140  IF (A5>.1) THEN 5230
5150  V2=.5*(T[7]-2400)/C1
5160  C2=V2
5170  S1=V2
5180  T[9]=0
5190  C3=0
5200  GOTO 5880
5210  REM- THIS PREVENTS FURTHER INCREMENT OF VENOUS RETURN
5220  REM- WHEN VENOUS PRESSURE IS LESS THAN -2.
5230  T[5]=-2 MAX T[5]
5240  V2=T[5]
5260  REM- NOW FOLLOWS ITERATIVE SOLUTION OF EQUATIONS. VARIOUS VALUES 
5270  REM- FOR VENOUS PRESSURE ARE TRIED. STARTING AT LOWEST 
5280  X7=T[5]
5290  X8=T[6]+T[5]
5300  C3=0
5310  C2=V2
5320  S1=V2
5330  C4=(T[7]-2400)/C1
5340  IF (C4 <= V2) THEN 5880
5350  REM- IF NOT SOLVED YET, VENOUS PRESSURE AT END OF ASCENDING
5360  REM- LIMB OF FUNCTION
5365  V2=X8
5370  FOR I8=1 TO 25
5380  C2=(T[7]-2400)/C1-V2
5390  IF (C2<0) THEN 5460
5400  C3=A5*T[6]
5410  IF (C3>.1) THEN 5430
5420  C3=0
5430  S1=C3*D1+C2+B2
5440  C4=C2-C3*T[8]
5450  IF (C4>V2) THEN 5470
5460  IF (I8>1) THEN 5800
5465  GOTO 5510
5470  V2=V2+1.5
5475  NEXT I8
5480  GOTO 5800
5490  REM- SINCE V2 MUST LIE BETWEEN T(5) AND T(6) THIS FINDS THE RIGHT
5500  REM- VALUE BY INTERPOLATION
5510  FOR I8=1 TO 15
5520  G1=(X7+X8)/2
5530  V2=G1
5540  REM- 2400 IS THE FIXED BLOOD VOLUME, C2 DEEMED TO REPRESENT THE
5550  REM- CAPACITANCE OF THE VENOUS SYSTEM ON WHICH A PRESSURE MADE 
5560  REM- UP BY THE AVERAGE OF VENOUS AND CAPILLARY PRESSURES OPERATES
5570  C2=(T[7]-2400)/C1-V2
5580  REM- CARDIAC OUTPUT DEPENDS ON THE GRADIENT BETWEEN CAPILLARIES
5590  REM- AND VEINS, DIVIDED BY THE VENOUS RESISTANCE 
5610  T[8]=.01 MAX T[8]
5620  C3=(C2-V2)/T[8]
5630  REM- JUST TO MAKE SURE THE WORKS ARE NOT FOULED UP BY
5635  REM- NEGATIVE VALUES
5650  C3=0 MAX C3
5660  REM- ARTERIAL PRESSURE = PRODUCT OF CARDIAC OUTPUT AND PERIPHERAL
5670  REM- ARTERIAL RESISTANCE, PLUS CAPILLARY PRESSURE
5680  S1=C3*D1+C2+B2
5690  REM- CALCULATED CENTRAL VENOUS PRESSURE CAN BE BACK CALCULATED BY
5700  REM- THE STARLING CARDIAC FUNCTION CURVE INTRATHORACIC PRESSURE
5710  REM- MUST BE ADDED TO IT.
5720  C4=C3/A5+T[5]
5730  IF (C4=V2) THEN 5800
5735  IF (C4<V2) THEN 5760
5740  X7=G1
5750  GOTO 5770
5760  X8=G1
5770  NEXT I8
5780  REM- THIS RATHER RAPIDLY LOWERS CARDIAC FUNCTION AS S1 RISES 
5790  REM- ABOVE 160.
5800  IF (S1<160) THEN 5880
5810  IF (S1 >= 250) THEN 5840
5820  F1=-.00012*S1*S1+.0386*S1-2.11
5830  GOTO 5860
5840  F1=0
5850  S1=250
5860  T[9]=T[9]*(99+F1)/100
5870  C3=C3*(.1+F1)/1.1
5880  T[5]=B2
5890  RETURN 
6050  REM
6060  REM-THIS IS A ROUTINE WORKED OUT BY A TYPICAL BARORECEPTOR 
6070  REM-SENSITIVITY CURVE WHOSE VALUES ARE INITIALIZED BY THE DATA 
6080  REM-STATEMENTS IN THE MAIN PROGRAM. A SMOOTH CURVE RELATION
6090  REM-IS OBTAINED BY AN INTERPOLATION METHOD.
6100  REM-GIVEN ARTERIAL PRESSURE THE SUBROUTINE WORKS OUT THE 
6110  REM-CORRESPONDING AND APPROPIATE ARTERIAL RESISTANCE 
6115  REM-(VARIABLE COMPONENT) 
6120  REM
6150  IF S1 <= A[1,1] THEN 6260
6160  IF S1 >= A[1,20] THEN 6280
6170  N1=0
6180  FOR I7=1 TO 20
6190  N1=N1+1
6200  IF (S1=A[1,I7]) THEN 6300
6205  IF (S1<A[1,I7]) THEN 6220
6210  NEXT I7
6220  M3=N1-1
6230  R1=(S1-A[1,N1])/(A[1,N1]-A[1,M3])
6240  A3=A[2,M3]+R1*(A[2,N1]-A[2,M3])
6250  GOTO 6310
6260  A3=A[2,1]
6270  GOTO 6310
6280  A3=A[2,20]
6290  GOTO 6310
6300  A3=A[2,N1]
6310  RETURN 
7070  REM- THIS MAKES A PLOT OF ARTERIAL PRESSURE GRAPHWISE
7075  REM- IN TIME 
7090  DIM B$[72]
7100  B$="                                                                        "
7120  REM-P4 IS A FUNCTION OF ARTERIAL PRESSURE(S1) AND OF STROKE VOLUME
7130  REM-(CARDIAC OUTPUT DIVIDED BY HEART RATE).....THIS FUNCTION 
7140  REM PRODUCES APPROPIATE SPREAD OF DIASTOLIC AND SYSTOLIC PRESSURE
7145  REM- ABOUT THE MEAN
7150  P4=(S1*.09+((C3*1000/H1)-20)*.1)*A2+2
7160  I6=S1+2*P4
7190  I4=S1-P4
7200  P2=I6*.25
7205  P2=P2 MIN 69
7210  P1=I4*.25
7220  REM- ADJUST HEART RATE BY A FACTOR DIRECTLY RELATED TO ARTERIAL
7230  REM- RESISTANCE AND THUS REFLECTING OVERALL SYMPATHO-ADRENAL 
7235  REM- ACTIVITY. 
7240  H1=13.1*(4.2+A5)*(T[1]*(.0006*T[1]+.01)+.6)
7241  A9=.95
7242  B9=1.05
7243  GOSUB 8070
7245  IF P9=2 THEN 7395
7250  P3=H1*.25*Z9
7255  P3=P3 MIN 69
7256  IF P2 >= P1 THEN 7260
7257  B$[1,1]="X"
7258  GOTO 7360
7260  B$[P1,P2]="XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX"
7360  B$[P3,P3]="."
7365  GOSUB 8500
7370  PRINT J3;TAB(6);B$
7390  RETURN 
7395  PRINT  USING 3235
7400  RETURN 
8070  REM-TRI. RANDOM NUNBER GEN.
8130  Z9=.5*(B9-A9)*(RND(0)+RND(0))+A9
8160  RETURN 
8500  P8=((P1 MAX P2) MAX P3)
8505  REM STRIPS TRAILING BLANKS
8510  B$=B$[1,P8]
8520  RETURN 
9999  END 
