ALBERT

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 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 The parameter domain for which the quasi-steady state assumption is valid can be considerably extended merely by a simple change of variable. This is demonstrated for a variety of biologically significant examples taken from enzyme kinetics, immunology and ecology.

Notes: Abstract The ability of random fluctuations in selection to maintain genetic diversity is greatly increased when generations overlap. This result has been derived previously using genetic models with very special assumptions about the population age structure. Here we explore its robustness in more realistic population models, with very general age structure or physiological structure. For a range of genetic models (haploid, diploid, single and multilocus) we find that the condition for maintaining genetic diversity generalizes almost without change. Genetic diversity is maintained by selection if a product of the form (generation overlap)×(selection intensity)×(variability in the selection regime) is sufficiently large, where the generation overlap is measured in units of Fisher's reproductive value. This conclusion is based on a local evolutionary stability analysis, which differs from the standard “protected polymorphism” criterion for the maintenance of genetic diversity. Simulation results match the predictions from the local stability analysis, but not those from the protected polymorphism criterion. The condition obtained here for maintaining genetic diversity requires fitness fluctuations that are substantial but well within the range observed in many studies of natural populations.

Notes: Abstract Premised on relatively simple assumptions, mathematical models like those of Monod, Pirt or Droop inadequately explain the complex transient behavior of microbial populations. In particular, these models fail to explain many aspects of the dynamics of aTetrahymena pyriformis-Escherichia coli community. In this study an alternative approach, an individual-based model, is employed to investigate the growth and interactions ofTetrahymena pyriformis andE. coli in a batch culture. Due to improved representation of physiological processes, the model provides a better agreement with experimental data of bacterial density and ciliate biomass than previous modeling studies. It predicts a much larger coexistence domain than rudimentary models, dependence of biomass dynamics on initial conditions (bacteria to ciliate biomasses ratio) and appropriate timing of minimal bacteria density. Moreover, it is found that accumulation ofE. coli sized particles andE. coli toxic metabolites has a stabilizing effect on the system.

Notes: Abstract At the core of contemporarymorphometrics—the quantitative study of biological shape variation—is a synthesis of two originally divergent methodological styles. One contributory tradition is the multivariate analysis of covariance matrices originally developed as biometrics and now dominant across a broad expanse of applied statistics. This approach, couched solely in the linear geometry of covariance structures, ignores biomathematical aspects of the original measurements. The other tributary emphasizes the direct visualization of changes in biological form. However, making objective the biological meaning of the features seen in those diagrams was always problematical; also, the representation of variation, as distinct from pairwise difference, proved infeasible. To combine these two variants of biomathematical modeling into a valid praxis for quantitative studies of biological shape was a goal earnestly sought though most of this century. That goal was finally achieved in the 1980s when techniques from mathematical statistics, multivariate biometrics, non-Euclidean geometry and computer graphics were combined in a coherent new system of tools for the complete regionalized quantitative analysis oflandmark points together with the biomedical images in which they are seen. In this morphometric synthesis, correspondence of landmarks (biologically labeled geometric points, like “bridge of the nose”) across specimens is taken as a biomathematical primitive. The shapes of configurations of landmarks are defined as equivalence classes with respect to the Euclidean similarity group and then represented as single points in David Kendall'sshape space, a Riemannian manifold with Procrustes distance as metric. All conventional multivariate strategies carry over to the study of shape variation and covariation when shapes are interpreted in the tangent space to the shape manifold at an average shape. For biomathematical interpretation of such analyses, one needs a basis for the tangent space compatible with the reality of local biotheoretical processes and explanations at many different geometric scales, and one needs graphics for visualizing average shape differences and other statistical contrasts there. Both of these needs are managed by thethin-plate spline, a deformation function that has an unusually helpful linear algebra. The spline also links the biometrics of landmarks to deformation analysis of the images from which the landmarks originally arose. This article reviews the history and principal tools of this synthesis in their biomathematical and biometrical context and demonstrates their usefulness in a study of focal neuroanatomical anomalies in schizophrenia.

Notes: Abstract A competition model describing tumor-normal cell interaction with the added effects of periodically pulsed chemotherapy is discussed. The model describes parameter conditions needed to prevent relapse following attempts to remove the tumor or tumor metastasis. The effects of resistant tumor subpopulations are also investigated and recurrence prevention strategies are explored.

Notes: Abstract Increasing attention is being paid to the configuration and development of vascular structures and their possible correlations with physiological events. The study of angiogenesis in normal and pathological states as well as in the embryo and adult has provided new insights into the mechanism of vessel growth and organization of the vasculature. Various mathematical branching models have been developed. These constructions are mainly geometrical and only involve a branching phenomenon. We propose the use of a deterministic non-linear model based on physiological laws and hydrodynamics. Growth, branching and anastomosis, the three actual main events occurring in vascular growth, are included in this model. Space growth, including cells and vessels, is defined by a decreasing transformation. Space density and the length of new sprouts are controlled by a set of parameters. The conditions on these parameters are well established, which allows the production of realistic patterns.

Notes: Abstract The quasi-stationary distribution of a population within a system of interacting populations is approximated by a stochastic logistic process. The parameters of this process can be expressed in the parameters of the full system. Using the diffusion approximation, an expression for the expected extinction time is derived from this logistic process. Since the expected extinction time is expressed in the parameters of the full system, the effect of these parameters on the extinction risk can be easily evaluated, which may be of use for studies in ecology, conservation biology and epidemiology. The outcome is compared with simulation results for the case of a prey-predator system.

Notes: Abstract The cytokines are the information superhighway of the immune system. They are an important component of the integrated behavior of the system. In order to be able to have a good understanding of the immune system, we must be able to model the effect of cytokines and their combined effect. This work is a step in that direction. We study the combined effect of two cytokines: interleukin-2 (IL-2) and interleukin-4 (IL-4) on some cells of the immune system. Interleukin-2 and interleukin-4 are important growth and differentiation factors for B and T cells. Interleukin-4 antagonizes the effect of interleukin-2 on B cells and some T cells while it synergizes with interleukin-2 on other T cells. We build a mathematical model of the interaction of both cytokines on T and B cells as a building block toward a model of the Th1/Th2 cross-regulation. The response of a given cell to the combination of interleukin-2 and interleukin-4 is shown to involve competing dynamical effects which can lead to either antagnostic or synergistic combined effect.

Notes: Abstract We propose a mathematical approach to the modelling of self-organizing hierarchies in animal societies. This approach relies on a basic positive feedback mechanism that reinforces the ability of a given individual to win or to lose in a hierarchical interaction, depending on how many times it won or lost in previous interactions. Motivated by experiments carried out on primitively eusocial waspsPolistes, the model, is based on coupled differential equations supplemented with a small stochastic term. Numerical integrations allow many different hierarchical profiles to be obtained depending on the model parameters: (1) the particular form of the probability for an individual to win or lose a fight given its history, (2) the probability of interaction between two individuals, (3) the forgetting strength, which determines the rate at which events in the past are forgotten and no longer influence the force of an individual and (4) two individual recognition parameters, which set the contribution of individual recognition in the process of hierarchical genesis. We compare the results, expressed in terms of a hierarchical index or of the Landau number that describes the degree of linearity of the hierarchy, with various experimental results.

Notes: Abstract The normal process of dermal wound healing fails in some cases, due to fibro-proliferative disorders such as keloid and hypertrophic scars. These types of abnormal healing may be regarded as pathologically excessive responses to wounding in terms of fibroblastic cell profiles and their inflammatory growth-factor mediators. Biologically, these conditions are poorly understood and current medical treatments are thus unreliable. In this paper, the authors apply an existing deterministic mathematical model for fibroplasia and wound contraction in adult mammalian dermis (Olsenet al., J. theor. Biol. 177, 113–128, 1995) to investigate key clinical problems concerning these healing disorders. A caricature model is proposed which retains the fundamental cellular and chemical components of the full model, in order to analyse the spatiotemporal dynamics of the initiation, progression, cessation and regression of fibro-contractive diseases in relation to normal healing. This model accounts for fibroblastic cell migration, proliferation and death and growth-factor diffusion, production by cells and tissue removal/decay. Explicit results are obtained in terms of the model processes and parameters. The rate of cellular production of the chemical is shown to be critical to the development of a stable pathological state. Further, cessation and/or regression of the disease depend on appropriate spatiotemporally varying forms for this production rate, which can be understood in terms of the bistability of the normal dermal and pathological steady states—a central property of the model, which is evident from stability and bifurcation analyses. The work predicts novel, biologically realistic and testable pathogenic and control mechanisms, the understanding of which will lead toward more effective strategies for clinical therapy of fibro-proliferative disorders.

Notes: Abstract Populations often exhibit abrupt changes in abundance associated with a smooth, continuous change in some component of their environment, with the abruptness usually attributed to inter-specific interactions or physical extremes. This paper presents a spatially explicit single-species population model in which intra-specific interactions alone are responsible for such an abrupt change. The essential mechanism involves cooperation in both colonization (through enhanced recruitment near other individuals) and mortality (protection through a “safety-in-numbers” interaction). Large fluctuations in population density would likely be observable near the transition region.

Notes: Abstract A mathematical model of the nitrogen transformation cycle in an aquatic environment is studied. Using Pontryagin's maximum principle, a preferential utilization of ammonium to nitrate by phytoplankton is explained and verified by experimental data. A multiparameter bifurcation is given. The model was found to have four types of equilibrium sets. It is shown that a Hopf bifurcation may occur.

Notes: Abstract Parallel computation employing a domain decomposition method was used to calculate precisely without approximations the spatio-temporal distribution of Ca2+ in nerve terminals. The results showed, contrary to expectations, that for equal admitted Ca2+ currents at low (one channel open) and high (four channels open) depolarization, the average Ca2+ concentration at the release area is higher at the low depolarization. These calculations provide additional support for the Ca2+-voltage hypothesis for neurotransmitter release.

Notes: Abstract Type I hypersensitivity, which functions to protect the organism from parasites, is caused by binding of antigen to IgE antibodies pre-attached to the cell surface of tissue mast cells and their circulating counterparts, the basophils. In “allergy,” type I hypersensitivity is inappropriately induced by protein-based foreign substances (such as pollen) or protein components of insect stings, which in the normal course of events would be cleared from the organism without causing any damage. Paradoxically, a successful clinical treatment of allergy involves repeated immunization of allergic persons with low doses of the allergen—immunotherapy. Investigation of the available experimental evidence leads to the conclusion that the phenomena of immunotherapy are best addressed in terms of the interplay among the mechanism(s) of immune memory—Th1/Th2 cross-regulation—and the physical compart-mentalization of the immune system. These conclusions are illustrated with a numerical simulation.

Notes: Abstract In the presence of seasonal forcing, predator-prey models with quadratic interaction terms and weak dissipation can exhibit infinite numbers of coexisting periodic attractors corresponding to cycles of different magnitude and frequency. These motions are best understood with reference to the conservative case, for which the degree of dissipation is, by definition, zero. Here one observes the familiar mix of “regular” (neutrally stable orbits and tori) and chaotic motion typical of non-integrable Hamiltonian systems. Perturbing away from the conservative limit, the chaos becomes transitory. In addition, the invariant tori are destroyed and the neutrally stable periodic orbits becomes stable limit cycles, the basins of attraction of which are intertwined in a complicated fashion. As a result, stochastic perturbations can bounce the system from one basin to another with consequent changes in system behavior. Biologically, weak dissipation corresponds to the case in which predators are able to regulate the density of their prey well below carrying capacity.

Notes: Abstract The standard method for measuringin vitro antibiotic efficacy is based on a point observation of bacterial activity 18 hours after inoculation. The method, while simple, forgoes significant information by ignoring the dynamics of the interations between antibiotic and bacteria. This paper proposes a simple dynamic model describing these interactions. The model consists of two non-linear differential equations of the S-system type. Its parameter values are estimated, through the minimization of residual errors, from data on the effect of the carbapenem antibiotic imipenem onPseudomonas aeruginosa. The model adequately describes the dynamic behavior of the bacterial populations in the presence of the antibiotic: beginning with drug administration, then through the decline of the bacterial population and possibly ending with bacterial resurgence.

Notes: Abstract We have studied an ecological system of two species, which we denotestrong andweak, respectively, that compete for a single food resource. This system is modelled as a three component reaction-diffusion process. In the presence of a solitary pulse of increased resources, the weaker competitor can diffuse toward this surplus, gaining a competitive advantage and hence persisting in contraposition with the classical Lotka-Volterra result. An exact analytical solution has been found through a quantum mechanical analogy. A stability analysis of this solution against changes in different parameters has been carried out.

Notes: Abstract Collapsible-tube flow with self-excited oscillations has been extensively investigated. Though physiologically relevant, forced oscillation coupled with self-excited oscillation has received little attention in this context. Based on an ODE model of collapsible-tube flow, the present study applies modern dynamics methods to investigate numerically the responses of forced oscillation to a limit-cycle oscillation which has topological characteristics discovered in previous unforced experiments. A devil's staircase and period-doubling cascades are presented with forcing frequency and amplitude as control parameters. In both cases, details are provided in a bifurcation diagram. Poincaré sections, a frequency spectrum and the largest Lyapunov exponents verify the existence of chaos in some circumstances. The thin fractal structure found in the strange attractors is believed to be a result of high damping and low stiffness in such systems.

Notes: Abstract In this work, we show that a one-dimensional model of the blood flow across the lungs can reproduce the evolution of a bolus versus the time. Solving the differential equation governing the bolus concentration in the framework of this model, we determine the solution which fulfills Gaussian initial boundary conditions. An effective parameter related to the ratio of a diffusion coefficient to the square of the mean speed of the flow is defined. The determination of its numerical values following a semi-empirical approach enables us to know accurately the mean transit time and the cardiac output. The results have been compared to other methods, and were found in good agreement. Such an approach could be of interest in all studies where the knowledge of flow—including micro-circulation—is needed.

Notes: Abstract How two species interact during and after colonization influences which of them will be present in each stage of succession. In the tolerance model of ecological succession in a patchy environment, empty patches can be colonized by any species, but the ability to tolerate reduced resource levels determines which species will exclude the other. Here, we analyze a meta-population model of the possible roles of competition in colonization and succession, using non-linear Markov chains as a mathematical framework. Different kinds of competition affect the final equilibrial, abundances of the species involved in qualitatively different ways. An explicit criterion is given to determine which interactions have stronger effects on the final equilibrial levels of the weaker, species. Precise conditions are stated for the co-existence of both species. Both species are more likely to co-exist in the presence of an intermediate disturbance frequency.

Notes: Abstract This work addresses the design of LoG filters in the frequency domain within a structure formed by the cascade of quasi-Gaussian and discrete Laplacian filters. The main feature of such a structure is that it requires half the number of convolutions of the classical structure in which the LoG transfer function is expressed as the sum of two separable transfer functions of 1-D Gaussian and LoG type. Such a perspective allows one to rephrase the design of IIR and FIR filters for edge detection as a frequency domain approximation problem solvable by standard digital filter design tools. The zero-phase IIR solutions have a good performance at low orders and approximation errors practically independent of the aperture parameter. The characteristics of the nearly linear-phase IIR filters solving the problem suggest the consideration of linear-phase FIR filters with zeros constrained on the unit circle. The use of such filters leads to remarkable computational savings with respect to the filters designed by impulse response sampling. The agreement between the edge values obtained by the filters designed according to the scheme proposed in this work and those obtained by standard techniques is very good.

Notes: Abstract We study boundary value problems for the time-harmonic form of the Maxwell equations, as well as for other related systems of equations, on arbitrary Lipschitz domains in the three-dimensional Euclidean space. The main goal is to develop the corresponding theory for Lp-integrable bounday data for optimal values of p's. We also discuss a number of relevant applications in electromagnetic scattering.

Notes: Abstract Beurling's algebra $A^*=\{f:\sum_{k=0}^{\infty} \sup_{k\le |m|} |\hat f (m)| 〈 \infty \}$ is considered. A* arises quite naturally in problems of summability of the Fourier series at Lebesgue points, whereas Wiener's algebra A of functions with absolutely convergent Fourier series arises when studying the norm convergence of linear means. Certainly, both algebras are used in some other areas. A* has many properties similar to those of A, but there are certain essential distinctions. A* is a regular Banach algebra, its space of maximal ideals coincides with $[-\pi,\pi],$ and its dual space is indicated. Analogs of Herz's and Wiener-Ditkin's theorems hold. Quantitative parameters in an analog of the Beurling-Pollard theorem differ from those for A. Several inclusion results comparing the algebra A* with certain Banach spaces of smooth functions are given. Some special properties of the analogous space for Fourier transforms on the real axis are presented. The paper ends with a summary of some open problems.

Notes: Abstract We construct an algebra of left-invariant pseudodifferential operators on SU(2). We require only that the symbols be homogeneous and C2. For Fourier-bandlimited symbols, we derive the expected formulae for composition and commutators and construct an orthonormal basis of common approximate eigenvectors that could be used to study spectral theory. Some remarks on applications to matrices of operators are made.

Notes: Abstract A novel and complete method is presented for constructing the generalized Bezout identity in polynomial form of the transfer matrix from its state-space representation of a reachable and observable system. As a result, one can easily apply well-developed synthesis and analysis theories to systems described in the frequency domain.

Notes: Abstract The aim of this paper is to obtain the uncertain value set in the complex plane for systems with real and complex parameters that are known to lie inside a ball in a weightedl p-norm. It generalizes previously available results and may be used to test the robust stability of polynomials whose coefficients lie in a weighted lp-ball.

Notes: Abstract This paper describes a new approach—a Fixed-Level Least Squares (FLLS) method for linear-phase FIR filter design. It is intended for rejection of the Gibbs phenomenon through the introduction of a set of equally spaced fixed levels in the transition band and subsequent redefinition of the approximated and weighted functions. Detailed mathematical solutions of the problem as well as many examples are given. The results in graphical form are shown as an output of the FLLS software model.

Notes: Abstract In this paper we examine order reduction of parabolic systems using modal truncation. The parabolic distributed system is first approximated using the Galerkin method. The system matrices have a special structure that allows us to find the approximate spectrum of the parabolic system. To do this we compute approximate inverses of tridiagonal, diagonally dominant symmetric matrices. The approximation leads to algorithms of orderO(n), as opposed to traditional algorithms of orderO(n), wheren is the order of the system. Finally, an example is presented to illustrate the proposed algorithm.

Notes: Abstract We give results concerning the problem of approximating the input-output maps of nonlinear discrete-time approximately finite-memory systems. Here the focus is on the linear dynamical parts of the approximating structures, and we give examples showing that these linear parts can be derived from asingle prespecified function that meets certain conditions. This is done in the context of an approximation theorem in which attention is focused on what we call “basic sets.” We also consider the related but very different problem of approximating nonlinear functionals using lattice operations or the usual linear ring operations. For this problem we give criteria, not just sufficient conditions, for approximation on compact subsets of reflexive Banach spaces (any Hilbert space is a reflexive Banach space).

Notes: Abstract New calibration methods are presented for bearing estimation with sensor location uncertainties, which are based on eigenvalue decomposition of the sample covariance matrix and three or more nondisjoint sources in known bearings. The methods can be applied to arbitrary arrays, including linear arrays, require fewer computations, and are suitable for low signal-to-noise ratio (SNR) cases. Computer simulations are given to illustrate the performance of the proposed methods.

Notes: Abstract A doubly recursive algorithm for time domain convolution with a piecewise linear weighting function is presented that combines the speed of a recursive (IIR) digital filter with the flexibility and ease of design of a nonrecursive (FIR) digital filter. The approach approximates the desired FIR weighting function by a sum-of-triangles weighting function. ForL triangles (or triangle pairs for a linear phase filter) the algorithm is of orderLN. The approximation improves with the number of triangles. A significant advantage of the algorithm compared to FFT filtering or direct convolution is that there is no necessity of a tradeoff between frequency response accuracy and computation time per output point as the data spacing decreases in the filtered signal. The computational complexity is dependent on the number of triangles chosen, not the width of the weighting function, so the algorithm is especially effective for filters with an inherently wide FIR weighting function.

Notes: Abstract A solution of the problem of correlation analysis of stochastic amplitude-angle modulated signals is presented. The solution is obtained on the basis of the theory of linear stochastic processes for the general case of nonstationary correlated Gaussian noises in both modulation channels. These noises are regarded as the results of a linear filtration of Gaussian white noise processes. As an example of the analysis a case of band-pass modulation noises is considered.

Notes: Abstract Two-dimensional (2-D) periodically shift-variant (PSV) digital filters are considered. These filters have applications in processing video signals with cyclostationary noise, scrambling digital images, and 2-D multirate signal processing. The filters are formulated in the form of a Givone-Roesser (GR) state-space model with periodic coefficients. This PSV model is then presented in block form as a shift-invariant system that also has the same GR state-space form. This block form has reduced computations and ease of analysis. An algorithm is developed that transforms any given 2-D PSV GR system to its equivalent shift-invariant model. Invertibility of this model is an important consideration, especially in image scrambling applications. A condition is established for the invertibility of the shift-invariant model of the 2-D PSV system. Also, the inverse system can be easily computed from the original. The established results are illustrated with an example.

Notes: Abstract In a recent paper a method is described for constructing certain approximations to a general element in the closure of the convex hull of a subset of an inner product space. This is of interest in connection with neural networks. Here we give an algorithm that generates simpler approximants with somewhat less computational cost.

Notes: Abstract In this paper, we gave an exponential estimate for the solutions of singular systems of differential equations with a delay and establish a comparison inequality for the systems. Based on this we estabish the criteria for asymptotic stability for the composite singular systems.

Notes: Abstract This paper presents the VLSI architectures for three-level correlator design based on 1-μm CMOS technology. The architecture performs very high speed, real-time, three-level cross-correlation of signals. Two architectures, one for serial incoming samples of signals (serial data) and the other for stored signal samples (parallel data), are described in the paper.

Notes: Abstract This paper gives the existence of and bounds for the branch voltages and currents in a given infinite network, which on each branch are related by a maximal monotone resistance function. We give a simpler proof than that of Minty, with smaller bounds, in the finite network case.

Notes: Abstract An efficient procedure is presented for computing the frequency-response templates of uncertain rational transfer functions depending on two independent interval polynomials. It consists of a modified Cohen-Sutherland algorithm for the intersection test of two rectangles and a pivoting algorithm. Because the proposed procedure traces out the boundary of a frequency-response template directly, it is, in general, more efficient than the previously proposed algorithms.

Notes: Abstract In this paper, a novel VLSI algorithm for the computation of a two-dimensional discrete cosine transform is proposed. The 2D-DCT equation can be expressed by the sum of high order cosine functions, and the algorithm can be realized by combining a highly efficient first order recursive structure with some simplified matrix multiplications, which results in highly regular hardware architecture and simple routing. The algorithm has temporal and spatial locality of connection and can be segmentized for pipeline operations, so the computation time is greatly reduced. Owing to the simplicity in hardware structure, it is especially good for VLSI implementation.

Notes: Abstract We examine a generalized matched-filter problem in which the interference is a nonstationary process generated by passing white noise through a general linear time-varying filter. First a matched filter is constructed by transforming the problem into an equivalent formulation involving stationary interference and a time-varying propagation channel. Whereas the response of a time-invariant matched filter is sampled at its peak, the response of this time-varying matched filter is normalized before sampling to account for variations in the signal power. Next a matched filter is constructed using a spectral characterization of the nonstationary interference. This construction is then used to formulate a simplified solution for the case where the rate of variation in the nonstationary interference is sufficiently small. The different solutions are illustrated by a numerical example.

Notes: Abstract The purpose of this paper is to investigate the approximation capability of elliptic basis function (EBF) neural networks. The main results are: (1) A necessary and sufficient condition for a functionS′ (R 1) to be qualified as an activation function in the hidden layer of an EBF neural network is given. (2) Every nonzero functionG ε L 2(R n ) is qualified to be an activation function in the elliptic neural network to approximate any function in L2(Rn). (3) As applications, we give new proofs of the theorems concerning the approximation capability of affine basis function (ABF) neural networks and generalized radial basis function (GRBF) neural networks (including radial basis function neural networks) with arbitrary activation functions. In particular, we solve the problem of the approximation capability of sigma-pi neural networks.

Notes: Abstract The frequency response of linear interval systems under PID controllers is analyzed. A two-stage elimination algorithm is presented to describe the boundary of the response in the complex plane. Illustrative examples are included.

Notes: Abstract This article deals with a new intelligent analog circuit design system. Unlike previous design methods, this approach introduces a formal circuit representation for both the numerical and heuristic knowledge of the design system. The predicate logic circuit representation is proposed as a new formal analog circuit description language. The language's syntax and semantics provide a precise symbolic description of analog circuit functionality at higher levels of hierarchy and network component connectivities, together with CMOS transistor sizes as the lowest level of hierarchy. It is shown how sentence conversion rules of language grammer can be used to derive transistor level circuits from input performance specifications. Language sentences have the form of clauses. In addition to clause representation, a frame representation to reflect VLSI design hierarchy is also introduced. An original bidirectional inference mechanism with elements of hypothesis has been introduced to infer designs from the knowledge in clause and frame representations. The unique feature of this circuit knowledge representation is its ability to automate the analog CMOS circuit design process. The design methodology is described in detail. The proposed iterative closed loop design system adopts an expert system approach to provide topological synthesis, uses a SPICE circuit simulator to evaluate the circuit performance, and uses a new diagnostic expert system to provide advice on how to improve the design. The implementation of the methodology and associated experimental results for simple CMOS operational amplifier designs are also presented.

Notes: Conclusion For all these reasons, the rediscovery of the equivalence theorem, first stated by David Ricardo in 1817, must be rejected as an adequate basis for policy, just as Ricardo had denied its applicability to the real world. Correspondingly, the concern with the adverse consequences of unfunded social insurance wealth for the supply of national saving, capital intensity, and living standards remains well founded. p ]If, as a practical matter, public pension and social security programs will never be funded actuarially in the United States and most other postindustrial countries, then government-supervised substitution of private for public retirement plans is the only way to achieve at least partial funding. If such substitution follows the British model of allowing employers to contract out of the earnings-related part of the state scheme if equivalent pensions are provided by the company plan, payroll taxes and social security wealth decline so as to reduce their adverse impact on capital formation.

Notes: Summary: Social Policy in the Italian Economy Favourable social and economic conditions constitute the essential framework for a stable development of savings. Saving in the form of insurance becomes advantageous for the individual, and private insurance can thus extend its activity, when social attitudes and the economic situation favour the propensity to save. If conditions change, the State can take over the coverage of risks through social insurance. By means of this institutions, an anti-cyclical policy can be pursued: the amount of social security contributions, for instance, can be increased during the expansion of the cycle and the amounts thus accumulated can be used to grant benefits during the recession period, when contributions can be fixed at a lower percentage of wages. Another type of policy can be pursued by government authorities: that of adjusting social security contributions to industrial profits, thereby directing the subsequent effects on economic growth. Inflation cab cause instability in decisional policies taken by private insurance companies. A solution to the unbalanced increase of costs can be found in index-linking. Life policies of this kind, for instance, cab be closely related to investments in houses, to be bought by the insured themselves in price-linked instalments. After a reference to present developments regarding risk instability and to possibilities of new forms of insurance, this paper considers the competition resulting from the opening of the EEC insurance markets as an opportunity for the Italian market to strengthen its structures.

Notes: Abstract Non-steady-state equations for kidney models are stated. General conservation relations for these equations are derived. Transient equations for the central core model of the renal medulla are developed. Solution of the equations by Laplace transform methods for time invariant volume flows is discussed. The general theory of solving models with time dependent flows by finite difference methods is developed.

Notes: Abstract The behavior of large systems of randomly-interacting variables is examined using an intentionally simplified model. The stable positive solutions are found to exhibit to a significant degree some well-known properties of ecological systems. This resemblance (including for example the predominance of “predator-prey” interactions) is all the more striking in view of the lack of biological “data” at the input end. The findings suggest it advisable to distinguish two kinds of properties in ecosystems. One kind would depend on specifically biological mechanisms; the other would characterize a wide class of persistent systems, and arise from the need for a dynamic balance between positive and negative feedback.

Notes: Abstract For the tumor model of Skipper and Zubrod, which has been analyzed previously for the theoretical FLM function and the effect of chemotherapy against tumors of known or assumed kinetic characteristics, the theoretical continuous labeling (CL) function is derived by considering an equivalent tumor (in terms of unlabeled cell populations) in which the density function of phase duration of cells inS-phasef 2(a 2)=δ(a 2−∞) and the loss functionL 2(t)=0. This mathematical concept of blocking is applied to the analysis of synchronization in tumor growth and blocking effects in cancer chemotherapy. These effects of chemical agents on the cell cycle progression are being incorporated into a previously written computer simulation program for cancer chemotherapy. Whereas, a program is written and used to simulate the CL functions for L1210 leukemia, and primary and metastatic Lewis lung carcinoma.

Notes: Abstract Global stability is established in a class of prey-predator models. This includes a prey-predator model in which the predator has Type 2 functional response and no intraspecific interactions. Two simple examples demonstrate that Kolmogoroff’s theorem does not apply to some members of this class of models.

Notes: Abstract The skeletal muscle is regarded as a periodically deformed plate. An equation of energy for the biological tissue in the system is used in the Lagrange variables. With the heat exchange through the above tissues to the exterior media and also with account of some relations between the physiological factors the approximate analytical solution for this equation is obtained and compared with the physiological data.

Notes: Abstract One of Bobisud's (1976) models for the evolution, of cannibalism is discussed. His analysis is criticised for not being based on the principle of individual selection. Assuming the operation of that principle, we show by simulating his model that cannibals may establish themselves in a noncannibal population. This will happen both in cases where Bobisud concluded cannibalism to be optimal and in cases where he concluded cannibalism not to be optimal.

Notes: Abstract The D'Arcy Thompson concept of biological transformations is developed in a form analogous to such physical concepts as the Law of Corresponding States in thermodynamics, and the Principles of Similitude found in engineering. We find that such concepts depend on a distinction between fundamental and derived quantities, in which the values assigned to the fundamental quantities set the natural scales for the derived ones. Among other things, we see that critical phenomena, such as phase transitions, arise as an immediate consequence of this distinction. In a biological context, we explore the implications of Thompson's hypothesis that closely related organisms are phenotypically similar, assuming that the organisms we see are the result of selection processes operating on phenotypes.

Notes: Abstract If information is coded in the combination of activities of many neurons operating in parallel, then information present in a network can be defined by the correlation of present network activity with the activity which had been elicited by a stimulus in the past; a high correlation indicates the presence of the previously encoded stimulus. Information is distributed in the network because coding is dependent upon the activities of all cells. A model based on Hartline-Ratliff lateral inhibition with a time delay shows that lateral inhibition can distribute information across a parallel network, reduce output noise, and also briefly store information. With no changes in model parameters, and the use of a correlation measure for recognition, the model can stimulate psychophysical results in eleven variations of the metacontrast masking paradigm.

Notes: Abstract A theory is derived to calculate the regional and total deposition of aerosol particles in the nasal passages during inhalation. The particle size studied range from 0.2 to 10.0 μm diameter. The deposition is calculated in five regions; (I) the region filled with nasal hair, (II) the nasal valve, (III) the expansion region, (IV) the turbinate region and (V) the posterior bend. Equations are derived to determine the deposition caused by direct impaction on the nasal hairs and bends of the passages. The calculations show the deposition due to direct impaction does not account for the amount or location of deposited particles measured in experiments. Secondary flows have been speculated to exist in the expansion region after the nasal valve and an equation is derived to estimate the deposition caused by the secondary flows. The calculated deposition, due to direct impaction and secondary flows, shows general agreement with the experiment as to the predicted amount and location of deposited particles.

Notes: Abstract The paper reviews a series of models for circadian clocks and discusses their conclusions and predictions. Attention is focused on Pittendrigh's empirical model, two mathematical models by the author and Winfree's work.

Notes: Abstract From a model of facilitated ionic transport across axonal membranes proposed by McIlroy (1975), the equipotentials and electric field distributions in the vicinity of a potassium conducting pore and of a sodium conducting pore are computed and presented as two-dimensional mappings. The model is then extended to include the effect of impurity ions in the conducting pores viz. of potassium ions in a sodium pore and of sodium ions in a potassium pore. The ionic selectivities and permeabilities of the transported species are discussed in relation to the extended model. Bounds are deduced for the ionic selectivity coefficients for both the sodium and potassium current-carrying systems in squid giant axons and the electric-field distributions in the vicinities of the pores are computed for the extended model and compared with the impurity-free fields first calculated. Finally the permeability coefficients defined in terms of the extended model are shown to reconcile the results of attempts to measure permeability by means of radioactive tracer techniques, with the classical description of the resting nerve.

Notes: Abstract The partial differential equations describing transcapillary exchange and subsequent removal of solutis from an idealized liver sinusoid are amenable to solution by similarity analysis. The power and utility of this technique, which is not widely appreciated as a method for solving biological models, is illustrated here for a system whose Laplace transforms is difficult to invert.

Notes: Abstract The phenomenological theory for the chemotaxis and consumption of oxygen by motile aerobic bacteria is shown to yield a remarkably simple one-dimensional steady-state solution for a congregation of bacteria close to the surface of an oxygen-depleted aqueous medium.

Notes: Abstract In this paper biological compartmental models are considered which take into account the intrinsic randomness of the transport rate parameters and their possible variability in time. An identification procedure is presented for the estimation of the stochastic processes representing the transport rate parameters of a compartmental model from a noisy input-output experiment. The problem is formulated in terms of nonlinear filtering. A simple model is discussed for the case in which the transport rate parameters are independent of each other. The possibility of testing possible important features of the behavior of the transport rate parameters is also evidenced.

Notes: Abstract For any essentially nonlinear system of reaction-diffusion equations of the generic form ∂ci/∂t=Di∇2ci+Qi(c,x,t) supplemented with Robin type boundary conditions over the surface of a closed bounded three-dimensional region, it is demonstrated that all solutions for the concentration distributionn-tuple function c=(c 1(x,t),...,c n (x,t)) satisfy a differential variational condition. Approximate solutions to the reaction-diffusion intial-value boundary-value problem are obtainable by employing this variational condition in conjunction with a Galerkin-Ritz procedure. It is shown that the dynamical evolution from a prescribed initial concentrationn-tuple function to a final steady-state solution can be determined to desired accuracy by such an approximation method. The variational condition also admits a systematic Galerkin-Ritz procedure for obtaining approximate solutions to the multi-equation elliptic boundary-value problem for steady-state distributions c=−c(x). Other systems of phenomenological (non-Lagrangian) field equations can be treated by Galerkin-Ritz procedures based on analogues of the differential variational condition presented here. The method is applied to derive approximate nonconstant steady-state solutions for ann-species symbiosis model.

Notes: Abstract The results of an earlier effort to provide a geometrical analysis of Hutchinsonian niche space are extended. The concept of diversity of a species in niche space is introduced and the maximization of this diversity provides a rationale for a within-niche fitness distribution which is Gaussian. Niche expansion is seen as a consequence of diffusion in niche space, and an evolutionary version of the Volterra competition equations is proposed as a way to relate niche geometry with population dynamics. Applications to topics in community evolution, species packing and the statistical fitting of species abundance data are mentioned.

Notes: Abstract In this paper we use marginal probabilities to derive expressions for the means, variances and covariances ofm-compartment systems. We also present an efficient algorithm for the estimation of the parameters of the system using time series data when measurements are available fromk of them compartments. An application of the analysis and parameter estimation procedure for a model representing the results of a cancer treatment follow-up study is given.

Notes: Abstract The self-organizing properties of an ensemble of interconnected units are studied by linear stability analyses. Small perturbations of a uniform steady-state may result in bifurcations to other solutions that exhibit spatial or temporal order. We show that increasing the number of connections that a unit makes with its neighbors changes the nature of these solutions and tends to destroy spatiotemporal patterns. If an unconnected system is orginally stable, the formation of multiple interconnections can never induce temporal periodicity but may, under certain circumstances, allow the emergence of stationary spatial patterns. We have verified the predictions of the linear stability analysis on a model system and comment on the implications of these results for multicellular ensembles.

Notes: Abstract A performance criterion and weighting factors for the optimal cardiac assistance are investigated by applying Tellegen's network theorem and tolerance analysis on animal experimental data for left ventricular (LV) bypass on the failing heart. Two major factors with respect to cardiac assistance (total power delivered to the peripheral circulatory system, and changes in temporal pattern of ventricular contraction) are represented by two performance criteria,J 1 andJ 2 whereJ 1 relates to the sum of LV and pump power, andJ 2 relates to the “peakedness” factor of LV power. The total performance index (J) is determined as the weighted sum ofJ 1 andJ 2;J=w 1J1+w2J2. The weighting factors,w 1 andw 2, are computed as inverses of the tolerance in the performance contours with respect to improvement of stroke work per minute from pre- to post-bypass condition.

Notes: Abstract In treating the Volterra-Verhulst prey-predator system with time dependent coefficients, we ask how far this deterministic system represents or approximates the dynamics of the population evolving in a realistic environment which is stochastic in nature. We consider a stochastic system withsmall Gaussian noise type fluctuations. It is shown that the higher moments of the deviation of the deterministic system from the stochastic approach zero as the strength δ of the perturbation decays to zero. For any δ〉0 and allT〉0, ε〉0, the sample population paths that stay within ε distance from the deterministic path during [0,T] form a collection of positive probability. In comparing the stationary distributions of the two systems, we show that the weak limits of those of the stochastic system form a subset of those of the deterministic system. This is in analogy with a result of May connected with the stability of the two systems. Plant and rodent populations possess periodic parameters andexhibit periodic behaivor. We establish theoretically this periodicity under periodicity conditions on the coefficients and perturbing random forces. We also establish a central limit property for the prey-predator system.

Notes: Abstract Adaptation of repetitively firing sensory neurons and nerve models is correlated with specific inhibitory feedback phenomena—an electrogenic sodium pump, and post synaptic self inhibition. The quality of the adaptive responses depends on the excitation properties of the neuron in the interspike interval, or the description of these properties by the underlying impulse encoder model. THis model dependence is demonstrated by comparisons of the behavior of two classes of models; the “leaky integrator models” which assume a passive neural membrane, and the “variable-γ models”, for which the neural state of excitation varies according to first order differential equations. The complexity inherent in the variable-γ models is effectively boiled down to mathematically simple relationships which are derived from studies of the neural- and model frequency responses to small amplitude sinusoidal stimuli. It is argued, and supported with examples, that these relationships hold for impulse frequency transients resulting from more general stimulus conditions. Expressions are then derived which permit feedback parameters to be determined from impulse frequency data. In this connection, recent studies of neural dynamics are brought to bear to resolve ambiguities in data interpretation.

Notes: Abstract A setQ(n) of noncongruent connected shapes, constructed from a fixed numbern of congruent regular hexagons laid edge to edge, is defined. Various measures of shape are applied toQ(n) and the results are numerically analyzed in the special casesn≦9. The concept of spatial entropy is introduced which affords a measure of the “complexity” of shape. Possible biological applications include the analysis of ecological cover,stigmergy or the complex nest-building activities of social insects and morphogenesis. Essentially any comparative study of shape, regardless of the specific application, might be carried out along the lines suggested in the paper.

Notes: Abstract This paper uses optimal control theory in conjunction with a Gompertzian type model for cellular growth to determine the optimal method of administering cycle non-specific chemotherapy or more generally the optimal durations of treatment and rest periods during chemotherapy. The performance critera employed to determine the relative merits of the therapy include not only the destruction of malignant cells, but also the sparing of a critical normal tissue. Since these criteria are at odds with one another, the solutions are found which satisfy the Pareto optimality conditions.

Notes: Abstract Oscillations governed by generalized Volterra-Gause-Witt equations to include retardation effects in population dynamics are considered. Models containing either small or significant time delay are discussed. The Krylov-Bogoliubov-Mitropolskii perturbation method and its extension for differential equations with retarded argument is used.

Notes: Abstract A general solution to the dynamical equation for the probability distribution associated withn interacting species is obtained by employing the author's generic canonical expression for the rate functions. Interacting species models with limit-cycle dynamics and no stable equilibrium points feature probability distributions that are asymptotic for large values oft to Dirac δ-distributions concentrated on the limit-cycles, as illustrated here for an analytically solvable two-species model. For ann-species Volterra model, a stationary or temporally-averaged probability distribution should generally be much more complicated than the specialized Poisson form studied by Kerner and others.

Notes: Abstract A set of coupled nonlinear differential equations, determining the concentration profiles and electric potentials valid for isothermal transport of ions and molecules across a diffusion barrier are formulated, using a correction to the limiting expression for chemical potential gradients and the molecular expression for frictional force. These differential equations are similar to Nernst-Planck equations and reduce to these under appropriate approximations. Solutions of these equations valid under specified conditions are presented. Expressions for permeability, concentration profiles of many ion systems are included.

Notes: Abstract The global flow equations of nonequilibrium thermodynamics for a single nonelectrolyte solute and water passing through a membrane are obtained by solving the local equations of motion. The method follows that developed for the general,n-solute case in the previous paper (Mikulecky, 1978). It is easily seen in this simple case that the passage from local interactions, formulated as position dependent frictional interactions in the equations of motion, to ghe global result involves a loss of any simple way of identifying particulars about local information. Two particular cases are analyzed in further detail: the case of no interaction within the pore and the case of constant interaction for both solute and solvent across the pore. In the former case, Onsager reciprocity survives in the global result if a self-consistent definition of the partial viscosity coefficients is used, while in the latter case, reciprocity is lost. Since, in many biologically interesting cases, the presence of interaction of the type considered here is likely to occur, the reciprocity condition should not automatically be assumed to hold.

Notes: Abstract This study is related to a model describing the behavior of barium-treatedAplysia neurons generating regular burst-plateau patterns. The model is represented by an autonomous dynamical system, defined inR 4 and depending on a small parameter. This paper is restricted to the qualitative study of three “reduced systems” deduced from the “complete system”. Part of the study is performed with the use of the qualitative theory of singular perturbations. The predicted behaviors are compared with experimental results.

Notes: Abstract Using linear stability analysis, the qualitative stability properties of open nonlinear chemical systems, in which reactions of any order may occur, will be studied. Systems will be classified in three fundamental classes: trees, cycles and loops, according to their knot graphs. The study of the Jacobian matrix for the kinetic equations of the system shows that the symmetrizability by a particular procedure (calledD-symmetrizability) is a sufficient condition for stability. It has been proved that tree-graphs always satisfy the above condition. For the cycle-graphs, theD-symmetrizability condition leads to a cyclic relation between forward and reverse steady state flows. The stability may be assured, even if the cyclic relation is not satisfied, providing that the “symmetry breaking” be lower than an upper bound; further alternative criteria for stability of cycles have been derived. All these results are independent of the number of diffusive exchanges with the environment.

Notes: Abstract The basis of an analytical description of the behaviour of large random nets of binary elements of the type first investigated in detail by S. A. Kauffman is presented. It is shown that information about the network dynamics can be deduced from quite general considerations of the properties of the state transition graph and matrix. An expression for the matrix elements of the state transition matrix in terms of the Boolean function specification of the net is derived. Using these ideas the distribution of limit cycle lengthsl for a completely random net is calculated and shown to bex 1/l, a result which agrees well with experimental data.

Notes: Abstract The state-transition matrix description of Kauffman binary networks described in the previous paper is further developed to obtain an analytical expression for the fraction of states involved in limit cycles as a function of the network size and connectivity. The result obtained for totally connected networks agrees with that derived from quite different considerations by other workers. For low connectivity networks the results are in qualitative agreement with the experimental data of Kauffman but there is a quantitative discrepancy which remains to be resolved.

Notes: Abstract Persistence-extinction in simple food chains modelled by Lotka-Volterra dynamics is governed by a single parameter which depends upon the interspecific interaction coefficients, the intraspecific interaction coefficients, and the length of the food chain. In persistent systems with nonzero carrying capacity, two new features predominate. Trophic level influence factors relate persistence on different trophic levels and determine, in conjunction with the persistence parameter, the magnitude of persistence. Equilibrium component ordering, which results in persistent systems, mandates once again that systems need to be studied on the complete ecosystem level; static field measurements reflect species location in the food chain, the total length of the food chain and assume characteristics according to these factors.

Notes: Abstract The geometrical nature of the Elementary catastrophes (Thom, 1969) is reviewed. Histories of the movement of catastrophe manifolds and bifurcation sets are presented for some of the space-equivalent unfoldings described by Wassermann (1975). These unfoldings provide descriptions of the variation with time of the stability of stationary states of associated potential energy functions. Identification of these stationary energy states with stationary states of a system therefore provides a description of its behavior with time. Qualitative descriptions of this type are particularly useful when the complexity of a system prevents a detailed quantitative description. Histories of bifurcation set movements suggest different types of system behavior at different space-like coordinates. This type of theory may be a useful model for the processes leading to differentiation of cells and to emergence of adult forms of a biological organism.

Notes: Abstract Three conjectures are given which predict the existence of unique stable limit cycle oscillations in a class of piecewise linear (PL) differential equations. The equations are appropriate to model biological or other complex systems in which there are switchlike interactions between the elements of the network. Methods are presented which can be used to develop mathematical models which are conjectured to display stable limit cycle oscillations, from qualitative experimental information about relative phases of activity in the dynamical systems. Several illustrative numerical examples are given, and one experimental example from neurobiology is discussed.

Notes: Abstract The formulation and results of the Kalman State Regulator Problem are applied to a mathematical model of the arterial system of a dog to obtain an optimal control for blood pressure. The criterion for optimality is minimum energy per cycle.

Notes: Abstract A mathematical model is developed in state variable form suitable for the study of the control of blood pressure and flow in the mammalian cardiovascular system. The applicability of the model to steady state, both mean and pulsatile, and transient phenomena is demonstrated by the agreement of the results with experimental data. This model was developed to study the neural and renal-endocrine-electrolyte control of cardiovascular functions.

Notes: Abstract As a contribution to the discussion of oscillatory models for interacting species it is shown that two-species Volterra models can never have limit cycles, and a complete enumeration is given of conditions which the parameters of these models must satisfy in order that a part of the phase space be filled with a family of closed curves; sketches of phase portraits are also given. These results complement and correct older results by Bautin and by Coppel on quadratic differential systems. The paper opens with a brief discussion of some more practical aspects of the ecological application of oscillatory models.

Notes: Abstract The purpose of this paper is to present a method for the determination of a solution for a coupled system of first order initial boundary-value problems arising from some biological systems. The physical problem is to determine the suspended and the superficial molecular concentrations of a traced substance passing through an organ containing a tangle of vessels, such as the kidney-ureter system. The approach to the problem is by successive approximation which leads to a recursion formula for the determination of the solution as well as error estimates for the approximations. The recursion formula involves only direct integration which indicates a promising possibility in obtaining numerical results by using a computer. In addition to the determination of a solution, some qualitative analysis of the solution is given. This includes the existence of a unique solution, the continuous dependency of the solution on the data, and the stability problem of a steady-state solution.

Notes: Abstract A global kinetic analysis of a model consisting of an autocatalytic zymogen-activation process, in which an irreversible inhibitor competes with the zymogen for the active site of the proteinase, and a monitoring coupled reaction, in which the enzyme acts upon one of its substrates, is presented. This analysis is based on the progress curves of any of the two products released in the monitoring reaction. The general solution is applied to an important particular case in which rapid equilibrium conditions prevail. Finally, we suggest a procedure to predict whether the inhibition or activation route dominates in the steady state of the system. These results generalize our previous analysis of simpler mechanisms.