with reference to fig1 the apparatus 10 of the present invention is schematically shown as comprising the elements of sensor 11 , computer 12 and pump 13 . sensor 11 is capable of generating computer input signals which are transmitted to computer 12 over line 14 . computer 12 is capable of generating output signals which are transmitted to pump 13 over line 15 . in the operation of this apparatus , blood samples are obtained from patient 16 and are conducted to sensor 11 through conduit 17 . pump 13 is connected to a supply 18 of insulin or glucose through conduit 19 and to the patient 16 through conduit 20 . the sensor 11 determines the blood glucose concentrations of blood samples obtained from patient 16 at convenient intervals , such as once each minute , and transmits such serial signals indicative of such blood glucose concentrations to computer 12 . the output signals from computer 12 are directly related to the serial blood glucose concentrations measured by sensor 11 and provide instructions to pump 13 to supply adequate amounts of insulin or glucose to bring the blood glucose concentration within the desired range . the particular novel aspects of the present invention relate to the algorithms or equations which the computer 12 employs to derive the proper output signals based upon the input signals from sensor 11 as well as to the flexibility of operation of such computer . it has been found through experimentation with blood glucose levels of diabetic patients that improved control over such blood glucose concentration can be achieved through operation in accordance with certain quadratic and biquadratic computer equations . the apparatus of the present invention is capable of operating in several modes at the choice of the operator . in a first mode , the insulin is infused to a patient at a rate dependent upon both the particular concentration of glucose in the blood and the rate of change of blood glucose concentration . in a second mode , the insulin is infused at a rate dependent only upon the particular concentration of glucose in the blood . in a third mode , the insulin is infused at a rate dependent only upon the rate of change of blood glucose concentration . in order to overcome a hypoglycemia condition ( undesirably low blood glucose concentration ) the apparatus also provides infusion of glucose at a rate dependent only upon the particular concentration of blood glucose whenever the blood glucose concentration is below a predetermined value . this latter situation is considered emergency situation to be used primarily when the patient initially is in a hypoglycemis condition . the operation in the first three modes will normally prevent such condition from occurring as a result of infusion of an excess amount of insulin . it could also be used if desired when the operator intends to maintain a &# 34 ; glucose clamping &# 34 ; condition in which the glucose level is maintained at an abnormally high level . in the above first mode of operation , the computer is capable of deriving output signals from serial input signals in accordance with an equation of the form ## equ1 ## wherein g = measured previous serial blood glucose reading and bi , ri and qi are preselected values . if the blood glucose concentration is rising , wherein a is the average change in blood glucose concentration over several previous serial readings and k r is a preselected value and wherein if the blood glucose concentration is falling , wherein a is the same as above and k f is a preselected value which is different from and less than k r , said gd being substracted from g in the above equation when the blood glucose concentration is falling . in the above equation ir and ri are expressed in milliunits of insulin per minute of infusion . bi , qi , g , gd and a are expressed in mg . percent of glucose . this can also be expressed as milligrams per deciliter ( mg . dl ). bi represents a selected basal blood glucose concentration which would be present in a normal patient under resting conditions . ri represents the basal insulin infusion rate normally provided by the body of the normal patient under such conditions . when the apparatus of the present invention is used with a 70 kg . human patient , for example , ri is typically about 9 to 12 milliunits of insulin per min . and bi is about 80 mg . percent glucose . a typical qi is about 85 mg . percent glucose . k r generally is in the range from 60 to 200 and is typically about 125 . k f generally is in the range from 0 to 100 and is typically about 25 . the values for ri , bi , qi , k r and k f are selected by the operator depending upon the specific patient involved and upon previous experience . in the above second mode of operation , the computer is capable of deriving output signals from serial input signals in accordance with an equation of the form ## equ2 ## wherein ir , ri , g , bi and gi are discussed above . in the above third mode of operation , the computer is capable of deriving output signals from serial input signals in accordance with an equation of the form ## equ3 ## wherein ir , ri , g , gd , bi and qi are discussed above . when the apparatus is employed to infuse glucose , the computer is capable of deriving output signals from serial input signals in accordance with an equation of the form ## equ4 ## wherein g = measured previous serial blood glucose reading and rd , bd and qd are preselected values . in the above equation dr and rd are expressed in milligrams per minute of glucose infusion . bd , qd and g are expressed in mg . precent of glucose . bd represents a selected basal blood glucose concentration and rd represents a normal glucose infusion rate to obtain such basal concentration . when the apparatus of the present invention is used with a 70 kg . human patient , for example , rd is about 9 to 11 milligrams of glucose per minute and bd is about 70 mg . percent glucose . a typical qd is about 20 mg . percent glucose . it has been determined that in order to prevent any undersirable &# 34 ; overshooting &# 34 ; of desired blood glucose concentrations by excess infusion or glucose , the infusion rates for insulin and glucose should have controlled maximum values . a maximum insulin infusion rate of about 600 milliunits per minute and a maximum glucose infusion rate of about 200 milligrams per minute have been found suitable . typical operating conditions for computer calculated insulin infusion rates ( ir ) vs . measured serial blood glucose concentrations ( g ) for the second mode of operation are shown in fig2 . the ri and qi are kept constant at typical values of 10 milliunits per minute and 85 mg . percent , respectively . bi is shown for three typical values of 60 , 80 and 100 mg . percent . the calculated ir is then obtained from the appropriate curve 20 depending on the measured blood glucose reading ( g ) and the appropriate bi selected by the operator . comparable curves can also be obtained for different values of bi . it is noted that the curves are truncated at a maximum insulin infusion rate of 600 milliunits per minute . typical operating conditions for computer calculated insulin infusion rate ( ir ) vs . measured serial blood glucose concentrations ( g ) for the second mode of operation are shown in fig3 . the bi and qi are kept constant at typical values of 80 and 85 mg . percent , respectively . ri is shown for four typical values of 5 , 10 , 20 and 40 milliunits of insulin per minute . the calculated ir is then obtained from the appropriate curve 21 depending on the measured blood glucose reading ( g ) and the appropriate ri selected by the operator . comparable curves can also be obtained for different values of ri . typical operating conditions for computer calculated insulin infusion rates ( ir ) vs . measured serial blood glucose concentrations ( g ) for the second mode of operation are shown in fig4 . the bi and ri are kept constant at typical values of 80 mg . percent and 10 milliunits per minute , respectively . qi is shown for three typical values of 50 , 85 and 120 mg . percent . the calculated ir is then obtained from the appropriate curve 22 depending on the measured blood glucose reading ( g ) and the appropriate qi selected by the operator . comparable curves can also be obtained for different values of qi . typical operating conditions for computer calculated gd values to be used in the appropriate equations vs . rate of change of blood glucose concentration ( a ) for the first and third modes of operation are shown in fig5 . the effect of variations in k are known by the different curves 23 . if the blood glucose concentration is rising , k refers to k r . if the blood glucose concentration is falling , k refers to k f . comparable curves can also be obtained for different values of k . in the operation of the apparatus in the various modes for insulin infusion , the portions of the equations relating to ( g - bi qi ) + 1 and ( g + gd - bi qi ) + 1 must have values greater than zero . if either of these portions has values less than zero , the computer will stop all insulin infusion . in the event that the measured serial blood glucose concentrations drop below about 70 mg . percent , it is desirable to infuse glucose to bring the blood glucose concentration within the desired physiological range . in this case the computer will automatically switch to the glucose infusion mode described above . glucose infusion according to this mode can also be accomplished at any time at the discretion of the operator irrespective of the specific blood glucose concentration . typical operating conditions for computer calculated glucose infusion rates ( dr ) vs . measured serial blood glucose concentrations ( g ) are shown in fig6 . the qd and rd are kept constant at typical values of 35 mg . percent and 10 mg . per min . respectively . bd is shown for three typical values of 60 , 80 and 100 mg . percent . the calculated dr is then obtained from the appropriate curve 24 depending on the measured blood glucose reading ( g ) and the appropriate bd selected by the operator . comparable curves can also be obtained for different values of bd . it is noted that the curves are truncated at a maximum glucose infusion rate of 200 milligrams per minute . typical operating conditions for computer calculated glucose infusion rates ( dr ) vs . measured serial blood glucose concentrations ( g ) are shown in fig7 . the bd and qd are kept constant at typical values of 80 and 35 mg . percent , respectively . rd is shown for four typical values of 5 , 10 , 20 and 40 mg . per min . the calculated dr is then obtained for the appropriate curve 25 depending on the measured blood glucose reading ( g ) and the appropriate rd selected by the operator . comparable curves can also be obtained for different values of rd . typical operating conditions for computer calculated glucose infusion rates ( dr ) vs . measured serial blood glucose concentrations ( g ) are shown in fig8 . the bd and rd are kept constant at typical values of 80 mg . percent and 10 mg . per min ., respectively . qd is shown for three typical values of 20 , 35 and 50 mg . percent . the calculated dr is then obtained for the appropriate curve 26 depending on the measured blood glucose reading ( g ) and the appropriate qd selected by the operator . comparable curves can also be obtained for different values of qd . when glucose is being infused , the equation portion relating to ( bd - g gd ) + 1 must have a value greater than zero . if this value is less than zero , the computer will stop all glucose infusion . it can be seen from the above discussion that the apparatus of the present invention provides considerable flexibility to the operator to select not only the specific desired mode of operation but also to select desired control values in the computer operating equation for the selected operating mode . it is possible , when desired to have the apparatus simultaneously infusing both insulin and glucose according to separate different operating modes if this is necessary to maintain proper control over the blood glucose concentration of the patient . the operator can thus program the apparatus to most closely match the physiological condition of the particular patient and most effectively control the blood glucose concentration of this patient at a predetermined level . the apparatus of the present invention has been used in conjunction with several diabetic patients and has successfully maintained their blood glucose concentrations at desired levels under both stressed and unstressed situations and with various diets . the degree of control achieved has been superior to that obtained by prior art apparatus . in the apparatus of fig1 any suitable sensor means can be employed to measure the blood glucose concentration of the blood samples obtained from the patient . likewise any suitable pump can be used to infuse insulin or glucose to the patient in response to signals from the computer . the particular details of the sensor 11 , the computer 12 or the pump 13 do not form a part of the present invention . it is understood that for convenience there will be separate pumps for infusion of insulin and for infusion of glucose .
apparatus is described for controlling blood glucose concentration in a subject by selected infusion of insulin and / or glucose , depending upon the general blood glucose concentration . this apparatus comprises in combination sensor means for measuring the blood glucose concentration and for providing computer input signals based on such measurement , computer means for receiving such input signals and being programmed to provide output signals based upon such input signals , and pump means responsive to such computer output signals to supply insulin and / or glucose to such subject at a rate determined by such output signals . the computer derives the output signals in accordance with specific quadratic and biquadratic equations employing certain independently selected values relating to basal blood glucose concentrations and basal insulin and glucose infusion rates at such basal blood glucose concentrations .