Search

The deadly coronavirus disease 2019 (COVID-19) has recently affected each corner of the world. Many governments of different countries have imposed strict measures in order to reduce the severity of the infection. In this present paper, we will study a mathematical model describing COVID-19 dynamics taking into account the government action and the individuals reaction. To this end, we will suggest a system of seven fractional deferential equations (FDEs) that describe the interaction between the classical susceptible, exposed, infectious, and removed (SEIR) individuals along with the government action and individual reaction involvement. Both human-to-human and zoonotic transmissions are considered in the model. The well-posedness of the FDEs model is established in terms of existence, positivity, and boundedness. The basic reproduction number (BRN) is found via the new generation matrix method. Different numerical simulations were carried out by taking into account real reported data from Wuhan, China. It was shown that the governmental action and the individuals' risk awareness reduce effectively the infection spread. Moreover, it was established that with the fractional derivative, the infection converges more quickly to its steady state.

We formulate a simple susceptible‐infectious‐recovery (SIR) model to describe the spread of the coronavirus under strict social restrictions. The transmission rate in this model is exponentially decreasing with time. We find a formula for basic reproduction function and estimate the maximum number of daily infected individuals. We fit the model to induced death data in Italy, United States, Germany, France, India, Spain, and China over the period from the first reported death to August 7, 2020. We notice that the model has excellent fit to the disease death data in these countries. We estimate the model's parameters in each of these countries with 95% confidence intervals. We order the strength of social restrictions in these countries using the exponential rate. We estimate the time needed to reduce the basic reproduction function to one unit and use it to order the quality of social restrictions in these countries. The social restriction in China was the strictest and the most effective and in India was the weakest and the least effective. Policy‐makers may apply the Chinese successful social restriction experiment and avoid the Indian unsuccessful one. [ FROM AUTHOR] Copyright of Mathematical Methods in the Applied Sciences is the property of John Wiley & Sons, Inc. and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full . (Copyright applies to all s.)

In the current article, we aim to study in detail a novel coronavirus (2019-nCoV or COVID-19) mathematical model for different aspects under Caputo fractional derivative. First, from analysis point of view, existence is necessary to be investigated for any applied problem. Therefore, we used fixed point theorem's due to Banach's and Schaefer's to establish some sufficient results regarding existence and uniqueness of the solution to the proposed model. On the other hand, stability is important in respect of approximate solution, so we have developed condition sufficient for the stability of Ulam-Hyers and their different types for the considered system. In addition, the model has also been considered for semianalytical solution via Laplace Adomian decomposition method (LADM). On Matlab, by taking some real data about Pakistan, we graph the obtained results. In the last of the manuscript, a detail discussion and brief conclusion are provided.

COVID-19 pandemic has affected all aspects of people's lives and disrupted the economy. Forecasting the number of cases infected with this virus can help authorities make accurate decisions on the interventions that must be implemented to control the pandemic. Investigation of the studies on COVID-19 forecasting indicates that various techniques such as statistical, mathematical, and machine and deep learning have been utilized. Although deep learning models have shown promising results in this context, their performance can be improved using auxiliary features. Therefore, in this study, we propose two hybrid deep learning methods that utilize the statistical features as auxiliary inputs and associate them with their main input. Specifically, we design a hybrid method of the multihead attention mechanism and the statistical features (ATT_FE) and a combined method of convolutional neural network and the statistical features (CNN_FE) and apply them to COVID-19 data of 10 countries with the highest number of confirmed cases. The results of experiments indicate that the hybrid models outperform their conventional counterparts in terms of performance measures. The experiments also demonstrate the superiority of the hybrid ATT_FE method over the long short-term memory model.

In this paper, we consider the propagation of waves in an open waveguide in ℝ$^{2}$ where the index of refraction is a local perturbation of a function which is periodic along the axis of the waveguide (which we choose to be the x$_{1}$ axis) and equal to one for |x$_{2}$| > h$_{0}$ for some h$_{0}$ > 0. Motivated by the limiting absorption principle (proven in an earlier paper by the author), we formulate a radiation condition which allows the existence of propagating modes and prove uniqueness, existence, and stability of a solution under the assumption that no bound states exist. In the second part, we determine the order of decay of the radiating part of the solution in the direction of the layer and in the direction orthogonal to it. Finally, we show that it satisfies the classical Sommerfeld radiation condition and allows the definition of a far field pattern.

In the paper mentioned in the title, it is proved the boundedness of the Riesz potential operator of variable order 𝛼(x) from variable exponent Morrey space to variable exponent Campanato space, under certain assumptions on the variable exponents p(x) and 𝜆(x) of the Morrey space. Assumptions on the exponents were different depending on whether 𝛼(x)p(x)−n+𝜆(x) p(x) takes or not the critical values 0 or 1. In this note, we improve those results by unifying all the cases and covering the whole range 0 ⩽ 𝛼(x)p(x)−n+𝜆(x) p(x) ⩽ 1. We also provide a correction to some minor technicality in the proof of Theorem 2 in the aforementioned paper. info:eu-repo/semantics/publishedVersion

The present paper studies a fractional-order SEIR epidemic model for the transmission dynamics of infectious diseases such as HIV and HBV that spreads in the host population. The total host population is considered bounded, and Holling type-II saturation incidence rate is involved as the infection term. Using the proposed SEIR epidemic model, the threshold quantity, namely basic reproduction number R0, is obtained that determines the status of the disease, whether it dies out or persists in the whole population. The model’s analysis shows that two equilibria exist, namely, disease-free equilibrium (DFE) and endemic equilibrium (EE). The global stability of the equilibria is determined using a Lyapunov functional approach. The disease status can be verified based on obtained threshold quantity R0. If R0 < 1, then DFE is globally stable, leading to eradicating the population’s disease. If R0 > 1, a unique EE exists, and that is globally stable under certain conditions in the feasible region. The Caputo type fractional derivative is taken as the fractional operator. The bifurcation and sensitivity analyses are also performed for the proposed model that determines the relative importance of the parameters into disease transmission. The numerical solution of the model is obtained by the generalized Adams- Bashforth-Moulton method. Finally, numerical simulations are performed to illustrate and verify the analytical results.

The paper formulates the nonlinear problem of steady-state heat conduction at the constant electric potential difference on the surfaces of a plane dielectric layer with the temperature-dependent heat conduction coefficient and electrical resistivity. A fixed temperature value is set on one of the layer surfaces, and the convective heat exchange with the ambient medium occurs on the opposite surface. The formulation of the problem is transformed into integral ratios, which allows the calculation of the temperature distribution over the layer thickness, governed both by the monotonic and nonmonotonic function. The quantitative assay of the temperature state of a layer of a polymer dielectric made of amorphous polycarbonate is given as an example, as well as the analysis of nonuniformity of the absolute value of electric field intensity over the thickness of this layer.

We study the process of two-phase flow in thin porous media domains of Brinkman-type. This is generally described by a model of coupled, mixed-type differential equations of fluids' saturation and pressure. To reduce the model complexity, different approaches that utilize the thin geometry of the domain have been suggested. We focus on a reduced model that is formulated as a single nonlocal evolution equation of saturation. It is derived by applying standard asymptotic analysis to the dimensionless coupled model, however, a rigid mathematical derivation is still lacking. In this paper, we prove that the reduced model is the analytical limit of the coupled two-phase flow model as the geometrical parameter of domain's width--length ratio tends to zero. Precisely, we prove the convergence of weak solutions for the coupled model to a weak solution for the reduced model as the geometrical parameter approaches zero.

This paper deals with the classification of transition phenomena in the most basic dissipative system possible, namely, the 1D reaction-diffusion equation. The emphasis is on the relation between the linear and nonlinear terms and the effect of the boundaries which influence the first transitions. We consider the cases where the linear part is self-adjoint with second-order and fourth-order derivatives which is the case which most often arises in applications. We assume that the nonlinear term depends on the unknown function and its first derivative which is basically the semilinear case for the second-order reaction-diffusion system. As for the boundary conditions, we consider the typical Dirichlet, Neumann, and periodic boundary settings. In all the cases, the equations admit a trivial steady state which loses stability at a critical parameter. We aim to classify all possible transitions and bifurcations that take place. Our analysis shows that these systems display all three types of transitions: continuous, jump and mixed. Moreover they exhibit transcritical, supercritical bifurcations with bifurcated states such as finitely many equilibria, circle of equilibria, and slowly rotating limit cycle. Many applications found in the literature are basically corollaries of our main results. We apply our results to classify the first transitions of the Chaffee-Infante equation, the Fisher-KPP equation, the Kuramoto-Sivashinsky equation, and the Swift-Hohenberg equation.

The purpose of this paper is to analyze a new kind of Hadamard fractional boundary value problem combining integral boundary condition and multipoint fractional integral boundary condition on an infinite interval. By the help of the Bai-Ge’s fixed point theorem, multiplicity results of positive solutions are derived for the Hadamard fractional boundary value problem. In the end, to illustrative the main result, an example is also presented.

In this paper, we prove two identities involving quantum derivatives, quantum integrals, and certain parameters. Using the newly proved identities, we prove new inequalities of Simpson's and Newton's type for quantum differentiable convex functions under certain assumptions. Moreover, we discuss the special cases of our main results and obtain some new and existing Simpson's type inequalities, Newton's type inequalities, midpoint type inequalities, and trapezoidal type inequalities. National Natural Science Foundation of ChinaNational Natural Science Foundation of China (NSFC) [11971241] National Natural Science Foundation of China, Grant/Award Number: 11971241 WOS:000710085000001 2-s2.0-85117687946

In this paper, we develop a higher order numerical method for the fractional Bagley-Torvik equation. The main tools used include a new fourth-order approximation for the fractional derivative based on the weighted shifted Grünwald-Letnikov difference operator and a discrete cubic spline approach. We show that the theoretical convergence order improves those of previous work. Five examples are further presented to illustrate the efficiency of our method and to compare with other methods in the literature.

This paper is devoted to the mathematical modeling of a combined effect of directional and bulk crystallization in a phase transition layer with allowance for nucleation and evolution of newly born particles. We consider two models with and without fluctuations in crystal growth velocities, which are analytically solved using the saddle-point technique. The particle-size distribution function, solid-phase fraction in a supercooled two-phase layer, its thickness and permeability, solidification velocity, and desupercooling kinetics are defined. This solution enables us to characterize the mushy layer composition. We show that the region adjacent to the liquid phase is almost free of crystals and has a constant temperature gradient. Crystals undergo intense growth leading to fast mushy layer desupercooling in the middle of a two-phase region. The mushy region adjacent to the solid material is filled with the growing solid-phase structures and is almost desupercooled. © 2021 The Authors. Mathematical Methods in the Applied Sciences published by John Wiley & Sons, Ltd. Russian Science Foundation, RSF: 21-79-10012 The authors gratefully acknowledge financial support from the Russian Science Foundation (project no. 21-79-10012). Open Access funding enabled and organized by Projekt DEAL.

This paper is devoted to study the discrete Sturm-Liouville problem $$ \left\{\begin{array}{ll} -\Delta(p(k)\Delta u(k-1))+q(k)u(k)=\lambda m(k)u(k)+f_1(k,u(k),\lambda)+f_2(k,u(k),\lambda),\ \ k\in[1,T]_Z,\\[2ex] a_0u(0)+b_0\Delta u(0)=0,\ a_1u(T)+b_1\Delta u(T)=0, \end{array}\right. $$ where $\lambda\in\mathbb{R}$ is a parameter, $f_1, f_2\in C([1,T]_Z\times\mathbb{R}^2, \mathbb{R})$, $f_1$ is not differentiable at the origin and infinity. Under some suitable assumptions on nonlinear terms, we prove the existence of unbounded continua of positive and negative solutions of this problem which bifurcate from intervals of the line of trivial solutions or from infinity, respectively.

In this paper, we deal with a new approximation of a diffusive predator--prey model with Leslie--Gower term and social behavior for the prey subject to Neumann boundary conditions. A new approach for a predator-prey interaction in the presence of prey social behavior has been considered. Our main topic in this work is to study the influence of the prey's herd shape on the predator-prey interaction in the presence of Leslie--Gower term. First of all, we examine briefly the system without spatial diffusion. By analyzing the distribution of the eigenvalues associated with the constant equilibria, the local stability of the equilibrium points and the existence of Hopf bifurcation have been investigated. Then, the spatiotemporal dynamics introduced by self diffusion was determined, where the existence of the positive solution, Hopf bifurcation, Turing driven instability, Turing-Hopf bifurcation point have been derived. Further, the effect of the prey's herd shape rate on the prey and predator equilibrium densities as well as on the Hopf bifurcating points has been discussed. Finally, by using the normal form theory on the center manifold, the direction and stability of the bifurcating periodic solutions have also been obtained. To illustrate the theoretical results, some graphical representations are given.

In recent years various probability models have been proposed for describing lifetime data. Increasing model flexibility is often sought as a means to better describe asymmetric and heavy tail distributions. Such extensions were pioneered by the beta-G family. However, efficient goodness-of-fit (GoF) measures for the beta-G distributions are sought. In this paper, we combine probability weighted moments (PWMs) and the Mellin transform (MT) in order to furnish new qualitative and quantitative GoF tools for model selection within the beta-G class. We derive PWMs for the Fr\’{e}chet and Kumaraswamy distributions; and we provide expressions for the MT, and for the log-cumulants (LC) of the beta-Weibull, beta-Fr\’{e}chet, beta-Kumaraswamy, and beta-log-logistic distributions. Subsequently, we construct LC diagrams and, based on the Hotelling’s $T^2$ statistic, we derive confidence ellipses for the LCs. Finally, the proposed GoF measures are applied on five real data sets in order to demonstrate their applicability.

This paper studies the inhomogeneous nonlinear Schrödinger equations, which may model the propagation of laser beams in nonlinear optics. Using the cross-constrained variational method, a sharp condition for global existence is derived. Then, by solving a variational problem, the strong instability of solitary waves of this equation is proved.

In this paper, we first prove three identities for twice partially differentiable functions. Then, by using these equalities, we obtain several trapezoid and Ostrowski type inequalities via generalized fractional integrals for functions with bounded partial derivatives. Moreover, we present some results for Riemann-Liouville fractional integrals by special choice of main results. WOS:000687664600001 2-s2.0-85113769051

The aim of this paper is to give an alternative technique for the derivation of a prior energy estimate. Consequently, this allows to define a natural energy norm of the long-term well-posedness result established by S. Israwi in [2] but for the original system, in which the partial operator ∇× is not involved.

In this paper, we study a fractional order time delay for nonlinear financial system. By using Laplace transformation, stability and Hopf bifurcation analysis have been done for the model. Furthermore, numerical simulation has been carried out for better understanding of our results.

This paper describes the foundations of a differential geometry of a quaternionic curves. The Frenet-Serret equations and the evolutes and evolvents of a particular quaternionic curve are accordingly determined. This new formulation takes benefit of the quaternionic structure and the results are much simpler than the present formulations of quaternionic curves., Accept by Mathematical Methods in the Applied Sciences

In this paper, a model that was recently derived in Reinken et al. [11] to describe the dynamics of microswimmer suspensions is studied. In particular, the global existence of weak solutions, their weak-strong uniqueness and a connection to a different model that was proposed in Wensink et al. [18] is shown., 18 pages

The first case of COVID-19 in India detected on January 30, 2020, after its emergence in Wuhan, China, in December 2019. The lockdown was imposed as anemergency measure by the Indian government to prevent the spread of COVID-19 but gradually eased out due to its vast economic consequences. Just 15?days after the relaxation of lockdown restrictions, Delhi became India's worst city in terms of COVID-19 cases. In this paper, we propose a variable-order fractional SIR (susceptible, infected, removed) model at state-level scale. We introduce a algorithm that uses the differential evolution algorithm in combination with Adam?Bashforth?Moulton method to learn the parameters in a system of variable-order fractional SIR model. The model can predict the confirm COVID-19 cases in India considering the effects of nationwide lockdown and the possible estimate of the number of infliction inactive cases after the removal of lockdown on June 1, 2020. A new parameter p is introduced in the classical SIR model representing the fraction of infected people that get tested and are thereby quarantined. The COVID-19 trajectory in Delhi, as per our model, predicts the slowing down of the spread between January and February 2021, touching a peak of around 5 lakh confirmed cases.

In this paper, we investigate the wellposedness of the non-isentropic compressible Navier-Stokes/Allen-Cahn system with the heat-conductivity proportional to a positive power of the temperature. This system describes the flow of a two-phase immiscible heat-conducting viscous compressible mixture. The phases are allowed to shrink or grow due to changes of density in the fluid and incorporates their transport with the current. We established the global existence and uniqueness of strong solutions for this system in 1-D, which means no phase separation, vacuum, shock wave, mass or heat or phase concentration will be developed in finite time, although the motion of the two-phase immiscible flow has large oscillations and the interaction between the hydrodynamic and phase-field effects is complex. Our result can be regarded as a natural generalization of the Kazhikhov-Shelukhin's result ([Kazhikhov-Shelukhin. J. Appl. Math. Mech. 41 (1977)]) for the compressible single-phase flow with constant heat conductivity to the non-isentropic compressible immiscible two-phase flow with degenerate and nonlinear heat conductivity.

In this paper, we investigate the numerical approximation of stochastic convection-reaction-diffusion equations using two explicit exponential integrators. The stochastic partial differential equation (SPDE) is driven by additive Wiener process. The approximation in space is done via a combination of the standard finite element method and the Galerkin projection method. Using the linear functional of the noise, we construct two accelerated numerical methods, which achieve higher convergence orders. In particular, we achieve convergence rates approximately $1$ for trace class noise and $\frac{1}{2}$ for space-time white noise. These convergences orders are obtained under less regularities assumptions on the nonlinear drift function than those used in the literature for stochastic reaction-diffusion equations. Numerical experiments to illustrate our theoretical results are provided

In this paper, we introduce the concepts of asymptotically f-statistical equivalence, asymptotically f-lacunary statistical equivalence, and strong asymptotically f-lacunary equivalence for non-negative two delta measurable real-valued functions defined on time scales with the aid of modulus function f. Furthermore, the relationships between these new concepts are investigated. We also present some inclusion theorems.

One of the main motivation points in studies on inequalities is to obtain generalizations and to introduce new approaches. In this direction, the generalized fractional integral operators defined within the scope of fractional analysis are quite functional. In this paper, some new integral inequalities have been proved by using generalized fractional integral operators and some classical inequalities for integrable functions. In the proofs of the main findings, the definitions of the generalized fractional integral operator, certain classical relations, and some classical inequalities are used., H.E.C. Pakistan under NRPU [7906]; H.E.C. Pakistan [7906], H.E.C. Pakistan under NRPU project 7906; H.E.C. Pakistan, Grant/Award Number: 7906

The aim of this paper is to construct a numerical scheme for solving a class of bilateral free boundaries problem. First, using a shape functional and some regularization terms, an optimal control problem is formulated, in addition, we prove its solution existence's. The first optimality conditions and the shape gradient are computed. the proposed numerical scheme is a genetic algorithm guided conjugate gradient combined with the finite element method, at each mesh regeneration, we perform a mesh refinement in order to avoid any domain singularities. Some numerical examples are shown to demonstrate the validity of the theoretical results, and to prove the robustness of the proposed scheme.

In this paper, we further Meirong Zhang, et al.'s work by computing the number of weighted eigenvalues for Sturm-Liouville equations, equipped with general integrable potentials and Dirac weights, under Dirichlet boundary condition. We show that, for a Sturm-Liouville equation with a general integrable potential, if its weight is a positive linear combination of $n$ Dirac Delta functions, then it has at most $n$ (may be less than $n$, or even be $0$) distinct real Dirichlet eigenvalues, or every complex number is a Dirichlet eigenvalue; in particular, under some sharp condition, the number of Dirichlet eigenvalues is exactly $n$. Our main method is to introduce the concepts of characteristics matrix and characteristics polynomial for Sturm-Liouville problem with Dirac weights, and put forward a general and direct algorithm used for computing eigenvalues. As an application, a class of inverse Dirichelt problems for Sturm-Liouville equations involving single Dirac distribution weights is studied., 23 pages

[EN] This paper presents a methodology to deal with the randomness associated to ecological modelling. Data variability makes it necessary to analyse the impact of random perturbations on the fitted model parameters. We conduct such analysis for the logistic growth model with a certain sigmoid functional form of the carrying capacity, which was proposed in the literature for the study of parasite growth during infection. We show how the probability distributions of the parameters are set via the maximum entropy principle. Then the random variable transformation method allows for computing the density function of the population., Spanish Ministerio de Economia y Competitividad, Grant/Award Number: PID2020-115270GB-I00; Agencia Estatal de Investigacion

In this paper, a procedure for the detection of the sources of industrial noise and the evaluation of their distances is introduced. The above method is based on the analysis of acoustic and optical data recorded by an acoustic camera. In order to improve the resolution of the data, interpolation and quasi interpolation algorithms for digital data processing have been used, such as the bilinear, bicubic, and sampling Kantorovich (SK). The experimental tests show that the SK algorithm allows to perform the above task more accurately than the other considered methods.

This paper is concerned with a parabolic-elliptic Keller-Segel system where both diffusive and chemotactic coefficients (motility functions) depend on the chemical signal density. This system was originally proposed by Keller and Segel in \cite{KS-1971-JTB2} to describe the aggregation phase of {\it Dictyostelium discoideum} cells in response to the secreted chemical signal cyclic adenosine monophosphate (cAMP), but the available analytical results are very limited by far. Considering system in a bounded smooth domain with Neumann boundary conditions, we establish the global boundedness of solutions in any dimensions with suitable general conditions on the signal-dependent motility functions, which are applicable to a wide class of motility functions. The existence/nonexistence of non-constant steady states is studied and abundant stationary profiles are found. Some open questions are outlined for further pursues. Our results demonstrate that the global boundedness and profile of stationary solutions to the Keller-Segel system with signal-dependent motilities depend on the decay rates of motility functions, space dimensions and the relation between the diffusive and chemotactic motilities, which makes the dynamics immensely wealthy.

The purpose of this paper is to give the delay control based on discrete-time state observations to stabilize highly nonlinear hybrid stochastic functional differential equations (SFDEs). It is considered that time lag generated by the controller in each discrete observation should be different. The new controlled hybrid SFDEs are affected the variable delay caused by the controller, the distributed delay and the superlinear coefficients of the systems itself, which makes the problem handling more complicated. Then, a series of criteria for the exponential stability of the controlled SFDEs are obtained, and an upper bound for the discrete observation interval and variable delay is given. Finally, numerical example illustrate the proposed theoretical results.

Massless field equations are fundamental in particle physics. In Clifford analysis, the Euclidean version of these equations has been dealt with but it is not clear, even in dimension 4, what should be the right analogue of massless field equations for fields with values in a general irreducible Spin(4)-module. The main aim of the paper is to explain that a good possibility is to take the so-called generalized Cauchy-Riemann equations proposed a long time ago by E. Stein and G. Weiss. For this choice of the equations, we show that their polynomial solutions form different irreducible Spin(4)-modules. This is an important step in developing the corrresponding function theory.

In this paper, the application of Ledley antecedence solutions in designing state feedback stabilizers of generic logic systems has been proposed. To make the method feasible, two modifications are made to the original Ledley antecedence solution theory: (i) the preassigned logical functions have been extended from being a set of equations to an admissible set; (ii) the domain of arguments has been extended from the whole state space to a restricted subset. In the proposed method, state feedback controls are considered as a set of extended Ledley antecedence solutions for a designed iterative admissible sets over their corresponding restricted subsets. Based on this, an algorithm has been proposed to verify the solvability, and simultaneously to provide all possible state feedback stabilizers when the problem is solvable. All stabilizers are optimal, which stabilize the logic systems from any initial state to the destination state/state set in the shortest time. The method is firstly demonstrated on Boolean control networks to achieve point stabilization. Then, with some minor modifications, the proposed method is also proven to be applicable to set stabilization problems. Finally, it is shown that in $k$-valued and mix-valued logical systems, the proposed method remains effective., Comment: 10 pages, 3 figures

© 2021 The Authors. The aim of this paper is to develop a layer potential theory in L2-based weighted Sobolev spaces on Lipschitz bounded and exterior domains of Rn , n ≥ 3, for the anisotropic Stokes system with L∞ viscosity tensor coefficient satisfying an ellip- ticity condition for symmetric matrices with zero matrix trace. To do this, we explore equivalent mixed variational formulations and prove the well-posedness of some transmission problems for the anisotropic Stokes system in Lipschitz domains of Rn, with the given data in L2-based weighted Sobolev spaces. These results are used to define the volume (Newtonian) and layer potentials and to obtain their properties. Then, we analyze the well-posedness of the exterior Dirichlet and Neumann problems for the anisotropic Stokes system with L∞ symmetrically elliptic tensor coefficient by representing their solutions in terms of the obtained volume and layer potentials. EPSRC grant EP/M013545/1: "Mathematical Analysis of Boundary-Domain Integral Equations for Nonlinear PDEs"; Babeş-Bolyai University research grant AGC35124/31.10.2018; Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy-EXC 2075-390740016.