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 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 The entropy budget of a white-tailed deer (50kg) on a maintenance diet and a full-feed diet in a standing posture in an open field under clear nocturnal skies with an air temperature of −20°C is investigated based on the energetics given by Moen. Entropy inflow into a white-tailed deer due to infra-red radiation and entropy outflows from a deer due to infra-red radiation, convection, evaporation of water and conduction to ingested food are calculated. Also the entropy production due to metabolic heat production is estimated. Net entropy flow into a deer from its environment becomes negative. On the assumption that a white-tailed deer is in a steady state in entropy, the total entropy production in a deer on a maintenance diet becomes +0.46 J/sec/K. Positiveness of the entropy production shows that the Second Law of Thermodynamics certainly holds in a white-tailed deer. The entropy production per effective radiating surface area of a deer on a maintenance diet is 0.32×10−4 J/cm2/sec/K. On the other hand, the entropy production in a deer on a full-feed diet is 0.59 J/sec/K and that per effective surface area is 0.41×10−4 J/cm2/sec/K. Uptake of 1 g of food produces 22 J/K of entropy within the body of a white-tailed deer. Comparison is made with the results for entropy production in a lizard and in plant leaves.

Notes: Abstract Biological adaptability has been proved to be analysable by means of the Maximum Entropy Formalism (MAXENT) in some cases of non-interacting systems. This formalism is extended to the biomass statistical structures of populations exhibiting internal interactions (i.e. predatorprey effects).

Notes: Abstract The temporal behaviours of the nonlinear substructure of a self-organized compartmental model of calcium metabolism were investigated. The order-two autocatalytic process included in this simple two-dimensional model is compared to some secondary nucleation mechanisms which should take place at the extracellular fluid-bone interface. The model gives rise to complex dynamic behaviours, and multistability properties, involving up to two stable periodic regimes (birhythmicity), were established in different topological configurations. The bifurcations occurring on the boundaries between regions of different qualitative behaviour have been determined. These properties are discussed in relation to the dynamical behaviour of other two-variable models, especially those including the same nonlinearity.

Notes: Abstract A linearized oscillation theorem due to Kulenović, Ladas and Meimaridou (1987,Quart. appl. Math. XLV, 155–164) and an extension of it are applied to obtain the oscillation of solutions of several equations which have appeared in population dynamics. They include the logistic equation with several delays, Nicholson's blowflies model as described by Gurney, Blythe and Nisbet (1980,Nature, Lond. 287, 17–21) and the Lasota-Wazewska model of the red blood cell supply in an animal. We also developed a linearized oscillation result for difference equations and applied it to several equations taken from the biological literature.

Notes: Abstract A theoretical approach to the explanation of the structural design of metabolic pathway is presented. It is based on the hypothesis that due to natural selection during evolution the cellular metabolism of present-day organisms may be characterized by optimal properties. Two cardinal terms enter the theory: (i) the efficiency of a metabolic pathway and (ii) the evolutionary effort for the change of the kinetic parameters of enzymes by mutations of the corresponding genes. For both quantities simple mathematical expressions are proposed. While the efficiency is related to the reaction rates of the enzymes constituting the metabolic pathway, the evolutionary effort is considered to be a monotonically increasing function of the parameter values. By maximizing the efficiency under the constraint of a fixed evolutionary effort the theory allows the calculation of the optimal parameter distribution as the outcome of evolution processes. The methods developed are applied to the following systems: (a) linear reaction sequences with very low affinities of the enzymes towards substrates, (b) linear sequences consisting of saturable enzymatic reactions, (c) branched metabolic pathways consisting of segments of linear chains and (d) glycolysis of erythrocytes. The conclusion is derived that the optimal distribution of kinetic constants depends strongly on the equilibrium constants of the reactions as well as on the total osmolarity of the metabolic intermediates. Without osmotic constraints the evolutionary effort is mainly spent on the enzymes at the beginning of the chain. Using Michaelis-Menten equations the optimal state is characterized by a decrease of the maximal activities of the enzymes towards the end of the chain. These results are modified if osmotic constraints are taken into account. At the investigation of branched pathways the following results were obtained: firstly, if a certain end product may be synthesized along different pathways those which are thermodynamically more unfavourable (e.g. characterized by a small change of free energy) are eliminated in the course of evolution; secondly, if a branched pathway leads to several important end products those reaction segments which are thermodynamically unfavourable are characterized by a higher evolutionary effort. The application of the theory to a realistic model of glycolysis of erythrocytes leads to a correct description of various functionally important properties of the system, such as the ratio between fluxes through different branches and the ATP/ADP ratio, whereas the theory cannot predict the strong separation of time constants observed in the real glycolytic system. It is concluded that the improvement of the predictive power of the theory necessitates the use of more complex functionals for the efficiency which take into account not only the fluxes but also other system properties such as the stability of the pathway or homoeostatic effects.

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 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 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 Linear birth and death processes are used to derive simple expressions for sequential extinction times and gene fixation probabilities in asexual populations.

Notes: Abstract We are here concerned with the functionf which assigns to each pointP of an object the numberf(P) which is the shortest distance fromP to the border. This function appears in various guises in diverse biological studies. The functionf(P) is itself a measure of shape—or more precisely, an infinite set of measures, one for each point (and hence, in view of its geometric definition, usually in a form inconvenient for use). Thus in this paper we sought a reasonable representative of this infinite set of measures, namely themean of the numbersf(P) asP ranges over all points of the entity. Computability studies are developed for various classes of shapes. For example, (1) the mean for a lamina bounded by a polygon circumscribable about a circle of radiusr isr/3; (2) the mean for a domain bounded by a polyhedron circumscribable about a sphere of radiusr isr/4. The transition from pointwise to piecewisef(P), especially in the non-convex case, requires working with inequalities.

Notes: Abstract A mathematical model for traveling bands of motile and chemotactic bacteria in the presence of cell growth and death is examined. It is found that asymptotic traveling wave solutions exist in the absence of chemotaxis, due to the balance of growth, death and random motility. Thus random motility confers the ecological advantage of population propagation through migration into nutrient-rich regions. The presence of chemotaxis amplifies this advantage by moving more cells into higher nutrient concentration regions, resulting in larger and faster bands. Therefore there seem to be two types of traveling bands that can be attained by chemotactic bacteria in the presence of growth and death: (1) these growth/death/motility bands; and (2) pure chemotactic ‘Keller-Segel'-type bands. Comparison to experimental observations by Chapman in 1973 indicate that the latter seem to be formed. The relationship between these two types of solution is at present uncertain. The growth/death/motility bands may have relevance on longer time or distance scales characteristic of microbial ecological systems.

Notes: Abstract Based on the principle of minimum power, a mathematical model of the pathological functional state of the oxygen transport system is presented. The model is used to determine the optimal functional parameters of the oxygen transport system in hyperthyroidism, anemia and hypertension. Theoretical results are compared with clinical data.

Notes: Abstract A modified SIRS model is developed as a framework for the study of epizootiological dynamics in an insect-pathogen system. Linearized stability analysis reveals that the system with one immune and one susceptible host class can exhibit stable, periodic or unstable behavior depending on model parameters. In general, high pathogenicity, short pathogen propagule lifespan and high host reproductive rate are stabilizing influences. Pathogen transmissibility and propagule production/host do not influence local stability. The effect of seasonal host reproduction is studied because most insect hosts are seasonal in temperate climates. The basic stability dependence on model parameters holds except as modified by the length of the reproduction interval. The results of this study are compared with the recent work of Anderson and May.

Notes: Abstract The random predator-prey type interactions of the population species in a random varying environment have been investigated. The Fokker-Planck equation for the transition probability, corresponding to the stochastic differential equations established from Lotka-Volterra equations by the introduction of randomness and variability, has been integrated in the form of a path integral. The transition probabilities for extinction or survival of one or several species have been approximately evaluated and investigated.

Notes: Abstract If two microbial populations compete for a single resource in a homogeneous environment with time invariant inputs they cannot coexist indefinitely if the resource competed for is not renewed by biological activity within the system. Mathematical studies have shown that in a predator-prey system, where the resource (prey) is self-renewing, the two competitors (predators) can coexist in a limit cycle. This suggests that if the resource competed for is renewed by biological activity within the system coexistence can occur in any microbial system provided that it exhibits the same features as, but without being, a predator-prey one. A food chain involving commensalism, competition and amensalism is presented here. Two subcases are considered. It is only when maintenance effects are taken into account that coexistence, in limit cycles, can occur for this system. Limit cycle solutions for the system are demonstrated with the help of computer simulations. Some necessary conditions for coexistence are presented, as are some speculations regarding the possible physical explanations of the results.

Notes: Abstract The branching structure of the mammalian arterial tree has been known to be close to that of an optimal conduit system of the minimum work model characterized as the branch system of constant wall shear rate. The physiological mechanism producing such construction was considered to be based on the local response of arterial caliber induced by the wall shear stress (shear rate × blood viscosity) and thereby maintaining this stress constant, which was previously observed at the canine common carotid artery shunted to the external jugular vein. The stress levels at various parts of the arterial system estimated from available data fell within ±50% of the mean (15 dyn/cm2), which was consistent with the value predicted from the model. Theoretical analyses on the cost function of the model indicated that the suspected variation of shear rate levels in the arterial tree due to the anomalous changes in blood viscosity which might bring about 3- to 4-fold differences between the minimum and maximum shear rates would cause less than 10% increase in the total energy cost. It was concluded that a local adaptive response to wall shear stress is the mechanism which effectively optimizes the design of the arterial tree.

Notes: Abstract The structure of the genetic code is related to a Gray code, which is a plausible theoretical model for an amino acid code. The proposed model implies that the most important factor in shaping the code was the effects of mistakes in translation, not effects of mutations. Another possible implication is that the preservation of stiffness and flexibility at appropriate places in a protein chain is as important in protein structure as the appropriate placement of hydrophilic (external) and hydrophobic (internal) residues. Other results are a simple conceptualization of the relationships among the 20 amino acids and their relations to their codons. The detailed relationships are summarized in the following ‘similarity alphabet’: ala, thr, gly, pro, ser; asp, asn, glu, gln, lys; his, arg, trp, tyr, phe; leu, met, ile, val, cys; (ATGPS DNEQK HRWYF LMIVC in the one-letter code). This alphabet falls into four groups of amino acids: small, external, large, internal. The approximate relation of the groups to their codons is expressed as: the first base of a codon controls size—a purine means a small amino acid, a pyrimidine means large; the middle base controls cloisterednes—purine means external, pyrimidine means internal. These relationships express the minimum change principle upon which the code appears to be founded.

Notes: Abstract A theoretical study of the Brusselator model with non-uniform distribution of component A and a concentration-dependent diffusion coefficient has been performed. Numerical simulation reveals that a variable diffusion coefficient alters the bifurcation pattern and the stability properties of the steady-state as well as periodic solutions. A simple approximate method, based on one-point collocation, has been proposed to analyze the bifurcation phenomena for the case of fixed boundary conditions and low system size.

Notes: Abstract In this paper the effects of the occurrence of cut trees in the topological analysis of branching patterns have been studied. It is assumed that branches are removed at random from the trees. We prove that, for both the segmental and terminal growth models, the probability distributions of the cut trees are identical to those of complete trees.

Notes: Abstract The goal of this work is an examination of capillary exchange models as mathematical operators. The concentration function relations for the Krogh cylinder of a single capillary, basic to many organ models, are studied via the theory of operators on the Lebesgue normed spacesL p[0,∞], (1〈-p〈-∞). A discussion is included of theL p -normsvis-à-vis the coefficient of variation currently used in finding capillary parameters and evaluating parameter searches. The capillary model determines two operators on the space of locally integrable functions: O K (relating extravascular concentration to intravascular) and K a, k (relating intravascular concentration to input), wherek is the ratio of permeabilitysurface area (PS) to extravascular volume, and α is the ratio of PS to flow. These operators are shown to induce contractive (‖O K ‖ p 〈-1, ‖K a, k ‖ p 〈-1), isotone, linear operators onL p . The uniform convergence relation $$K_{a,k} = \mathop {\lim _{(p)} }\limits_{N \to \infty } \left( {\sum\limits_{n = 0}^N {P_n (a)O_k^n } } \right)$$ (as operators onL p) is derived, whereP n (a) is the Poisson probabilitye −a a n /n!. For the important special cases ofp=∞, 1, 2 the norms are found (‖Ok‖=‖Ka,k‖p=1). Consideration is also given to the norms and operators when the functions involved are limited to a finite interval of time.

Notes: Abstract One particular kind of structure offers possible explanations, for long-term memory, efficient consolidation of stored information from the environment, clustering of data strings and multimodal functioning. It is a possible model for pieces of neural structure and its use offers a uniform method for both studying and constructing an extensive class of mechanisms.

Notes: Abstract Levenshtein dissimilarity measures are used to compare sequences in application areas including coding theory, computer science and macromolecular biology. In general, they measure sequence dissimilarity by the length of a shortest weighted sequence of insertions, deletions and substitutions required, to transform one sequence into another. Those Levenshtein dissimilarity measures based on insertions and deletions are analyzed by a model involving valuations on a partially ordered set. The model reveals structural relationships among poset, valuation and dissimilarity measure. As a consequence, certain Levenshtein dissimilarity measures are shown to be metrics characterized by betweenness properties and computable in terms of well-known measures of sequence similarity.

Notes: Abstract It is pointed out that the methane flux measured experimentally for certain ponds and swamps is quantitatively consistent with a commensal dependence of Methanobacteria on O2-chemotactic motile aerobic bacteria. The Methano species is thereby shielded from oxygen and provided with carbon dioxide for the anaerobic production of methane.

Notes: Abstract A sufficient condition is given for stochastic boundedness persistence of a top predator in generalized Lotka-Volterra-type stochastic food web models in arbitrary bounded regions of state space. The main result indicates that persistence in the corresponding deterministic system is preserved in the stochastic system if the intensities of the random fluctuations are not too large.

Notes: Abstract One of Bobisud's models for the evolution of cannibalism is reanalyzed by applying the method of finding evolutionarily stable strategies (or ESS's). It is demonstrated that ‘no cannibalism’ never will be an ESS if the initial rate of cannibalism is too large. It is further demonstrated that individual selection may even result in the evolution of cannibalism during food abundance. Some empirical case studies are briefly discussed in relation to this model.

Notes: Abstract A computer algorithm is presented which equiprobably generates any member of the set of all directed trees withk labeled terminal nodes and unlabeled interior nodes. The algorithm requires roughlyk 2 /2 storage locations. The one-time initialization requiresO(k 2 ) time, while generating each tree requiresO(k) time.

Notes: Abstract The comparison of several sequences is central to many problems of molecular biology. Finding consensus patterns that define genetic control regions or that determine structural or functional themes are examples of these problems. Previously proposed methods, such as dynamic programming, are not adequate for solving problems of realistic size. This paper gives a new and practical solution for finding unknown patterns that occur imperfectly above a preset frequency. Algorithms for finding the patterns are given as well as estimates of statistical significance.

Notes: Abstract A new development is introduced here in the use of dynamic programming in finding pattern similarities in genetic sequences, as was first done by Needleman and Wunsch (1969). A condition of pattern similarity is defined and an algorithm is given which scans any set of similarities and screens out those which fail to meet the condition. When the set to be scanned contains every pair of segments, one from each of two given sequences of lengthsm andn (i.e. every possible location for a pattern similarity), then it completes the scan in a number of computational steps proportional tom·n, leaving those pairs of segments which satisfy the similarity condition. The algorithm is based on the concept of match density, as suggested by Goad and Kanehisa (1982).

Notes: Abstract This paper concerns sequences of letters in which certain “distinguished” words are of interest. Such sequences arise as data in numerous fields including genetics and neuroscience. A probability distribution is given for the number of occurrences of a chosen word in a randomized sequence of letters. Such words are considered “favored” if they occur more than expected at random. Favored words have been discovered in nerve impulse trains and may reflect a neural coding scheme.

Notes: Abstract As the volume of protein sequence data grows, rapid methods for searching the protein sequence database become of primary importance. Rigorous comparison of sequences is obtained with the well-known dynamic programming algorithms. However, these algorithms are not rapid enough to use for routinely searching the entire database. In this paper we discuss some methods that can be used for rapid searches.

Notes: Abstract We give algorithms for computing the extent of similarity between two or three sequences of letters. The similarity measures we consider include a penalty for inserting gaps within the sequence in order to enhance similarity. The magnitude of the penalty for gaps is assumed to be independent of their size in order to accommodate certain biological applications. Our algorithm for three sequence comparisons, which is based on solving a system of recursive equations, improves upon the efficiency of existing methods. Although the system of recursive equations utilized by the algorithm is quite complicated as it stands, it has none the less been simplified by appeal to combinatorial considerations.

Notes: Abstract Well-known dynamic programming algorithms exist for comparing two finite sequences inO(N 2) time and storage, whereN is the common sequence length. Extensions to the comparison ofM finite sequences requireO((2N) M) time and storage, making such algorithms difficult even forM=3. A simple generalization of the sequences makes it possible to obtain some results about the geometry of sequence alignments. These ideas suggest heuristic approaches to problems of comparing several sequences. IfM sequences are known to be related by a binary tree, they can be aligned inO(MN 2) time andO(N 2+NM) storage.

Notes: Abstract Protein sequences of the Dayhoff databank of 1984 have been analyzed to evaluate the occurrences of the 400 dipeptides and 8000 tripeptides. Expected values and standard deviations for the di- and tripeptides were determined by Monte Carlo and binomial approximation. A condensed format containing this information, labeled a uniqueness diagram, is presented and made available in the form of a microfiche.

Notes: Abstract In two-state sliding filament models of muscle contraction a partial differential equation must be solved to find the cross-bridge distribution functionn(x, t). In this paper the analytical form of this function is obtained by integration along the characteristic line and special cases are presented in which the explicit expression forn(x, t) can be completely determined. These analytical solutions provide a direct mathematical connection between the microscopic contraction parameters contained in the kinetic theories and macroscopic muscle dynamics and are thus used to investigate what parameters influence the transient contractile tension in typical experimental conditions. The results of this investigation are consistent with relevant aspects of muscle physiology.

Notes: Abstract Experiments on random binary, ternary, etc. (P=2, 3,…, 10) switching nets are reported. Behavioral cycle lengths are examined as functions of output variety,P, input connectance,K, and net size,N. Overall, output variety appears an influential, well-behaved net property. Strong, but well-behaved interactions appear among net variables. In high connectance nets, median cycle length grows approx. asP N/2. Other factors constant, one-connected nets show the shortest cycles, and connectance effects appear to converge asymptotically aroundN. Data for cycle length as a function of net size suggest a concavity not compatible with the Kauffman “square root law” (Kauffman, 1969). Evidence of a positive relationship between cycle length and run-in length is found in two-input nets; weaker evidence is obtained that in higher connectance nets this relationship becomes negative in sign. The “modular complexity” ofP〉2 nets is examined briefly.

Notes: Abstract The influence of topographical situation on the spread of infection is studied. The investigation is based on a multigroup model. The population under consideration is thought to be divided into subpopulations living in regions that are separated from each other by natural barriers (mountains). Infection is carried from one region to another by migrating infectives. Migration is possible only along the river system so that the structure of the epidemiological network is that of a symmetric tree. The results allow comparison of the velocity of propagation of the epidemic for different geographical situations and allow quantification of the “channel-effect”, according to which mountainous regions are channels rather than barriers to the spread of an epidemic.

Notes: Abstract Bayesian image processing formalisms which incorporatea priori information about valued-uncorrelated and valued-correlated (patterned) source distributions are introduced and the corresponding iterative algorithms are derived using the EM technique. Striking improvement in image processing is demonstrated when applying these algorithms to Poisson and Gaussian randomized data in one-dimensional cases.

Notes: Abstract In the framework of the neural network theory effects similar to hypnotic displays are constructed. They are based on the associative paradigm involving non-linear interaction of excitatory and inhibitory channels with synaptic memory. The non-linearity of long-term memorizing processes may cause effects exhibited by blind spots, which are interpreted as the first stage of hypnosis. More complicated phenomena are discussed in terms of a two-layer network.

Notes: Abstract Mutation is introduced into autocatalytic reaction networks. The differential equations obtained are neither of repliator-type nor can they be transformed straightway into a linear equation. Examples of low dimensional dynamical systems —n=2, 3 and 4 — are discussed and complete qualitative analysis is presented. Error thresholds known from simple replication-mutation kinetics with frequency independent replication rates occur here as well. Instead of cooperative transitions or higher order phase transitions the thresholds appear here as supercritical or subcritical bifurcations being analogous to first-order phase transitions.

Notes: Abstract The non-linear behavior of a differential equations-based predator-prey model, incorporating a spatial refuge protecting a consant proportion of prey and with temperature-dependent parameters chosen appropriately for a mite interaction on fruit trees, is examined using the numerical bifurcation code AUTO 86. The most significant result of this analysis is the existence of a temperature interval in which increasing the amount of refuge dynamically destabilizes the system; and on part of this interval the interaction is less likely to persist in that predator and prey minimum population densities are lower than when no refuge is available. It is also shown that increasing the amount of refuge can lead to population outbreaks due to the presence of multiple stable states. The ecological implications of a refuge are discussed with respect to the biological control of mite pests.

Notes: Abstract In many applications of control theory on plant growth models biomass maximization is postulated to avoid analytically unsolvable problems while fruit maximization is commonly considered to be a more realistic criterion. In a special case, we are able to compare these criteria. Iwasa and Roughgarden (1984,Theor. Pop. Biol. 25, 78–105) have investigated a certain class of plant growth models using a fruit maximization criterion. They proved that, in the vegetative growth period, the organs follow a certain path of balanced growth. We show that this path remains optimal when biomass maximization is postulated. This underlines the importance of the balanced growth path found by Iwasa and Roughgarden. Furthermore, our result suggests that in the vegetative growth period the biomass maximization criterion is a good approximation of fruit maximization. In another theoretical control investigation, Schultzeet al. (1983,Oecologia 58, 169–177) derived a different type of balanced growth path. We apply the theory of Iwasa and Roughgarden to an improved version of the model of Schulzeet al. This leads to a new description of balanced growth between root and shoot that reflects non-linearities in the water uptake process and constitutes an interesting hypothesis for further experimental testing.

Notes: Abstract In this paper the effects of changing the ion concentration in and around a sample of soft tissue are investigated. The triphasic theory developed by Laiet al. (1990,Biomechanics of Diarthrodial Joints, Vol. 1, Berlin, Springer-Verlag) is reduced to two coupled partial differential equations involving fluid ion concentration and tissue solid deformation. These equations are given in general form for Cartesian, cylindrical and spherical geometries. After solving the two equations quantities such as fluid velocity, fluid pressure, chemical potentials and chemical expansion stress may be easily calculated. In the Cartesian geometry comparison is made with the experimental and theoretical work of Myerset al. (1984,ASME J. biomech. Engng,106, 151–158). This dealt with changing the ion concentration of a salt shower on a strip of bovine articular cartilage. Results were obtained in both free swelling and isometric tension states, using an empirical formula to acount for ion induced deformation. The present theory predicts lower ion concentrations inside the tissue than this earlier work. A spherical sample of tissue subjected to a change in salt bath ion concentration is also considered. Numerical results are obtained for both hypertonic and hypotonic bathing solutions. Of particular interest is the finding that tissue may contract internally before reaching a final swollen equilibrium state or swell internally before finally contracting. By considering the relative magnitude, and also variation throughout the time course of terms in the governing equations, an even simpler system is deduced. As well as being linear the concentration equation in the new system is uncoupled. Results obtained from the linear system compare well with those from the spherical section. Thus, biological swelling situations may be modelled by a simple system of equations with the possibility, of approximate analytic solutions in certain cases.

Notes: Abstract Many models of immune networks have been proposed since the original work of Jerne [1974,Ann. Immun. (Inst. Pasteur) 125C, 373–389]. Recently, a limited class of models (Weisbuchet al., 1990,J. theor. Biol. 146, 483–499) have been shown to maintain immunological memory by idiotypic network interactions. We examine generalizations of these models when the networks are both large and highly connected to study their memory capacity, i.e. their ability to account for immunization to a large number of random antigens. Our calculations show that in these minimal models, random connectivities with continuously distributed affinities reduce the memory capacity to essentially nil.

Notes: Abstract A kinetic model is proposed to delineate the factors that determine the coronary reactive hyperemic response (RHR) to transient ischemia. The model comprises of myocardial-interstitial (M) and vascular (V) compartments. Vasodilator metabolites (VM) are produced in the M compartment during the interval of coronary occlusion. The rate of VM production is dependent on the flow rate during the ischemic period, the ratio of excess flow above the control level (R) to the loss of flow during occlusion period (D), the amount of oxygen stored and the degree of vasodilation in the V compartment prior to occlusion. Following a complete release of occlusion, VM are transported from the M to V compartment and are washed out or degraded with time. The time course of RHR is determined by the coronary patency which is proportional to VM concentration in the V compartment. Based on a set of numerical constants, the model is tested by simulating RHR to the various occlusion manoeuvres: a pair of 10 sec occlusions separated by brief release, a 15 sec release followed by a second brief occlusion, a brief release of an occlusion followed by restriced inflow and a period of restricted inflow after occlusion. The simulated results fit the experimental R/D and RH durations data of canine hearts. Factors that determine the impairment of RH capacity in coronary stenosis are suggested in terms of the model scheme.

Notes: Abstract In the present paper a kinetic study is made of the behaviour of a Michaelis-Menten enzyme-catalysed reaction in the presence of irreversible inhibitors rendered unstable in the medium by their reaction with the product of enzymatic catalysis. A general mechanism involving competitive, non-competitive, uncompetitive and mixed irreversible inhibition with one or two steps has been analysed. The differential equation that describes the kinetics of the reaction is non-linear and computer simulations of its dynamic behaviour are presented. The results obtained show that the systems studied here present kinetic co-operativity for a target enzyme that follows the simple Michaelis-Menten mechanism in its action on the substrate, except in the case of an uncompetitive-type inhibitor.

Notes: Abstract The relative contributions of mitochondrial β-oxidation and peroxisomal β-oxidation and peroxisomal ω-oxidation to the oxidation of a given fatty acidin vivo can be quantitated by an isotopic method. The approach requires infusion of a fatty acid labelled on two specific carbon atoms (e.g. [1-14C] and [11-14C] palmitate) to an isotopic steady state, with subsequent isolation and degradation of an acetylated conjugate as a product of the liver cytosolic acetyl CoA pool and of ketone bodies as a product of the liver mitochondrial acetyl CoA pool.

Notes: Abstract Pancreatic β-cells in intact islets of Langerhans perfused with various glucose concentrations exhibit periodic bursting electrical activity (BEA) consisting of active and silent phases. The fraction of the time spent in the active phase is called the plateau fraction and appears to be strongly correlated with the rate of release of insulin from islets as glucose concentration is varied. Here this correlation is quantified and a theoretical development is presented in detail. Experimental rates of insulin release are correlated with “effective” plateau fractions over a range of glucose concentrations. There are a number of different models for BEA in pancreatic β-cells and a method is developed here to quantify the dependence of a glucose dependent parameter on glucose concentration. As an example, the plateau fractions computed from the Sherman-Rinzel-Keizer model are matched with experimental plateau fractions to obtain a relationship between the model's glucose-dependent parameter, β, and glucose concentration. Knowledge of the relationships between β and glucose concentration and between experimental measurements of rates of insulin release and plateau fractions permits the determination of theoretical rates of insulin release from the model.

Notes: Abstract When a suspension of bacterial cells of the speciesBacillus subtilis is placed in a chamber with its upper surface open to the atmosphere complex bioconvection patterns are observed. These arise because the cells: (1) are denser than water; and (2) usually swim upwards, so that the density of an initially uniform suspension becomes greater at the top than the bottom. When the vertical density gradient becomes large enough, an overturning instability occurs which ultimately evolves into the observed patterns. The reason that the cells swim upwards is that they are aerotactic, i.e. they swim up gradients of oxygen, and they consume oxygen. These properties are incorporated in conservation equations for the cell (N) and oxygen (C) concentrations, and these are solved in the pre-instability phase of development whenN andC depend only on the vertical coordinate and time. Numerical results are obtained for both shallow- and deep-layer chambers, which are intrinsically different and require different mathematical and numerical treatments. It is found that, for both shallow and deep chambers, a thin boundary layer, densely packed with cells, forms near the surface. Beneath this layer the suspension becomes severely depleted of cells. Furthermore, in the deep chamber cases, a discontinuity in the cell concentration arises between this cell-depleted region and a cell-rich region further below, where no significant oxygen concentration gradients develop before the oxygen is fully consumed. The results obtained from the model are in good qualitative agreement with the experimental observations.

Notes: Abstract We describe a classification scheme for bursting oscillations which encompasses many of those found in the literature on bursting in excitable media. This is an extension of the scheme of Rinzel (inMathematical Topics in Population Biology, Springer, Berlin, 1987), put in the context of a sequence of horizontal cuts through a two-parameter bifurcation diagram. We use this to describe the phenomenological character of different types of bursting, addressing the issue of how well the bursting can be characterized given the limited amount of information often available in experimental settings.

Notes: Abstract To ensure its sustained growth, a tumour may secrete chemical compounds which cause neighbouring capillaries to form sprouts which then migrate towards it, furnishing the tumour with an increased supply of nutrients. In this paper a mathematical model is presented which describes the migration of capillary sprouts in response to a chemoattractant field set up by a tumour-released angiogenic factor, sometimes termed a tumour angiogenesis factor (TAF). The resulting model admits travelling wave solutions which correspond either to successful neovascularization of the tumour or failure of the tumour to secure a vascular network, and which exhibit many of the characteristic features of angiogenesis. For example, the increasing speed of the vascular front, and the evolution of an increasingly developed vascular network behind the leading capillary tip front (the brush-border effect) are both discernible from the numerical simulations. Through the development and analysis of a simplified caricature model, valuable insight is gained into how the balance between chemotaxis, tip proliferation and tip death affects the tumour's ability to induce a vascular response from neighbouring blood vessels. In particular, it is possible to define the success of angiogenesis in terms of known parameters, thereby providing a potential framework for assessing the viability of tumour neovascularization in terms of measurable quantities.

Notes: Abstract The paper presents a simple method for solving the optimal (LQG) control problem on an infinite time horizon for linear plants described by proper rational transfer function matrices. Since the class of proper plants discussed in the paper is more general than the one commonly used, additional properties of solutions are presented. Two numerical examples illustrate the theoretical considerations. The examples have been performed on an IBM PC using the proposed algorithm. This algorithm is a part of a public ASMCS package developed by the author for analysis, synthesis, and simulations of multi-input, multi-output (MIMO) control systems.

Notes: Abstract We study the behavior of a Hilbert network (i.e., a finite or countably infinite network whose variables are in a Hilbert space and in which the associated total energy is finite) whose elements are affected by perturbations. More specifically, we will give estimates for a change of the current distribution caused by (a) perturbations of the elements of the nominal network when the voltage sources are fixed, and (b) a change of voltage sources in a network whose elements are perturbed. The conditions given in our theorems imply insensitivity and robust stability of the nominal network. The applications of the results are illustrated by an example of an infinite network.

Notes: Abstract A new lattice filter structure to model two-dimensional (2-D) autoregressive (AR) fields is proposed. The proposed structure utilizes and extracts the information contained in the backward prediction error fields and their delayed versions. The main idea is to use two sets of reflection coefficients corresponding to two quadrant filters and to increase the number of reflection coefficients with the order of the lattice filter. Increasing the number of reflection coefficients at each stage produces a sufficient number of independent parameters to model AR fields up to order three, which is an improvement over the existing 2-D lattice filter structures. The improvement is confirmed by computer simulations. In addition, a relationship between the reflection coefficients and the AR coefficients is derived. It is also shown that the entropy contained in the backward prediction error field vector of the proposed structure is closer to the input entropy when compared to those contained in existing 2-D lattice filters.

Notes: Abstract This paper addresses the problem of global asymptotic stability in recursive first hyperquadrant causalm-D digital filters. A set of simple-to-check 1-D conditions necessary for stability of them-D system is given. Generalizations tov-dimensional (v〈m) subsystems are also provided. These derived conditions are useful in determining asymptotic convergence ofm-D digital filters implemented in fixed or floating point arithmetic. The set of 1-D conditions can be considered the analogous result to practical bounded input bounded output (BIBO) stability for linearm-D systems.

Notes: Abstract A new efficient algorithm for calculating the weighting coefficients of maximally linear, FIR digital differentiators is presented. Simple closed-form explicit and recursive formulas are derived in a very straightforward manner. Moreover, a simple recursive equation is established, relating coefficients of two digital differentiators of adjacent ranks.