Magistro darbe nagrinėjami paprastųjų tiesinių antros eilės diferencialinių lygčių sistemos su ypatuma kraštiniai uždaviniai. Ištirta sprendinių asimptotika ypatingojo taško aplinkoje. Surasti lygčių sistemos sprendinio įverčiai. Nagrinėjamų diferencialinių uždavinių sprendiniai, konstruojami naudojant integralinius operatorius. Įrodyta sprendinių vienatis. / In this paper the some boundary value problems for a system of ordinary second order differential equations with singularity are considered. The asymptotic of solutions in the neighborhood of singular point are discussed. The estimates of the solution are given. The solutions of considered differential problems are constructed using some integral operators. The uniqueness of the solutions is proved.

Mathematical models provide unique insights to real-world problems. Within the context of infectious diseases, models are used to explore the dynamics of infections and control mechanisms. Human papillomavirus (HPV) globally infects about 630 million people, many of these infections develop into cancers and genital warts. Vaccines are available to protect against the most prevalent and devastating strains of HPV. The introduction of this vaccine as part of a national immunization program in Canada is a complex decision for policy-makers in which mathematical models can play a key role. We use the current recommendations provided by the World Health Organization to explore the integral role mathematical models have in the decision to incorporate the HPV vaccine within a national immunization program. We then provide a review of the literature discussing the role of mathematical models in the decision to include a vaccine in a national immunization program within the context of the HPV vaccine. Next, we evaluate the current standing of mathematical models used within the context of HPV immunization, to highlight the types of models used, underlying assumptions and general recommendations made about these immunization programs. Then, we create and analyze a model to explore the possibility of bettering the current HPV vaccine strategy in Canada. We focus on the effects of the grade of vaccination and the number of doses required to eradicate the targeted strains of HPV.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-07-06T14:21:44Z No. of bitstreams: 1 MOISÉS DANTAS DOS SANTOS - DISSERTAÇÃO PPGMAT 2005..pdf: 524230 bytes, checksum: a8d6e23eaf6da89e369eb29e888e7a1a (MD5) / Made available in DSpace on 2018-07-06T14:21:44Z (GMT). No. of bitstreams: 1 MOISÉS DANTAS DOS SANTOS - DISSERTAÇÃO PPGMAT 2005..pdf: 524230 bytes, checksum: a8d6e23eaf6da89e369eb29e888e7a1a (MD5) Previous issue date: 2005-12 / Neste trabalho usaremos métodos variacionais para mostrar a existência de solução fraca para dois tipos de problema. O primeiro consiste num problema não-linear, O segundo, trata-se de uma Equação Diferencial Ordinária / In this work we use variational methods to show the existence of weak solutions for two types problems. The first of them is a nonlinear problem, The second, is related with a following Ordinary Differential Equations.

Interval analysis is an essential tool in the construction of validated numerical solutions of Initial Value Problems (IVP) for Ordinary (ODE) and Partial (PDE) Differential Equations. A validated solution typically consists of guaranteed lower and upper bounds for the exact solution or set of exact solutions in the case of uncertain data, i.e. it is an interval function (enclosure) containing all solutions of the problem. IVP for ODE: The central point of discussion is the wrapping effect. A new concept of wrapping function is introduced and applied in studying this effect. It is proved that the wrapping function is the limit of the enclosures produced by any method of certain type (propagate and wrap type). Then, the wrapping effect can be quantified as the difference between the wrapping function and the optimal interval enclosure of the solution set (or some norm of it). The problems with no wrapping effect are characterized as problems for which the wrapping function equals the optimal interval enclosure. A sufficient condition for no wrapping effect is that there exist a linear transformation, preserving the intervals, which reduces the right-hand side of the system of ODE to a quasi-isotone function. This condition is also necessary for linear problems and "near" necessary in the general case. Hyperbolic PDE: The Initial Value Problem with periodic boundary conditions for the wave equation is considered. It is proved that under certain conditions the problem is an operator equation with an operator of monotone type. Using the established monotone properties, an interval (validated) method for numerical solution of the problem is proposed. The solution is obtained step by step in the time dimension as a Fourier series of the space variable and a polynomial of the time variable. The numerical implementation involves computations in Fourier and Taylor functoids. Propagation of discontinuo~swaves is a serious problem when a Fourier series is used (Gibbs phenomenon, etc.). We propose the combined use of periodic splines and Fourier series for representing discontinuous functions and a method for propagating discontinuous waves. The numerical implementation involves computations in a Fourier hyper functoid. / Mathematical Sciences / D. Phil. (Mathematics)

In this thesis we address the analytic and numerical aspects of isospectral flows. Such flows occur in mathematical physics and numerical linear algebra. Their main structural feature is to retain the eigenvalues in the solution space. We explore the solution of Isospectral flows and their stochastic counterpart using explicit generalisation of Magnus expansion. \par In the first part of the thesis we expand the solution of Bloch--Iserles equations, the matrix ordinary differential system of the form $ X'=[N,X^{2}],\ \ t\geq0, \ \ X(0)=X_0\in \textrm{Sym}(n),\ N\in \mathfrak{so}(n), $ where $\textrm{Sym}(n)$ denotes the space of real $n\times n$ symmetric matrices and $\mathfrak{so}(n)$ denotes the Lie algebra of real $n\times n$ skew-symmetric matrices. This system is endowed with Poisson structure and is integrable. Various important properties of the flow are discussed. The flow is solved using explicit Magnus expansion and the terms of expansion are represented as binary rooted trees deducing an explicit formalism to construct the trees recursively. Unlike classical numerical methods, e.g.\ Runge--Kutta and multistep methods, Magnus expansion respects the isospectrality of the system, and the shorthand of binary rooted trees reduces the computational cost of the exponentially growing terms. The desired structure of the solution (also with large time steps) has been displayed. \par Having seen the promising results in the first part of the thesis, the technique has been extended to the generalised double bracket flow $ X^{'}=[[N,X]+M,X], \ \ t\geq0, \ \ X(0)=X_0\in \textrm{Sym}(n),$ where $N\in \textrm{diag}(n)$ and $M\in \mathfrak{so}(n)$, which is also a form of an Isospectral flow. In the second part of the thesis we define the generalised double bracket flow and discuss its dynamics. It is noted that $N=0$ reduces it to an integrable flow, while for $M=0$ it results in a gradient flow. We analyse the flow for various non-zero values of $N$ and $M$ by assigning different weights and observe Hopf bifurcation in the system. The discretisation is done using Magnus series and the expansion terms have been portrayed using binary rooted trees. Although this matrix system appears more complex and leads to the tri-colour leaves; it has been possible to formulate the explicit recursive rule. The desired structure of the solution is obtained that leaves the eigenvalues invariant in the solution space.

Being part of a research group on functional differential equations (FDEs, for short), due to my formation in non-absolute integration theory and because certain kinds of FDEs can be expressed as integral equations, I was motivated to investigate the latter. The purpose of this work, therefore, is to develop the theory of integral equations, when the integrals involved are in the sense of Kurzweil- Henstock or Kurzweil-Henstock-Stieltjes, through the correspondence between solutions of integral equations and solutions of generalized ordinary differential equations (we write generalized ODEs, for short). In order to be able to obtain results for integral equations, we propose extensions of both the Kurzweil integral and the generalized ODEs (found in [36]). We develop the fundamental properties of this new generalized ODE, such as existence and uniqueness of solutions results, and we propose stability concepts for the solutions of our new class of equations. We, then, apply these results to a class of nonlinear Volterra integral equations of the second kind. Finally, we consider a model of population growth (found in [4]) that can be expressed as an integral equation that belongs to this class of nonlinear Volterra integral equations. / Sendo parte de um grupo de pesquisa em equações diferenciais funcionais (escrevemos EDFs), por causa de minha formação em teoria de integração não absoluta e porque certos tipos de EDFs podem ser escritas como equações integrais, decidi estudar esse último tipo de equações. O objetivo desse trabalho, portanto, é desenvolver a teoria de equações integrais, quando as integrais envolvidas são no sentido de Kurzweil-Henstock ou Kurzweil-Henstock-Stieltjes, através da correspondência entre soluções de equações integrais e soluções de equações diferenciais ordinárias generalizadas (ou EDOs generalizadas). A fim de obter resultados para estas equações integrais, propomos extensões de ambas a integral de Kurzweil e as EDOs generalizadas (encontradas em [36]). Desenvolvemos propriedades fundamentais dessa nova EDO generalizada, como resultados de existência e unicidade de solução, e propomos conceitos de estabilidade para as soluções de nossa nova classe de equações. Nós, então, aplicamos esses resultados a uma classe de equações integrais de Volterra não lineares de segunda espécie. Finalmente, consideramos um modelo de crescimento de populações (encontrado em [4]) que pode ser escrito como uma equação integral pertencente a essa classe de equações integrais de Volterra não lineares.

Submitted by JULIANA FROES FORSAN (juliana_ff@hotmail.com) on 2017-09-24T22:08:30Z No. of bitstreams: 1 1_Juliana_Froes_Forsan_Dissertacao_Mestrado_Profmat.pdf: 48728620 bytes, checksum: ca33b6f6fe78c0be5e8eea5d84e2c681 (MD5) / Approved for entry into archive by Monique Sasaki (sayumi_sasaki@hotmail.com) on 2017-09-27T19:55:50Z (GMT) No. of bitstreams: 1 forsan_jf_me_rcla.pdf: 48728620 bytes, checksum: ca33b6f6fe78c0be5e8eea5d84e2c681 (MD5) / Made available in DSpace on 2017-09-27T19:55:50Z (GMT). No. of bitstreams: 1 forsan_jf_me_rcla.pdf: 48728620 bytes, checksum: ca33b6f6fe78c0be5e8eea5d84e2c681 (MD5) Previous issue date: 2017-09-01 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Neste trabalho, apresentamos um estudo sobre equações diferenciais ordinárias com o objetivo de compreender alguns modelos matemáticos abordados a nível de ensino médio, como o modelo de crescimento populacional segundo Malthus, Lei do decaimento radioativo, lei de resfriamento de Newton e sistema massa mola ideal. Neste sentido, selecionamos alguns conceitos e resultados matemáticos sobre cálculo diferencial e integral, introduzindo o estudo sobre limite, derivada e, de forma breve, sobre integração. Apresentamos o número de Euler (e) e trabalhamos alguns exercícios envolvendo os modelos citados. Para o desenvolvimento dos assuntos e das atividades propostas, procuramos abordar conceitos em física e utilizamos os softwares GeoGebra e Modellus bem como simuladores disponíveis na internet. / In this work, we present a study on ordinary differential equations with the objective of understanding some mathematical models addressed at secondary level, such as Malthus' theory of population growth, Law of radioactive decay, Newton's law of cooling and ideal mass spring system. In this sense, we have selected some concepts and mathematical results on differential and integral calculus, introducing the study on limit, derivative, and, briefly, on integration. We present the Euler's number (e) and we work some exercises involving the mentioned models. For the development of the subjects and the proposed activities, we try to approach concepts in physics and we use the software GeoGebra and Modellus as well as simulators available on the internet. / CAPES: 5512038

Made available in DSpace on 2015-03-03T15:32:43Z (GMT). No. of bitstreams: 1 MarcioLN_DISSERT.pdf: 2785073 bytes, checksum: 8f894388b11b263c73e967b4b680c52f (MD5) Previous issue date: 2013-08-13 / Coordena??o de Aperfei?oamento de Pessoal de N?vel Superior / This paper has two objectives: (i) conducting a literature search on the criteria of uniqueness of solution for initial value problems of ordinary differential equations. (ii) a modification of the method of Euler that seems to be able to converge to a solution of the problem, if the solution is not unique / O presente trabalho tem dois objetivos: (i) a realiza??o de uma pesquisa bibliografifica sobre os crit?rios de unicidade de solu??o para problemas de valor inicial de equa??es diferenciais ordin?rias. (ii) Introduzir uma modifica??o do m?todo de Euler que parece ser capaz de convergir a uma das solu??es do problema, caso a solu??o n?o seja ?nica

abstract: In the field of infectious disease epidemiology, the assessment of model robustness outcomes plays a significant role in the identification, reformulation, and evaluation of preparedness strategies aimed at limiting the impact of catastrophic events (pandemics or the deliberate release of biological agents) or used in the management of disease prevention strategies, or employed in the identification and evaluation of control or mitigation measures. The research work in this dissertation focuses on: The comparison and assessment of the role of exponentially distributed waiting times versus the use of generalized non-exponential parametric distributed waiting times of infectious periods on the quantitative and qualitative outcomes generated by Susceptible-Infectious-Removed (SIR) models. Specifically, Gamma distributed infectious periods are considered in the three research projects developed following the applications found in (Bailey 1964, Anderson 1980, Wearing 2005, Feng 2007, Feng 2007, Yan 2008, lloyd 2009, Vergu 2010). i) The first project focuses on the influence of input model parameters, such as the transmission rate, mean and variance of Gamma distributed infectious periods, on disease prevalence, the peak epidemic size and its timing, final epidemic size, epidemic duration and basic reproduction number. Global uncertainty and sensitivity analyses are carried out using a deterministic Susceptible-Infectious-Recovered (SIR) model. The quantitative effect and qualitative relation between input model parameters and outcome variables are established using Latin Hypercube Sampling (LHS) and Partial rank correlation coefficient (PRCC) and Spearman rank correlation coefficient (RCC) sensitivity indices. We learnt that: For relatively low (R0 close to one) to high (mean of R0 equals 15) transmissibility, the variance of the Gamma distribution for the infectious period, input parameter of the deterministic age-of-infection SIR model, is key (statistically significant) on the predictability of the epidemiological variables such as the epidemic duration and the peak size and timing of the prevalence of infectious individuals and therefore, for the predictability these variables, it is preferable to utilize a nonlinear system of Volterra integral equations, rather than a nonlinear system of ordinary differential equations. The predictability of epidemiological variables such as the final epidemic size and the basic reproduction number are unaffected by (or independent of) the variance of the Gamma distribution for the infectious period and therefore for the choice on which type of nonlinear system for the description of the SIR model (VIE's or ODE's) is irrelevant. Although, for practical proposes, with the aim of lowering the complexity and number operations in the numerical methods, a nonlinear system of ordinary differential equations is preferred. The main contribution lies in the development of a model based decision-tool that helps determine when SIR models given in terms of Volterra integral equations are equivalent or better suited than SIR models that only consider exponentially distributed infectious periods. ii) The second project addresses the question of whether or not there is sufficient evidence to conclude that two empirical distributions for a single epidemiological outcome, one generated using a stochastic SIR model under exponentially distributed infectious periods and the other under the non-exponentially distributed infectious period, are statistically dissimilar. The stochastic formulations are modeled via a continuous time Markov chain model. The statistical hypothesis test is conducted using the non-parametric Kolmogorov-Smirnov test. We found evidence that shows that for low to moderate transmissibility, all empirical distribution pairs (generated from exponential and non-exponential distributions) for each of the epidemiological quantities considered are statistically dissimilar. The research in this project helps determine whether the weakening exponential distribution assumption must be considered in the estimation of probability of events defined from the empirical distribution of specific random variables. iii) The third project involves the assessment of the effect of exponentially distributed infectious periods on estimates of input parameter and the associated outcome variable predictions. Quantities unaffected by the use of exponentially distributed infectious period within low transmissibility scenarios include, the prevalence peak time, final epidemic size, epidemic duration and basic reproduction number and for high transmissibility scenarios only the prevalence peak time and final epidemic size. An application designed to determine from incidence data whether there is sufficient statistical evidence to conclude that the infectious period distribution should not be modeled by an exponential distribution is developed. A method for estimating explicitly specified non-exponential parametric probability density functions for the infectious period from epidemiological data is developed. The methodologies presented in this dissertation may be applicable to models where waiting times are used to model transitions between stages, a process that is common in the study of life-history dynamics of many ecological systems. / Dissertation/Thesis / Ph.D. Applied Mathematics for the Life and Social Sciences 2014

Afin de faire face à différents types de stress et s'adapter à de nouveaux environnements, les bactéries ont développé de nombreux mécanismes génétiquement régulés. La compétence pour la transformation naturelle est un processus qui favorise le transfert horizontal de gènes. Si les espèces phylogénétiquement éloignées partagent des mécanismes conservés d'intégration et de remaniement de l'ADN, les circuits de régulation de la compétence ne sont toutefois pas universels mais adaptés au mode de vie de chaque espèce. Chez les bactéries Gram-positives, les cascades de régulation de Streptococcus pneumoniae et Bacillus subtilis sont les mieux documentées. Si de nombreux modèles mathématiques ont été établis pour étudier différents aspects de la régulation des compétences chez B. subtilis, un seul modèle à échelle de population a été développé pour S. pneumoniae, il y a plus de dix ans, sur la base d'hypothèses contestées par de nouvelles données expérimentales. Nous avons développé, chez S. pneumoniae, un modèle fondé sur la connaissance de la régulation de la compétence qui intègre les éléments biologiques essentiels connus à ce jour. La cohérence structurelle de la topologie du réseau est confirmée par le formalisme des réseaux de Petri. Le réseau est ensuite transformé en un ensemble d'équations différentielles ordinaires pour étudier son comportement dynamique. La cinétique des protéines a été estimée en utilisant des données de luminescence et l'estimation des paramètres a été contrainte à partir des connaissances disponibles. Après avoir testé des modèles alternatifs, nous avons proposé l'existence d'un produit de gène tardif supplémentaire pouvant inhiber l'action de ComW, l'activateur du facteur sx. Nous apportons également un nouvel éclairage sur cette cascade de régulation en prédisant la cinétique de composantes du système qui pourraient être impliquées dans des comportements spécifiques. Ce modèle consolide les connaissances expérimentales acquises sur la régulation de la compétence chez S. pneumoniae. De plus, il peut être appliqué aux autres espèces de streptocoques appartenant aux groupes mitis et anginosus puisqu'ils partagent le même circuit régulateur. À l'échelle populationnelle, la transition vers l'état de compétence se produit d'abord dans une sous-population de cellules et se propage ensuite dans toute la population par contact physique cellule à cellule. En permettant la simulation du comportement d'une cellule individuelle, le modèle pourra servir de module dans la conception d'un modèle d'une population bactérienne composée de cellules hétérogènes. / Bacteria have evolved many types of genetically induced mechanisms to face different types of stresses and to adapt to new environments. Competence for natural transformation is one such process that promotes horizontal gene transfer. If phylogenetically distant species share conserved uptake and processing apparatus, competence regulatory circuits are not universal but adapted to every species' lifestyle. In Gram-positive bacteria, Streptococcus pneumoniae and Bacillus subtilis regulatory cascades are the best documented. If many mathematical models have been established to study different aspects of competence regulation in B. subtilis, only one population-scaled model has been developed for S. pneumoniae, a decade ago, based on hypotheses that are challenged by new experimental data. We develop, in S. pneumoniae, a knowledge-based model of the competence regulation at cell level that integrates the enriched biological knowledge acquired to date. The structural consistency of the network topology is confirmed using Petri net formalism. The network is further turned into a set of ordinary differential equations to study its dynamics behavior. Protein kinetics are estimated using time-series luminescence data and other parameter estimations are constrained according to available knowledge. We point out some gap in competence shut-off knowledge, and, after testing alternative models, we predict the requirement of a yet unknown late com gene product inhibiting the action of ComW, the ?x factor activator. We also bring new insights into this regulatory cascade by predicting the system components that might be involved in specific experimental behavior. Our model consolidates the experimental knowledge acquired on competence regulation in S. pneumoniae. Moreover, it can be applied to the other streptococci species belonging to the mitis and anginosus groups since they shared the same regulatory circuit. In the population, the competence shift happens first in a subpopulation of cells and spreads into the whole population through cell to cell contact. Allowing simulation of individual cell behavior, our model will provide a brick for the design of a population-scale model composed of heterogeneous cells.