Mathematics
We are an active university Mathematics Department with a strong teaching and research reputation. We offer students the chance to study at undergraduate or postgraduate level on degree programmes leading to: BSc in Mathematics, BSc/BA joint courses in Mathematics or Applied Statistics and a wide range of other subjects. We have an active research group focusing on Computational Applied Mathematics, with research students studying for the degrees of MPhil and PhD, postdoctoral workers and associated collaborators from across the world.
This collection is licenced under a Creative Commons licence. The collection may be reproduced for noncommerical use and without modification, providing that copyright is acknowledged.
Recent Submissions

Developing A Highperformance Liquid Chromatography Method for Simultaneous Determination of Loratadine and its Metabolite Desloratadine in Human Plasma.Allergic diseases are considered among the major burdons of public health with increased prevalence globally. Histamine H1receptor antagonists are the foremost commonly used drugs in the treatment of allergic disorders. Our target drug is one of this class, loratadine and its biometabolite desloratadine which is also a non sedating H1 receptor antagonist with antihistaminic action of 2.5 to 4 times greater than loratadine. To develop and validate a novel isocratic reversedphase high performance liquid chromatography (RPHPLC) method for rapid and simultaneous separation and determination of loratadine and its metabolite, desloratadine in human plasma. The drug extraction method from plasma was based on protein precipitation technique. The separation was carried out on a Thermo Scientific BDS Hypersil C18 column (5µm, 250 x 4.60 mm) using a mobile phase of MeOH : 0.025M KH2PO4 adjusted to pH 3.50 using orthophosphoric acid (85 : 15, v/v) at ambient temperature. The flow rate was maintained at 1 mL/min and maximum absorption was measured using PDA detector at 248 nm. The retention times of loratadine and desloratadine in plasma samples were recorded to be 4.10 and 5.08 minutes respectively, indicating a short analysis time. Limits of detection were found to be 1.80 and 1.97 ng/mL for loratadine and desloratadine, respectively, showing a high degree of method sensitivity. The method was then validated according to FDA guidelines for the determination of the two analytes in human plasma. The results obtained indicate that the proposed method is rapid, sensitive in the nanogram range, accurate, selective, robust and reproducible compared to other reported methods. [Abstract copyright: Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.]

Double Bordered Constructions of SelfDual Codes from Group Rings over Frobenius RingsIn this work, we describe a double bordered construction of selfdual codes from group rings. We show that this construction is effective for groups of order 2p where p is odd, over the rings F2 + uF2 and F4 + uF4. We demonstrate the importance of this new construction by finding many new binary selfdual codes of lengths 64, 68 and 80; the new codes and their corresponding weight enumerators are listed in several tables

On the behavior of the solutions for linear autonomous mixed type difference equationA class of linear autonomous mixed type difference equations is considered, and some new results on the asymptotic behavior and the stability are given, via a positive root of the corresponding characteristic equation.

Gcodes over Formal Power Series Rings and Finite Chain RingsIn this work, we define $G$codes over the infinite ring $R_\infty$ as ideals in the group ring $R_\infty G$. We show that the dual of a $G$code is again a $G$code in this setting. We study the projections and lifts of $G$codes over the finite chain rings and over the formal power series rings respectively. We extend known results of constructing $\gamma$adic codes over $R_\infty$ to $\gamma$adic $G$codes over the same ring. We also study $G$codes over principal ideal rings.

A Modified Bordered Construction for SelfDual Codes from Group RingsWe describe a bordered construction for selfdual codes coming from group rings. We apply the constructions coming from the cyclic and dihedral groups over several alphabets to obtain extremal binary selfdual codes of various lengths. In particular we find a new extremal binary selfdual code of length 78.

Composite Constructions of SelfDual Codes from Group Rings and New Extremal SelfDual Binary Codes of Length 68We describe eight composite constructions from group rings where the orders of the groups are 4 and 8, which are then applied to find selfdual codes of length 16 over F4. These codes have binary images with parameters [32, 16, 8] or [32, 16, 6]. These are lifted to codes over F4 + uF4, to obtain codes with Gray images extremal selfdual binary codes of length 64. Finally, we use a buildingup method over F2 + uF2 to obtain new extremal binary selfdual codes of length 68. We construct 11 new codes via the buildingup method and 2 new codes by considering possible neighbors.

Numerical methods for solving space fractional partial differential equations by using Hadamard finitepart integral approachWe introduce a novel numerical method for solving twosided space fractional partial differential equation in two dimensional case. The approximation of the space fractional RiemannLiouville derivative is based on the approximation of the Hadamard finitepart integral which has the convergence order $O(h^{3 \alpha})$, where $h$ is the space step size and $\alpha\in (1, 2)$ is the order of RiemannLiouville fractional derivative. Based on this scheme, we introduce a shifted finite difference method for solving space fractional partial differential equation. We obtained the error estimates with the convergence orders $O(\tau +h^{3\alpha}+ h^{\beta})$, where $\tau$ is the time step size and $\beta >0$ is a parameter which measures the smoothness of the fractional derivatives of the solution of the equation. Unlike the numerical methods for solving space fractional partial differential equation constructed by using the standard shifted Gr\"unwaldLetnikov formula or higher order Lubich'e methods which require the solution of the equation satisfies the homogeneous Dirichlet boundary condition in order to get the first order convergence, the numerical method for solving space fractional partial differential equation constructed by using Hadamard finitepart integral approach does not require the solution of the equation satisfies the Dirichlet homogeneous boundary condition. Numerical results show that the experimentally determined convergence order obtained by using the Hadamard finitepart integral approach for solving space fractional partial differential equation with nonhomogeneous Dirichlet boundary conditions is indeed higher than the convergence order obtained by using the numerical methods constructed with the standard shifted Gr\"unwaldLetnikov formula or Lubich's higer order approximation schemes.

Numerical Approximation of Stochastic TimeFractional DiffusionWe develop and analyze a numerical method for stochastic timefractional diffusion driven by additive fractionally integrated Gaussian noise. The model involves two nonlocal terms in time, i.e., a Caputo fractional derivative of order $\alpha\in(0,1)$, and fractionally integrated Gaussian noise (with a RiemannLiouville fractional integral of order $\gamma \in[0,1]$ in the front). The numerical scheme approximates the model in space by the standard Galerkin method with continuous piecewise linear finite elements and in time by the classical Gr\"unwaldLetnikov method, and the noise by the $L^2$projection. Sharp strong and weak convergence rates are established, using suitable nonsmooth data error estimates for the deterministic counterpart. One and twodimensional numerical results are presented to support the theoretical findings.

A discrete mutualism model: analysis and exploration of a financial applicationWe perform a stability analysis on a discrete analogue of a known, continuous model of mutualism. We illustrate how the introduction of delays affects the asymptotic stability of the system’s positive nontrivial equilibrium point. In the second part of the paper we explore the insights that the model can provide when it is used in relation to interacting financial markets. We also note the limitations of such an approach.

An Altered Four Circulant Construction for SelfDual Codes from Group Rings and New Extremal Binary Selfdual Codes IWe introduce an altered version of the four circulant construction over group rings for selfdual codes. We consider this construction over the binary field, the rings F2 + uF2 and F4 + uF4; using groups of order 4 and 8. Through these constructions and their extensions, we find binary selfdual codes of lengths 16, 32, 48, 64 and 68, many of which are extremal. In particular, we find forty new extremal binary selfdual codes of length 68, including twelve new codes with \gamma=5 in W68,2, which is the first instance of such a value in the literature.

Characteristic functions of differential equations with deviating argumentsThe material here is motivated by the discussion of solutions of linear homogeneous and autonomous differential equations with deviating arguments. If $a, b, c$ and $\{\check{\tau}_\ell\}$ are real and ${\gamma}_\natural$ is realvalued and continuous, an example with these parameters is \begin{equation} u'(t) = \big\{a u(t) + b u(t+\check{\tau}_1) + c u(t+\check{\tau}_2) \big\} { \red +} \int_{\check{\tau}_3}^{\check{\tau}_4} {{\gamma}_\natural}(s) u(t+s) ds \tag{\hbox{$\rd{\star}$}} . \end{equation} A wide class of equations ($\rd{\star}$), or of similar type, can be written in the {\lq\lq}canonical{\rq\rq} form \begin{equation} u'(t) =\DSS \int_{\tau_{\rd \min}}^{\tau_{\rd \max}} u(t+s) d\sigma(s) \quad (t \in \Rset), \hbox{ for a suitable choice of } {\tau_{\rd \min}}, {\tau_{\rd \max}} \tag{\hbox{${\rd \star\star}$}} \end{equation} where $\sigma$ is of bounded variation and the integral is a RiemannStieltjes integral. For equations written in the form (${\rd{\star\star}}$), there is a corresponding characteristic function \begin{equation} \chi(\zeta) ):= \zeta  \DSS \int_{\tau_{\rd \min}}^{\tau_{\rd \max}} \exp(\zeta s) d\sigma(s) \quad (\zeta \in \Cset), \tag{\hbox{${\rd{\star\star\star}}$}} \end{equation} %%($ \chi(\zeta) \equiv \chi_\sigma (\zeta)$) whose zeros (if one considers appropriate subsets of equations (${\rd \star\star}$)  the literature provides additional information on the subsets to which we refer) play a r\^ole in the study of oscillatory or nonoscillatory solutions, or of bounded or unbounded solutions. We show that the related discussion of the zeros of $\chi$ is facilitated by observing and exploiting some simple and fundamental properties of characteristic functions.

Quadruple Bordered Constructions of SelfDual Codes from Group RingsIn this paper, we introduce a new bordered construction for selfdual codes using group rings. We consider constructions over the binary field, the family of rings Rk and the ring F4 + uF4. We use groups of order 4, 12 and 20. We construct some extremal selfdual codes and nonextremal selfdual codes of length 16, 32, 48, 64 and 68. In particular, we construct 33 new extremal selfdual codes of length 68.

Analysis of transient RivlinEricksen fluid and irreversibility of exothermic reactive hydromagnetic variable viscosityThe study analysed unsteady RivlinEricksen fluid and irreversibility of exponentially temperature dependent variable viscosity of hydromagnetic twostep exothermic chemical reactive flow along the channel axis with walls convective cooling. The nonNewtonian HeleShaw flow of RivlinErickson fluid is driven by bimolecular chemical kinetic and unvarying pressure gradient. The reactive fluid is induced by periodic changes in magnetic field and time. The Newtons law of cooling is satisfied by the constant heat coolant convection exchange at the wall surfaces with the neighboring regime. The dimensionless nonNewtonian reactive fluid equations are numerically solved using a convergent and consistence semiimplicit finite difference technique which are confirmed stable. The response of the reactive fluid flow to variational increase in the values of some entrenched fluid parameters in the momentum and energy balance equations are obtained. A satisfying equations for the ratio of irreversibility, entropy generation and Bejan number are solved with the results presented graphically and discussed quantitatively. From the study, it was obtained that the thermal criticality conditions with the right combination of thermofluid parameters, the thermal runaway can be prevented. Also, the entropy generation can minimize by at low dissipation rate and viscosity.

SpaceTime Discontinuous Galerkin Methods for the '\eps'dependent Stochastic AllenCahn Equation with mild noiseWe consider the $\eps$dependent stochastic AllenCahn equation with mild space time noise posed on a bounded domain of R^2. The positive parameter $\eps$ is a measure for the inner layers width that are generated during evolution. This equation, when the noise depends only on time, has been proposed by Funaki in [15]. The noise although smooth becomes white on the sharp interface limit as $\eps$ tends to zero. We construct a nonlinear dG scheme with spacetime finite elements of general type which are discontinuous in time. Existence of a unique discrete solution is proven by application of Brouwer's Theorem. We first derive abstract error estimates and then for the case of piecewise polynomial finite elements we prove an error in expectation of optimal order. All the appearing constants are estimated in terms of the parameter $\eps$. Finally, we present a linear approximation of the nonlinear scheme for which we prove existence of solution and optimal error in expectation in piecewise linear finite element spaces. The novelty of this work is based on the use of a finite element formulation in space and in time in 2+1dimensional subdomains for a nonlinear parabolic problem. In addition, this problem involves noise. These type of schemes avoid any RungeKutta type discretization for the evolutionary variable and seem to be very effective when applied to equations of such a difficulty.

Numerical analysis of a twoparameter fractional telegraph equationIn this paper we consider the twoparameter fractional telegraph equation of the form $$\, ^CD_{t_0^+}^{\alpha+1} u(t,x) + \, ^CD_{x_0^+}^{\beta+1} u (t,x) \, ^CD_{t_0^+}^{\alpha}u (t,x)u(t,x)=0.$$ Here $\, ^CD_{t_0^+}^{\alpha}$, $\, ^CD_{t_0^+}^{\alpha+1}$, $\, ^CD_{x_0^+}^{\beta+1}$ are operators of the Caputotype fractional derivative, where $0\leq \alpha < 1$ and $0 \leq \beta < 1$. The existence and uniqueness of the equations are proved by using the Banach fixed point theorem. A numerical method is introduced to solve this fractional telegraph equation and stability conditions for the numerical method are obtained. Numerical examples are given in the final section of the paper.

Stabilizing a mathematical model of plant species interactionIn this paper, we will consider how to stabilize a mathematical model of plant species interaction which is modelled by using LotkaVolterra system. We first identify the unstable steady states of the system, then we use the feedback control based on the solutions of the Riccati equation to stabilize the linearized system. We further stabilize the nonlinear system by using the feedback controller obtained in the stabilization of the linearized system. We introduce the backward Euler method to approximate the feedback control nonlinear system and obtain the error estimates. Four numerical examples are given which come from the application areas.

Stability of a numerical method for a fractional telegraph equationIn this paper, we introduce a numerical method for solving the timespace fractional telegraph equations. The numerical method is based on a quadrature formula approach and a stability condition for the numerical method is obtained. Two numerical examples are given and the stability regions are plotted.

On hereditary reducibility of 2monomial matrices over commutative ringsA 2monomial matrix over a commutative ring $R$ is by definition any matrix of the form $M(t,k,n)=\Phi\left(\begin{smallmatrix}I_k&0\\0&tI_{nk}\end{smallmatrix}\right)$, $0<k<n$, where $t$ is a noninvertible element of $R$, $\Phi$ the compa\nion matrix to $\lambda^n1$ and $I_k$ the identity $k\times k$matrix. In this paper we introduce the notion of hereditary reducibility (for these matrices) and indicate one general condition of the introduced reducibility.

A high order numerical method for solving nonlinear fractional differential equation with nonuniform meshesWe introduce a highorder numerical method for solving nonlinear fractional differential equation with nonuniform meshes. We first transform the fractional nonlinear differential equation into the equivalent Volterra integral equation. Then we approximate the integral by using the quadratic interpolation polynomials. On the first subinterval $[t_{0}, t_{1}]$, we approximate the integral with the quadratic interpolation polynomials defined on the nodes $t_{0}, t_{1}, t_{2}$ and in the other subinterval $[t_{j}, t_{j+1}], j=1, 2, \dots N1$, we approximate the integral with the quadratic interpolation polynomials defined on the nodes $t_{j1}, t_{j}, t_{j+1}$. A highorder numerical method is obtained. Then we apply this numerical method with the nonuniform meshes with the step size $\tau_{j}= t_{j+1} t_{j}= (j+1) \mu$ where $\mu= \frac{2T}{N (N+1)}$. Numerical results show that this method with the nonuniform meshes has the higher convergence order than the standard numerical methods obtained by using the rectangle and the trapzoid rules with the same nonuniform meshes.

New SelfDual and Formally SelfDual Codes from Group Ring ConstructionsIn this work, we study construction methods for selfdual and formally selfdual codes from group rings, arising from the cyclic group, the dihedral group, the dicyclic group and the semidihedral group. Using these constructions over the rings $_F2 +uF_2$ and $F_4 + uF_4$, we obtain 9 new extremal binary selfdual codes of length 68 and 25 even formally selfdual codes with parameters [72,36,14].