ALBERT

Notes: Abstract It is observed that a dynamical continuity equation for biomass distribution yields the asymptotic steady-state exponential dependencen=A exp( $$ - m/\bar m$$ ) exhibited by certain fishery data, wherem is the biomass of an individual,n is the number of individuals per unit biomass interval, andA, $$\bar m$$ are positive constants. This dynamical approach to biomass distribution is an alternative to the global maximization principle proposed recently by Lurié and Wagensberg.

Notes: Abstract Biological tree-like structures, such as mammalian tracheobronchial airways, are complicated branching systems. One problem in modeling such systems is the reassignment of the number of segments at a given generation in the model being constructed. A hypothesis is proposed which has successfully been used in modeling mammalian lung airways.

Notes: Abstract The theory of complementary variational principles is used to obtain maximum and minimum principles for diffusion problems with Michaelis-Menten kinetics. In an illustrative calculation we obtain an extremely accurate variational solution in good agreement with the numerical solution of McElwain (1978).

Notes: Abstract Necessary and sufficient conditions for primitivity of a product of two Leslie matrices are given. Such a product could be used in modeling the growth of a population governed alternately by two different sets of fertility and survival parameters.

Notes: Abstract Zadeh's transfer function method for linear time-variable systems is used to apply frequency-domain analysis to a periodically time-varying elastance model of the left ventricle. Left ventricular pressure computed from the system function of the time-varying elastance and the phasors of aortic flow shows a typical waveform of the measured ventricular pressure.

Notes: Abstract By studying the behavior of various tracer species in the lungs, one can assess many important characteristics which distinguish normal and abnormal function. Quantitative evaluation of function depends on the use of an appropriate model in conjunction with experimental data. A multi-compartment model is derived from mass balances to describe dynamic as well as (breath-averaged) steady-state transport processes between the environment and pulmonary capillary blood. The breathing cycle is divided into three time periods (inspiration, expiration, and pause) so that the model equations are discrete in time. No other model of tracer species transport in the lungs deals simultaneously with species dynamics, variable breathing pattern, distribution inhomogeneities, and non-equilibrium between alveolar gas and capillary blood. Models currently in the literature are shown to be special cases of the model presented here.

Notes: Abstract Local stability seems to imply global stability for population models. To investigate this claim, we formally define apopulation model. This definition seems to include the one-dimensional discrete models now in use. We derive a necessary and sufficient condition for the global stability of our defined class of models. We derive an easily testable sufficient condition for local stability to imply global stability. We also show that if a discrete model is majorized by one of these stable population models, then the discrete model is globally stable. We demonstrate the utility of these theorems by using them to prove that the regions of local and global stability coincide for six models from the literature. We close by arguing that these theorems give a method for demonstrating global stability that is simpler and easier to apply than the usual method of Liapunov functions.

Notes: Abstract The asymptotic behaviour of a logistic equation with diffusion on a bounded region and a diffusionally coupled delay is investigated. An equivelent parabolic system is derived for certain types of delays. Using a Layapunov functional, sufficient conditions for the global asymptotic stability of the constant steady state are obtained. When the global stability is lost, using Hopf's bifurcation theory, existence of travelling waves is shown for ring-like and periodic one dimensional habitats.

Notes: Abstract It was hypothesized in an earlier work that sensory perception can occur only when the perceiving system is uncertain about the nature of the event being perceived. In the absence of any uncertainty, perception will not take place. The response of the sensory afferent neuron (impulse transmission rate) was calculated using Shannon's measure of uncertainty or entropy. It will now be shown that when the event being perceived is the position and momentum of a particle, Shannon's measure of uncertainty leads to the Heisenberg Uncertainty relationship.

Notes: Abstract It is not unusual for several classifications to be given for the same collection of objects. We present a method, called majority rule, which can be used to define a consensus of these classifications. We also discuss some mathematical properties of this consensus tree.

Notes: Abstract The dependence of the spatial concentration profiles of morphogens on a characteristic dimension is obtained by continuation techniques for Gierer and Meinhardt's activator-inhibitor model of morphogenesis. The study of the behaviour of the system during growth, where the linear and exponential increase of the characteristic dimension is considered, revealed that more complex patterns of morphogen spatial concentrations appear regularly in a reproducible way.

Notes: Abstract Computer models have been used by various authors to simulate both the growth of normal cellular tissue and the development of cancerous cells within normal tissue. As these models were the result of considerable idealization, the purpose of the present paper is to propose a model in which the degree of simplification is relaxed: the features of simultaneous growth, and cell growth whose rate depends on the free absorbing periphery of the cell are introduced. Simulation experiments have been conducted using the model, and the results are presented.

Notes: Abstract The theory of complementary variational principles is used to obtain maximum and minimum principles for a nonlinear model of heat conduction in the human head. Accurate variational solutions are obtained in illustrative calculations. The effect of nonlinearity is seen to be significant from a comparison with the linearized model.

Notes: Abstract A model for the nerve impulse due to Zeeman (1972) and based on catastrophe theory is compared with alternative models and criticisms of Zeeman's model by Sussmann and Zahler (1977, 1978) are assessed. The criticisms of Zeeman's motivation for his model are found to carry some weight. Sussmann and Zahler (1977, 1978) list numerous features of Zeeman's model which, they state, are not in agreement with experiment. These statements as they stand are largely erroneous, and the model still remains to be tested by a critical series of experiments. However, a detailed analysis reveals defects in Zeeman's model, not among those claimed by Sussmann and Zahler, showing that the explicit equations of the model cannot be correct. The possibility of a modified approach along similar lines and its ultimate adoption remains open.

Notes: Abstract The irreversible Michaelis-Menten reaction is studied by the use of the method of multiple scales. Three stages of the reaction are identified, one of which is studied in detail. The results are compared with those of two earlier analyses.

Notes: Abstract A numerical study of the coupled nerve fibre problem is given which verifies and extends the perturbation theory of Luzader. Pulses on adjacent fibres can couple together with two possible stable pulse separations.

Notes: Abstract A possible mechanism for effects of microwave radiation on the auditory system is the generation of field-induced forces at interfaces that divide materials of dissimilar electrical properties. A general expression for these “Maxwell stresses” is derived and then used to calculate the approximate magnitude of field-induced force within the organ of Corti during microwave exposure. Comparison of the results with data on the force needed to excite cochlear hair cells indicates auditory responses could be evoked by this mechanism at power densities near the threshold of rf hearing sensations.

Notes: Abstract A definition of homogeneity for neural networks is given which permits their construction as group quotients. The significance of this for neural dynamics is discussed.

Notes: Abstract The left ventricle is represented as a cylinder contracting both radially and longitudinally. A simple method is indicated to derive an expression for the rate of change of the kinetic energy of this three-dimensional model, which quantity can be used as an index for the study of the contractile behaviour of the myocardium. An application to the study of muscle mechanics is also indicated.

Notes: Abstract A perturbation method is proposed to calculate approximately the limit cycle type nonequilibrium steady-state resulting from periodic perturbation of coefficients of stable population systems; the two species Lotka-Volterra competition system is explicity studied and the results are formulated for general multi-species population systems. Avoidance of competitive or other types of exclusion of species in a periodic environment is indicated.

Notes: Abstract A theorem is proved, concerning expected values of a multitype branching process in a varying environment. The consequence of the theorem is that the branching process can be treated (in the sense of expected values) as a dynamical system with control terms. This is of importance in situations where the process serves as an abstract model of the dynamics of malignant cells for use in chemotherapy. A simple example of this kind is given.

Notes: Abstract A method is developed for a full nonlinear evaluation of all velocities and stresses represented in the Navier-Stokes equations and in the general stress tensor. The information required is essentially that for solution of linearized forms. The solution is analytical except for the calculation of the axial velocity, which requires computer assistance to step through time and space. The treatment of the problem, although directed towards solutions involving fluid flow in elastic vessels, is also adaptable to solid deformations (strain vs rate of strain) where the general stress tensor applies. A special case for the distorting ellipse is presented as well as a limited, spatially analytic, solution.

Notes: Abstract The stability characteristics and dynamical behavior of a system of mutually excitatory neurons in close spatial proximity are investigated with a mathematical model. The model predicts the existence of uniform, intermediate levels of activity other than those of no activity and maximal activity. The model also, yeilds a good explanation of data obtained from periglomerular neurons in the olfactory bulb of the cat.

Notes: Abstract A multicompartmental model in which particles enter the system from the environment and reproduce according to a Markov branching process has been considered. Explicit expressions have been obtained for the mean vector and the correlation structure for the numbers of particles in different compartments in different time points of the system. Growth rates of the mean vector and some special cases of the system are also discussed.

Notes: Abstract The study of bi-directionally coupled oscillators is relevant in biological modelling of such systems as gastro-intestinal electrical activity, cardiac pacemarkers, cardiovascular and respiratory interactions and circadian rhythms. Interconnecting pathways in biological systems often exhibit pure time-delay characteristics. In this paper the multiple-mode limit-cycle behaviour of such systems is analysed using the method of harmonic blance. It is shown that the coupling time delay radically affects the number, frequency and amplitudes of entrained limit-cycles.

Notes: Abstract The effect of keeping all the parameters constant, except the diffusion coefficients, in a pair of reaction-diffusion equations is studied. It is shown that the stability of the constant solution and the bifurcation points can be easily established by constructing a simple stability diagram. The possible qualitatively different diagrams are enumerated.

Notes: Abstract We have numerically examined more than one million Large Complex Systems (LCS) of interacting variables (interpretable as interacting populations) governed by Generalized Lotka-Volterra Equations (GLV), with self-regulation term. The scope was to have some insight on the stability-complexity relationship. We considered systems of prey-predator type, and we gave appropriate rules for constructing the model systems, rules that specify the behaviour of model systems in order to put them near the biological reality. The results show, among other things, a strict correlation between the stability and the prey-predator ratio (which, in our model, uniquely determines the connectedness of the system).

Notes: Abstract The rate-controlling process in the oxygenation of red blood cells is investigated using a Roughton-like model for oxygen diffusion and reaction with hemoglobin. The mathematical equations describing the model are solved using two independent techniques, numerical inversions of the Laplace transform of the equations and numerical solutions via an implicit-explicit finite difference form of the equations. The model is used to re-examine previous theoretical models that incorporate either a red cell membrane that is resistive to oxygen diffusion or an unstirred layer of water surrounding the cell. Although both models have been postulated to be equivalent, the results of the computer simulations demonstrate significant differences between the two models in the rate of oxygenation of the red cells, depending upon the values chosen for the diffusion coefficient for O2 in the membrane and the thickness of the water layer. The difference is apparently due to differences in the induction and transient periods of the water layer model relative to the membrane model.

Notes: Abstract The chromosomal theory of inbreeding based on a gametic interaction system lead us to define a depression coefficientD. Comparison of random, sib and half-sib matings (with inbreeding coefficientF=0, 1/4 and 1/8) shows thatD depends on the structure of the starting population and on values of the model parameters. This result accounts for responses of lines whose depression does not depend directly on the inbreeding coefficient and which theories of inbreeding based on increasing homozygosity fail to explain.

Notes: Abstract An idealization of chemical combination is formulated as a model of computability, and it is shown that this model has universal computational power just in case assembly has at least two-dimensional space in which to occur. It is also shown that this model, under reinterpretation, corresponds to a cellular automaton in which growth occurs by differentiation only (i.e., the state into which any cell is born is thereadfter fixed). Hence this latter model of growth is also computationally universal.

Notes: Abstract Kinetics of biological light emission processes do not mean what they seem to mean, because measured light intensity is not proportional to reactant concentration but to reaction rate. Therefore, the differential equation for light decay is usually different from that of concentration decay, so that mass action interpretations cannot be applied directly to light intensity decay. An observed second order light decay for Chlorella at 6.5°C, implies Elovich solid state reaction kinetics, which agrees with other evidence for solid state processes in photosynthesis. An observed 1.5 order light decay for Cholorella at 28°C implies second order liquid or solid state reaction kinetics. First ordere light decay implies first order reaction kinetics.

Notes: Abstract Analysis based on the integration of differential inequalities is employed to derive upper and lower bounds on the total populationN(t) = ∫ R θ(x 1,x 2,t) dx 1 dx 2 of a biological species with an area-density distribution function θ=θ(x 1,x 2,t) (≥0) governed by a reaction-diffusion equation of the form ∂θ/∂t =D∇2θ +fθ −gθ n+1 whereD (〉0),n (〉0),f andg are constant parameters, θ=0 at all points on the boundary ∂R of an (arbitrary) two-dimensional regionR, and the initial distribution (θ(x 1,x 2, 0) is such thatN(0) is finite. Forg≥0 withR the entire two-dimensional Euclidean space, a lower bound onN(t) is obtained, showing in particular thatN(∞) is bounded below by a finite positive quantity forf≥0 andn〉1. An upper bound onN(t) is obtained for arbitrary bounded or unbounded)R withn=1,f andg negative, and ∫ R θ(x 1,x 2, 0)2 dx 1 dx 2 sufficiently small in magnitude, implying that the population goes to extinction with increasing values of the time,N(∞)=0. Forg≥0 andR of finite area, the analysis yields upper bounds onN(t), predicting eventual extinction of the population if eitherf≤0 or if the area ofR is less than a certain grouping of the parameters in cases for whichf is positive. These results are directly applicable to biological species with distributions satisfying the Fisher equation in two spatial dimensions and to species governed by certain specialized population models.

Notes: Abstract We examine certain mathematical structures presented in Part I. The most important of these are the energy structures determined by the couple (ω×E, ψ) the space of causality defined by ψ-1(0) and the notion of collapsibility, i.e., the descent of a species from a higher to a lower equilibrium configuration as a result of energy loss.

Notes: Abstract This paper discusses a general stochastic model for a two-compartment reversible system with non-homogeneous Poisson inputs, arbitrary residence times at each of the compartments and time-dependent transition probabilities. The probability distributions of the number of particles in each compartment and in the system are obtained together with the number of particles which depart from the system. In addition, various covariance functions with a time lag are obtained. Some of the above obtained results are deduced for time-independent arrivals, exponential residence times and time-independent transition probabilities. Fluctuations of the particles present in the system are also analysed. Similar analysis is provided for the model into which some particles are initially introduced at the system. Some possible applications are discussed at the end.

Notes: Abstract The general multispecies prey-predator system with Gompertz's antisymmetric interactions is nonlinearly stable in the absence of dispersion and continues to remain stable with dispersion under both homogeneous reservoir and zero flux boundary conditions in a region containing the equilibrium state. It is proved that a general multispecies food-web model without antisymmetric interactions is stable in the absence of dispersion and remains stable with dispersion in the above-mentioned region.

Notes: Abstract We introduce a graphical approach in the study of the qualitative behavior ofm species predator-prey systems. We prove that tree graphs imply global stability for Volterra models and local stability for general models; furthermore, we derive sufficient conditions so that loop graphs imply stability and boundedness of the solutions.

Notes: Abstract A system of integro-differential equations is derived to describe epizootics of a fungal pathogen in an insect population. Because of piecewise continuous behavior under some parametric conditions, it is concluded that standard phase orbits can be misleading. Using a different analytic approach yields a simple system of finite difference equations. Both the continuous and discrete versions are compared to classical forms. The continuous version differs from a classical one in possessing a second derivative dependent on population density. The discrete version differs in maintaining positive, non-zero populations of both infectives and susceptibles in finite time.

Notes: Abstract This note is an attempt to demonstrate that hypothalamic pulsatile GnRH secretion is not the result of a short-term, negative steroid hormone feedback. Clarification of this point is of importance for further modelling the control of gonads.

Notes: Abstract A modern theory of the calculus of variations is used to form necessary and sufficient conditions for the existence of a Lagrangian representation of a system of first-order ordinary differential equations. There exists a theorem to the effect that when a system of ordinary differential equations is variationally self-adjoint, the fulfillment of such conditions is guaranteed. In addition, self-adjointness, allows establishement of an algorithm by which a Lagrangian for the system may be explicitly constructed. Examples in mathematical biology are given to illustrate the use of the stated theorem.

Notes: Abstract Engineering optimal control theory is applied to equations describing insulin and glucose interactions. The nature of the optimal controller is established. It is shown how the results can be utilized in a closed loop feedback control system.

Notes: Abstract An attempt is made to compare the conditions for the general error-optimality of linear systems developed by Kalman with the conditions for feasibility of linear models of neuromuscular and physiological control systems. Models of three actual physiological systems are tested for both the above criteria. Theoretical analysis presented here shows that there are no simple relationships between the two sets of conditions. Analysis carried out on the physiological systems models suggests the need for a general set of conditions for other optimality criteria, such as time and energy minimization, similar to Kalman's condition for error minimization.

Notes: Abstract The equilibrium probability distribution of the process level is studied for a general class of reversible stochastic reactions. A calculationally convenient approximation for equilibrium probabilities is derived and its accuracy is investigated over a range of values of the equilibrium constant. A method of approximating the equilibrium means and variance is developed and illustrated forQ th-order processes.

Notes: Abstract Beside the concept of material inputs and outputs of components of the representation of biological systems given to us by Rosen, the concept of energy is incorporated. The interaction of material and energy is represented by a cartesian product; and separate material and energetical mappings are considered as the new representation of components. These developments generate aMα category, and it is shown thatMα is isomorphic to theM category of previous developments.

Notes: Abstract This paper discusses the solution of a generaln-compartment system with time dependent transition probabilities utilizing the technique described by Cardenas and Matis (1975) (hereafter abbreviated (CM)). In addition, the cumulant generating function is derived for a special class of reversiblen-compartment systems where the time-dependent intensity coefficients corresponding to the migration and death rates are some multiple of each other. The immigration rates can be any integrable function of time. The moments are also obtained and the solution to the two-compartment system is presented explicitly. The solution is illustrated with a linear and a periodic function which forms have been widely reported in the literature.

Notes: Abstract The relations (inflow) = (dose)/(area under indicator curve), and (volume of distribution) = (throughflow) × (mean transit time) are derived by a matrix method for a system of interconnected subsystems, within which spatial indicator activity gradients may exist, and for compartments, within which the indicator activity is spatially uniform. The inflow theorem, is different from the outflow theorem. Equivalent labeling of multi-input systems reduces them formally to single input systems. Foreign indicator flow-volume kinetics are more general than, and include as a special case, tracer flux-mass (metabolic) kinetics. Volume of distribution in the indicator steady state may be different from the equilibrium volume of distribution.

Notes: Abstract Following arteriolar occlusion, tissue oxygen concentration decreases and anoxic tissue eventually develops. Although anoxia first appears in the region most distal to the capillary at the venous end, it eventually spreads throughout the entire region of supply. In this paper the changing oxygen concentration, from the time of occlusion until the tissue is entirely anoxic, is examined mathematically. The equations governing oxygen transport to tissue are solved by iterating a nonlinear integral equation. This solution is valid until anoxia first appears. After anoxia develops it is necessary to solve a moving boundary problem. This is done using the method of matched asymptotic expansions, and accurate solutions are obtained for a wide range of physiological conditions.

Notes: Abstract Perturbation methods are applied to a differential equation predator-prey model to find the approximate amplitudes and period of limit cycle solutions. In the model the feeding rate per unit predator per unit prey decreases as the prey become scare. The rigorous applicability of the perturbation technique depends on the assumptions that the limit cycle amplitude is relatively small and that near the equilibrium point the growth rate of each species is most sensitive to changes in the density of the other species. The second assumption is usually roughly satisfied in practice and examples are considered which suggest that the first assumption can be greatly relaxed.

Notes: Abstract Signal Detection Theory can be used to provide a mathematical model describing the choice of a predator trying to distinguish between a model and a Batesian mimic. The mathematical model yields a number of a deductions, in particular that it may or may not assist the mimic population if mimics more closely resemble their models. The assumptions underlying the analysis are discussed in some detail.

Notes: Abstract A new type of physical transition, denotedS→S *, has been detected in irradiated organic molecules (λ=546 nm) through their interaction with specific biological macromolecules. In a specific enzyme-substrate interaction, a clear enhancement of the reaction rate is observed, when the substrate is irradiated with sharply well defined times. These “efficient irradiation times” are always of the 5k sec type (k=1, 2, 3, …). They have been consistently revealed in a great number of specific biological interactions. The present note demonstrates an important property, i.e. that forevery irradiation time aS→S * transition is induced in organic molecules. It is shown that for any irradiation times different from the 5k sec type (k=1, 2, 3, …) states of theS * type may occur, but the biological macromolecules may “detect” only theS * states induced by irradiations of the 5k sec type.

Notes: Abstract A semi-empirical model applicable to the flow of blood and other particulate suspensions through narrow tubes has been developed. It envisages a central core of blood surrounded by a wall layer of reduced hematocrit. With the help of this model the wall layer thickness and extent of plug flow may be calculated using pressure drop, flow rate and hematocrit reduction data. It has been found from the available data in the literature that for a given sample of blood the extent of plug flow increases with decreasing tube diameter. Also for a flow through a given tube it increases with hematocrit. The wall layer thickness is found to decrease with increase in blood hematocrit. A comparison between the results of rigid particulate suspensions and blood reveals that the thicker wall layer and smaller plug flow radius in the case of blood may be attributed to the deformability of the erythrocytes.

Notes: Abstract Three dimensional laminar, viscid flow is developed for Newtonian fluids which provides absolute values for axial, radial and tangential velocity fields everywhere if the dimensions of the vessel are known and two simultaneous axial velocities e.g. on and off the central axis in the same plane, and the central axis axial velocity gradient are measured. In addition, normal and shear stresses are determinable. The equation set satisfies geometric and other known flow limiting conditions such as no slip at surfaces etc. and are amenable for inclusion in general, dynamic flow expressions. Alternatively they may be used alone for certain problems involving gradients and secondary flows. A range of illustrations are shown for a distorting vessel with elliptic cross-section and small axial taper (analogous to the pulmonary trunk during ejection).

Notes: Abstract A regulated left ventricular dynamics model is presented which involves interaction of the dynamics of the left ventricular and circulatory systems and their regulation by the central nervous system. On-line human parametric simulation (parameter estimation) and consequential prognostic implications (based on parametric values) are demonstrated. Model responses to simulated physiologic stresses help delineate tolerances of subjects. In order to have an estimate of the reliability of the model, the sensitivity of the model's responses to changes in the values of its intrinsic parameters is assessed. Also determined is the extent to which errors in measuring the pressure affect the calculated values of the model's simulation parameters and subsequently influence the values of other diagnostically useful variables (such as contractility, oxygen consumption rate, heart rate), when the model is used to determine the limiting physiological stress sustainable by the subject. A comparison of the model's composition with those of other similar cardio-circulatory models is included.

Notes: Abstract Then-stage harvesting strategy of Elizarov and Svirezhev is examined. As a result, some important new features appear. A discussion is presented on whether or not one should harvest a species at one time stage or wait until a later time. The paper is concerned with contributions which are primarily mathematical formulations and results for continuous, as well as discrete time, logistic growth of a single species being harvested. Age class structure is ignored.

Notes: Abstract In order to evaluate the effect of anatomic asymmetries on the gas concentration distribution in the pulmonary airways, a Monte Carlo simulation of combined bulk flow and molecular diffusion was carried out in a realistic distal airway model (Parkeret al., 1971). This airway model, composed of branches distal to the 0.5-ram diameter airways, contained an upper symmetric segment consisting of four generations of conducting airways and a lower asymmetric segment of alveolar ducts and sacs arranged in five transport paths of varying lengths. In accounting for the volume increases of these ducts and sacs occurring during normal respiration, uniform alveolar filling rates and a fixed length-to-diameter ratio of all airways were assumed. For a pulse injection of inert tracer gas, the simulation was employed to determine the longitudinal concentration profiles in the conducting airways. In the alveolated airways, not only were the longitudinal profiles determined along each path, but radial transport from the core to the periphery of the airways was considered. The results of the simulations indicate that geometric asymmetries alone contribute substantially to regional concentration variations in the distal airways. For example, when a gas bolus is injected at mid*inspiration, there are concentration differences as great as 40% between two points along different transport paths located equi-distant from the proximal end of the model. As viewed from the terminal end of the model (acinus), average concentration differences as large as 6-to-1 exist between the longest and shortest transport paths respectively for gas boli introduced near the end of inspiration. The results further indicate because of large radial diffusion rates, no significant concentration differences exist between the periphery a-ld the central core of alveolated airways. Simulation of the expired concentration profiles indicate that boll injected very late during inspiration exhibit a sloping tail, unlike the earlier injected boll whose tails are virtually horizontal. Through the use of superposition teehniqnes, it was found that these sloping tails correspond to an alveolar slope of 1.5 vol% between 750 and 1250 ml expired for a continuous washing of tracer. This result is in disagreement with other transport analyses which did not directly account for the effect of geometric asymmetries.

Notes: Abstract Assuming a spherical geometry for the left ventricle, passive elastic stiffness-stress relations have been obtained on the basis of linear elasticity theory and large deformation theory. Employing pressure-volume aata taken from rat hearts of various age groups, it is shown that young rat heart muscle (1 month) is stiffer than either adult (7 months) or old rat heart muscle (17 months). Although the qualitative results are similar for both elasticity theories, the large deformation theory gave results in closer agreement with those obtained from papillary muscle studies. These results imply that stiffness of muscleper se can be assessed from left ventricular pressure-volume data.

Notes: Abstract Deliberate evaluation of the quantum theory of nerve excitation is made by comparing it with Hill's theory in fitting the experimental data on threshold-frequency relation, optimum frequency (v0) for nerve excitation and strength-duration relation. Decrease of v0 and increase of all the time constants (Hill's λ andk, Wei'sT 2 and spike durationw) with decreasing temperature are interpreted on the basis of the dipole relaxation timeT 2 but inexplicable from Hill's theory or any other existing theory. The closeness ofk,T 2 andw values is explained. A variety of experimental results obtained by others is discussed. Finally, a comparison is made between the Hodgkin-Huxley equations and the quantum theory. Most of the facts (electrical and non-electrical) tend to support the thesis that nerve excitation is a macroscopic expression of quantum transitions of dipoles between energy states.

Notes: Abstract In the periodic Leslie model the asymptotic period of total population is a divisor of the asymptotic period of the population vector. Under reasonable circumstances these periods are identical.

Notes: Abstract A number of biological branching systems, such as the bronchial and pulmonary arterial trees, are being investigated in an ongoing study in order to define their physiological properties. The technique involves the description of branching trees by the use of hierarchical systems of ordering, especially those described by Horsfield and by Strahler. During this work some mathematical properties of branching trees were demonstrated and these are described in this paper.

Notes: Abstract Luteinizing hormone (LH) is secreted continuously from the anterior pituitary gland. The concentration in the blood of this gonadotropic hormone plays a regulatory role in the development of puberty in both sexes, in the induction of ovulation in females, and in the production of testosterone in males. The secretion of LH is in turn controlled by luteinizing hormone releasing hormone (LHRH) secreted by the hypothalamus. LH and LHRH are removed from the blood by degradation and excretion. This hormonal system is modelled by a system of ordinary differential equations based upon specific physiological and biochemical assumptions current among experimentalists in this field. The one exception is the assumption that LHRH may bind reversibly to a serum protein; an analysis of the data shows that this or a similar mechanism is a crucial specification. Data on the serum levels of LH and LHRH in two human subjects were fitted using the model. The data consist of the transients and subsequent decays created by a bolus intravenous injection of LHRH.

Notes: Abstract A thick-wall incompressible, elastic sphere was used as a model for the diastolic rat left ventricle. A model for myocardial nonhomogeneity was derived assuming that fiber (circumferential) stress was independent of position in the ventricular wall. The theoretical implications of the resulting constitutive relations together with the spherical model were analyzed in the context of large deformation elasticity theory. It was found that muscle stiffness at a given level of uniaxial stress increased monotonically from the endocardium to the epicardium. In addition, fiber stress was found to be essentially a linear function of transmural pressure above a pressure of 6 g/cm2. It was also shown theoretically that neglecting the nonhomogeneity of the myocardium resulted in a state of stress which differed significantly from that predicted by the nonhomogeneous model. For example, at a transmural pressure of 14 g/cm2, fiber stress in the nonhomogenous model was equal to 17 g/cm2 while fiber stress in the homogeneous model varied between 100 g/cm2 at the endocardial surface and 2 g/cm2 at the epicardial surface. The change in muscle stiffness with position which characterized the nonhomogeneous model also tended to linearize the highly curvilinear radial stress distribution predicted by the homogeneous model at a given transmural pressure.

Notes: Abstract The phenomenon of axonal transport has been well documented (Ochs, 971; Lasek, 1970; and Grafstein, 1967). In a previous paper, we showed how diffusion alone could not account for this process. In this report we show that convection or convection with diffusion can account for the observed build-up of material. By including a first-order catabolic sequestration term, we are able to offer an understanding of the several apparent rates of transport with the same underlying velocity and variable sequestration.

Notes: Abstract It is known that the Lotka-Volterra coupled nonlinear differential equations for a two-species prey-predator ecosystem possess a periodic solution, although its exact form is not yet obtained analytically. The conventional linearization approximation for solving these nonlinear equations leads to a harmonic oscillator whose frequency depends only on the intraspecific coefficients. We propose here a prescription for obtaining nonlinear correction to the linear frequency by using the Hamilton-Jacobi canonical formalism of classical mechanics. It is found that the first-order correction, which also involves interspecific parameters, exhibits the basic qualitative features of the nonlinearity.

Notes: Abstract Environmental safety testing typically requires procedures for extrapolating from the relatively high experimental to the very low use doses of potentially harmful substances. In the present paper, a stochastic mammillary compartmental model for environmental safety testing is proposed and extrapolation procedures based on its dose-response relationship are developed. The proposed model is a direct generalization of one of the basic safety models, the one-hit model, in that a harmful reaction is assumed to occur if at any time any of the peripheral compartments attains a specified threshold of particles. Consideration of a closed model yields an upper bound on the probability of attaining a certain threshold level, thus providing a conservative procedure for extrapolating to a low dose, while a lower bound obtained from a related open model provides a useful monitoring device as to the sharpness of the upper, bound. The extrapolation procedure is illustrated with simulated data and approximations for initial values are developed.

Notes: Abstract By using chromosome images as a framework, algorithms for finding most dissimilar images are presented and illustrated by examples. In terms of angles, a chromosome image consists of two exterior biangles and two interior biangles. Biangles are defined and classified into 180° biangles, 〉180° biangles and 〈180° biangles. The dissimilarity of biangles and its geometric interpretation together with various properties of biangles are also presented. The results may have useful applications in pattern recognition, scene analysis, information storage and retrieval, artificial intelligence and fuzzy set theory.

Notes: Abstract The Volterra equations which represent competitions between two species are utilized to examine the phenomenon of boundary formation between two species of plants. The set of stable stationary points for these equations is determined and is illustrated in a product space of parameters and dynamical variables. The stages of boundary appearance and succession are visualized by considering slow changes of the parameters as functions of time and space.

Notes: Abstract We postulate that the biomass distribution function for an ecological population may be derived from the condition that the biomas diversity functional is maximal subject to an energetic constraint on the total biomass. This leads to a biomass distribution of the form $$p(m) = \bar m^{ - 1} \exp ( - m/\bar m)$$ , where $$\bar m$$ is the mean biomass per individual. The same condition yields a unique value for the biomass diversity functional. These predictions are tested against fishery data and found to be in good agreement. It is argued that the existence of a unique value for biomass diversity may provide a preliminary theoretical foundation for the observed upper limit to species diversity.

Notes: Abstract Pigment distribution presages hydranth regeneration in the marine hydroidTubularia. We suggest that such a distribution could result from a reaction-diffusion system. A model system based on a practical reaction scheme is studied and spatial structures found which closely resemble this pigment distribution. Finite-amplitude spatial structures in reaction-diffusion systems are considered. Whereas in one spatial dimension the final structures are normally very similar to the transient patterns which emerge from a linear analysis, it is shown that in more than one dimension this is not necessarily the case. The reasons for this are discussed.

Notes: Abstract An analytical model is used to described the behavior of inhaled particulate matter in the human respiratory tract. Three different geometries, symmetric and asymmetric, are utilized to simultate the tracheobronchial (TB) tree. The suitability of each geometry for representing the human is evaluated by comparing calculated aerosol deposition probabilities with experimental data from inhalation exposure tests. A symmetric, dichotomously branching pattern is found to be a reliable description of the TB tree for studies of factors affecting aerosol deposition in the human lung. Calculations with the theoretical model are in excellent agreement with measured aerosol deposition efficiencies. Furthermore, the model accurately predicts experimentally observed features of inhalation exposure data, such as effects of inter-subject lung morphology differences and relative efficiencies of specific deposition mechanisms, on aerosol deposition patterns in the TB tree.

Notes: Abstract For certain environments, the Darwinian model allows unique prediction of a function that any surviving system adapted to such an environment has to perform. This is the case for those environments that determine a “survival functional” of position in space-time of known shape. Purely temporal survival functionals can be distinguished from spatial and mixed ones. In each case, there exists an optimum path in combined physical and (reduced) metabolic space. Dependent on the admissible error, approximate solutions of different complexity are sufficient. All solutions possess an afferent, a central, and an efferent part. Within this general frame, specific, “probably simplest”, solutions are proposed for adaptive chemotaxis, insect locomotion, lower vertebrates locomotion, higher vertebrates locomotion, chronobiological systems, and immune systems, respectively—or rather, for the underlying functionals.

Notes: Abstract Mathematical models afford a procedure of unifying concepts and hypotheses by expressing quantitative relationships between observables. The model presented indicates the roles of both insulin and glucagon as regulators of blood glucose, albeit in different ranges of the blood glucose concentrations. Insulin secretion is induced during hyperglycemia, while glucagon secretion results during hypoglycemia. These are demonstrated by simulations of a mathematical model conformed to data from the oral glucose tolerance test and the insulin infusion test in normal control subjects and stable and unstable diabetic patients. The model studies suggest the parameters could prove of value in quantifying the diabetic condition by indicating the degree of instability.

Notes: Abstract A model for the dynamics of a single-species population whose birth rate depends on densities of previous generations is introduced. A difference equation formulation is proposed and the solutions classified for the various parameter values. Data from an experimental population of mice growing in limited space is cited and compared with the model predictions.

Notes: Abstract All soft tissues are modeled as either one-dimensionalstrings, two-dimensionalmembranes, or three-dimensionalsolids. Attention is restricted to tissues in which one of the principal stress components is large and positive in comparison with the other negligible components. Results indicate the following: (1) If a deformed string isconstrained to lie on a surface and is free of tangential pressure, the tension is carried by rays which are geodesics of the surface. If a string or membrane isfree to deform in space without normal pressure, the tension rays are straight lines. If a membrane deforms without tangential surface loads, the tension rays are always geodesics on the deformed surface. If a solid deforms without body forces, the tension rays are straight lines. (2) The stress in a string is a constant if the string is free of tangential pressure and has constant cross-sectional area. The stress in flat tension fields free of tangential surface loads decays inversely with distance along a tension ray from the edge of regression. The stress in a spherically symmetric tension field free of body forces decays inversely with the square of the distance from the center of the sphere. (3) Stress singularities can occur in soft tissues, such as at the corners of a closed rectangular hole in a flat membrane strip. (4) The tension rays in the torsion of soft annular membranes are more steeply inclined from the radial direction than the tension rays for hard metals equally displaced.

Notes: Abstract In this paper we are concerned with problems of the long-term behavior for nonlinear systems in random environment. The general model is assumed to be given by an ordinary differential equation with random parameters or random input. The disturbance process can be taken from a fairly general class of Markov processes having a bounded state space. In terms of the system’s dynamics we give sufficient conditions for the existence and uniqueness of invariant probabilities. Finally, we apply these results to the two-dimensional biochemical model which is known as the Brusselator.

Notes: Abstract In various applications one faces the problem of estimating a signal from discontinuous observations. For example, in biomedical applications the signal may be the ‘state’ of a given organ and one observes through an external counter the amount of radioactivity sequestered by the organ after injection of a radioactive tracer. Here the problem is studied in the context of nonlinear filtering when the signal can be modelled as either a random variable or a diffusion process, and the observations have a continuous and a purely discontinuous component; both components may be affected by the signal. When the signal is a random variable an explicitly computable solution is obtained; for the diffusion case the solution is given as a sequence of approximating filters that can be computed recursively.

Notes: Abstract Eigenvalue problems arise in various biological models. We outline a useful comparison method and a technique using Lyapunov functions that can be applied in many cases. An application to lateral diffusion is discussed.

Notes: Abstract This paper reviews, up to their recent developments, two types of models of the cell cycle: those considering the size controls over the cycle events and the transition probability models. The distribution of inter-mitotic time and the sister-sister and motherdaughter correlations implied by the two approaches are discussed in view of some relevant experimental data.

Notes: Abstract The development of a blood cell line originating from a pluripotent stem cell pool is modelled by a chain of multidimensional branching processes in which the sojourn times of the cells in certain resting states depend on the size of the following subpopulation. The stability of such a model is discussed qualitatively and some considerations concerning a possible malignant degeneration are presented. The behaviour of models for normal and malignant cell production are illustrated by stochastic stimulations. The model presented here describes the development of a certain line of blood cells (e.g. erythrocytes, monocytes or granulocytes) originating from the pluripotent stem cell up to the functional cell in the blood (for related models see, e.g., Rubinow and Lebowitz,J. math. Biol. 1, 87–225;Biophys. J. 16, 897–910).

Notes: Abstract This paper reviews some recent advances in single population stochastic differential equation growth models. They are a natural way to model population growth in a randomly varying environment. The question of which calculus, Itô or Stratonovich, is preferable is addressed. The two calculi coincide when the noise term is linear, if we take into account the differences in the interpretation of the parameters. This clarifies, among other things, the controversy on the theory of niche limiting similarity proposed by May and MacArthur. The effects of correlations in the environmental fluctuations and statistical methods for estimating parameters and for prediction based on a single population trajectory are mentioned. Applications to fisheries, wildlife management and particularly to environmental impact assessment are now becoming possible and are proposed in this paper.

Notes: Abstract A survey is given of the application of (functions of) continuous-time Markov chains in the statistical analysis of behavioural time series.

Notes: Abstract This paper demonstrates that there is one and only one solution to a non-linear singular two-point boundary-value problem which describes oxygen diffusion in a spherical cell. Previous authors have calculated numerical results that differ substantially. Numerical computations using the multiple shooting method support the results of McElwain.

Notes: Abstract The mathematics of distance geometry constitutes the basis of a group of algorithms for revealing the structural consequences of diverse forms of information about a macromolecule's conformation. These algorithms are of proven utility in the analysis of experimental conformational data. This paper presents the basic theorems of distance geometry in Euclidean space and gives formal proofs of the correctness and, where possible, of the complexity of these algorithms. The implications of distance geometry for the energy minimization of macromolecules are also discussed.

Notes: Abstract A fully developed pulsatile flow in a circular rigid tube is analysed by a microcontinuum approach. Solutions for radial variation of axial velocity and cell rotational velocity across the tube are obtained using the momentum integral method. Simplified forms of the solutions are presented for the relevant physiological data. Marked deviations in the results are observed when compared to a Newtonian fluid model. It is interesting to see that there is sufficient reduction in the mass flow rate, phase lag and friction due to the micropolar character of the fluid.

Notes: Abstract A mathematical model of the transport of fluorescein across the blood-retina barrier in the transient state and the subsequent diffusion of fluorescein in the vitreous body is presented. The function of the barrier is lumped in a single parameter—the permeability. The sensitivity of this parameter due to changes in the other parameters of the model is given. This establishes the foundation for the quantitative assessment of the barrier function through vitreous fluorophotometry.

Notes: Abstract The objective of this preliminary study was to develop a new quantitative method of setting the initial insulin infusion patterns in treatment of diabetic patients. The method is based upon the mathematical estimation of the insulin profile required to maintain the glucose level within the normal range after glucose loading in diabetic patients. Using our previously developed equivalent circuit model of glucose kinetics and the reported data of an intravenous glucose tolerance test (IVGTT) in two groups of normal and diabetic patients, two important physiological parameters of the model (the peripheral tissue's insulin resistivity and the hepatic sensitivity to glucose level) were computed for two clinical groups. Then the insulin profile was obtained by computing the plasma insulin concentrations required to keep the total glucose utilization rate of the tissue and the liver in the diabetic group equal to that of the normal group. The simulation result indicated that the computed insulin profile produced a plasma glucose profile which was more closely matched to the normal group's glucose profile than with the case of emulating the normal group's insulin profile in the diabetic group.