## On solutions related to FitzHugh-Rinzel type model

Fabio De Angelis ∗ Monica De Angelis †

## Abstract

A ternary autonomous dynamical system of FitzHugh-Rinzel type is analyzed. The system, at start, is reduced to a nonlinear integro differential equation. The fundamental solution H x,t ( ) is explicitly determined and the initial value problem is analyzed in the whole space. The solution is expressed by means of an integral equation involving H x,t ( ). Moreover, adding an extra control term, explicit solutions are achieved.

## 1 Introduction

The FitzHugh-Rinzel (FHR) system [1-4] is a three-dimensional model deriving from the FitzHugh-Nagumo (FHN) model [5-12] developed to incorporate bursting phenomena of nerve cells. Indeed, a number of different cell types exhibit a behaviour characterized by brief bursts of oscillatory activity alternated by quiescient periods during which the membrane potential only changes slowly, and this behaviour is called bursting, see e.g. [13]. Accordingly, bursting oscillations are characterized by a variable of the system that changes periodically from an active phase of rapid spike oscillations to a silent phase. These phenomena are becoming increasingly important as they are being investigated in many scientific fields. Indeed, phenomena of bursting have been observed as electrical behaviours in many nerve and endocrine cells such as hippocampal and thalamic neurons, mammalian midbrain and pancreatic in β -cells, see e.g. [1] and references therein. Also, in the cardiovascular system, bursting oscillations are generated by the electrical activity of cardiac cells that excite the heart membrane to produce the contraction of ventricles and auricles [14]. Furthermore, bursting oscillations represent a topic of potential interest in dynamics and bifurcation mechanisms of devices and structures and in the analysis of nonlinear problems in mechanics [15-19]. Recent studies proved that the development of this field helps in the studying of the restoration of synaptic connections.

∗ Dept. of Structures for Engineering and Architecture, University of Naples Federico II, Via Claudio 21, 80125, Naples, Italy.

† Dept. of Mathematics and Applications 'R. Caccioppoli', University of Naples Federico II, Via Cinthia 26,80126, Naples, Italy.

,

modeange@unina.it

Indeed, it seems that nanoscale memristor devices have potential to reproduce the behaviour of a biological synapse [20,21]. This would lead in the future, also in case of traumatic lesions, to the introduction of electronic synapses to connect neurons directly.

The paper is organized as follows. In section 1.1 the mathematical problem is defined and the state of the art with the aim of the paper are discussed. In section 2, the explicit expression of the fundamental solution H x,t ( ) is achieved and some properties are proved. In section 3 the integral solution for the initial value problem is given. In section 4 the insertion of an extra term allows us to obtain explicit solutions for the model.

## 1.1 Mathematical considerations, state of the art and aim of the paper

Generally, denoting by D,ε,β,c constant parameters, the (FHN) model is a p.d.e. system such that

$$\begin{cases} \begin{array} {$$

where function f ( u ) depends on the reaction kinetics of the model. In the literature f ( u ) can assume a piecewise linear form, see, e.g. [22] and reference therein, or f ( u ) = u -u / 3 3 [12]. However, in general, one has [5,13]:

$$( 1. 2 ) \quad \ \ f ( u ) = u \left ( \, a - u \, \right$$

As for the FitzHugh-Rinzel model, most of the articles consider the following system characterized by three o.d.e.:

$$( 1. 3 ) \quad & \begin{cases} \begin{array} { c } \frac { \partial u$$

where I ext , ε, β, c, d, h, β, δ indicate arbitrary constants.

In this paper, in order to evaluate the contribute of a diffusion term, the following FitzHugh-Rinzel type system is considered:

$$( 1. 4 ) \quad & \begin{cases} \ \frac { \partial \, u } { \partial \$$

Indeed, the second order term with D &gt; 0 represents just the diffusion contribution and it can be associated to the axial current in the axon. It derives from the Hodgkin-Huxley (HH) theory for nerve membranes where, if d represents the axon diameter and r i is the resistivity, the spatial variation in the potential V gives the term ( d/ r 4 i ) V xx from which term Du xx derives [23].

Moreover it is also assumed β, d, ε, δ as positive constants that together with c, h, characterize the model's kinetic.

Model (1.4) can be considered as a two time-scale slow-fast system with two fast variables ( u, w ) and one slow variable ( y ). However, if for instance, ε = δ the system can be considered as a two time-scale with one fast variable u and two slow variables ( w,y ). Otherwise, if δ and ε have significant difference, it can also be considered as a three-time-scale system with the fast variable u, the intermediate variable and the slow variable [24].

As for function f ( u ) one considers the non-linear form expressed in formula (1.2). As a consequence, it results

$$( 1. 5 ) \quad \ f ( u ) \, = \, - \, a \, u \, +$$

Then, the system (1.4) becomes

$$( 1. 6 ) \begin{cases} \ \frac { \partial u } { \partial t } = \ D \, \frac { \partial ^ { 2 } \, u } { \partial x ^ { 2 } } - a u \ - \ w \ + y \ + \varphi ( u ) \\ \ \frac { \partial w } { \partial t } \ = \ \varepsilon ( - \beta w + c + u ) \\ \ \frac { \partial y } { \partial t } \ = \ \delta ( - u + h - d y ). \end{cases}$$

Indicating by means of

$$( 1. 7 ) \quad \ u ( x, 0 ) \, = \, u _ { 0 } \,, \quad w ( x, 0 ) \, = \, w _ { 0 } \quad \ y ( x, 0 ) = y _ { 0 }, \quad \quad \ ( \, x \, \in \, \Re \, )$$

the initial values, from (1 6) . 2 3 , it follows:

$$\begin{cases} \ w \, = \, w _ { 0 } \, e ^ { - \varepsilon \beta t } \, + \, \frac { c } { \beta } \, ( 1 - e ^ { - \varepsilon \beta t } ) \, + \varepsilon \int _ { 0 } ^ { t } e ^ { - \varepsilon \beta \, ( t - \tau ) } \, u ( x, \tau ) \, d \tau \\ \ y \, = \, y _ { 0 } \, e ^ { - \delta \, d t } \, + \, \frac { h } { d } \, ( 1 - e ^ { - \delta \, d t } ) \, - \delta \int _ { 0 } ^ { t } e ^ { - \delta \, d \, ( t - \tau ) } \, u ( x, \tau ) \, d \tau. \\ \text{Constant} \, \text{$the$random$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$that$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$to$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$and$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$for$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$of$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$THE$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$their$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$other$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$or$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$The$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$th$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$get$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$both$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$he$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$
$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$there$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$them$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$the$$

Consequently, denoting the source term by

$$( 1. 9 ) \ F ( x, t, u ) = \varphi ( u ) - w _ { 0 } ( x ) e ^ { - \varepsilon \beta t } + y _ { 0 } ( x ) e ^ { - \delta d t } - \frac { c } { \beta } ( 1 - e ^ { - \varepsilon \beta t } ) + \frac { h } { d } ( 1 - e ^ { - \delta d t } ),$$

problem (1.6)-(1.7) can be modified into the following initial value problem P :

$$( 1. 1 0 ) \quad \begin{cases} \ u _ { t } - D u _ { x x } + a u + \int _ { 0 } ^ { t } [ \varepsilon e ^ { - \varepsilon \beta ( t - \tau ) } + \delta e ^ { - \delta d ( t - \tau ) } ] u ( x, \tau ) d \tau = F ( x, t, u ) \\ \ u ( x, 0 ) = u _ { 0 } ( x ) \quad x \in \Re. \end{cases}$$

As for the state of art, mathematical considerations allow to assert that the knowledge of the fundamental solution H x,t ( ) related to the linear parabolic operator L :

$$L u \equiv u _ { t } - D u _ { x x } + a u + \int _ { 0 } ^ { t } [ \varepsilon \, e ^ { - \varepsilon \beta ( t - \tau ) } \, + \delta e ^ { - \delta d ( t - \tau ) } ] u ( x, \tau ) \, d \tau,$$

leads to determine the solution of P . Indeed, if F ( x, t, u ) verifies appropriate assumptions, through the fixed point theorem, solution can be expressed by means of an integral equation, see f.i. [25,26].

Moreover, according to [26], when operator L assumes a similar but simpler form, many properties and inequalities are achieved.

The aim of the paper is to explicitly determine the fundamental solution H x,t ( ) which involves naturally the diffusion constant D. Then, the initial value problem in all the space is analyzed and the solution is deduced by means of an integral equation. Moreover, using a method of travelling wave, solutions of a modified FitzHughRinzel type system have been explicitly determined pointing out the influence of the diffusion parameter D.

## 2 Fundamental solution and its properties

Indicating by T an arbitrary positive constant, let us consider the initialvalue problem (1.10) defined in the whole space Ω T :

$$\Omega _ { T } = \{ ( x, t ) \, \colon x \in \Re, \ \ 0 < t \leq T \},$$

and let us denote by

$$\hat { u } ( x, s ) \, = \int _ { 0 } ^ { \infty } \, e ^ { - s t } \, u ( x, t ) \, d t \,, \ \hat { F } ( x, s ) \, = \int _ { 0 } ^ { \infty } \, e ^ { - s t } \, F \left [ x, t, u ( x, t ) \right ] d t \,,$$

the Laplace transform with respect to t. If ˆ ( H x,s ) expresses the L t transform of the fundamental solution H x,t , ( ) from (1.10) it follows:

$$\hat { u } ( x, s ) \, = \, \int _ { \Re } \, \hat { H } \left ( \, x - \xi, s \, \right ) \left [ \, u _ { 0 } ( \, \xi \, ) \, + \, \hat { F } ( \xi, s ) \, \right ] d \xi \,,$$

and formally it follows that

$$u ( x, t ) = \int _ { \Re } \, H ( x - \xi, t ) \, u _ { 0 } ( \xi ) \, d \xi \, +$$

(2.13)

$$+ \, \int _ { 0 } ^ { t } d \tau \int _ { \Re } H ( x - \xi, t - \tau ) \ F \left [ \xi, \tau, u ( \xi, \tau ) \right ] \, d \xi.$$

So that, denoting by J 1 ( z ) the Bessel function of first kind and order 1 , let us consider the following functions:

$$H _ { 1 } ( x, t ) \, = \, \frac { e ^ { - \frac { x ^ { 2 } } { 4 D t } } } { 2 \sqrt { \pi D t } } \ e ^ { - a \, t } \, +$$

(2.14)

$$\cdot \quad ^ { \prime } - \frac { 1 } { 2 } \int _ { 0 } ^ { t } \frac { e ^ { - \frac { x ^ { 2 } } { 4 D y } - a y } } { \sqrt { t - y } } \ \frac { \sqrt { \varepsilon } \ e ^ { - \beta \varepsilon ( t - y ) } } { \sqrt { \pi \ D } } J _ { 1 } ( \, 2 \, \sqrt { \varepsilon \, y ( t - y ) } \ ) \ \} d y,$$

$$( 2. 1 5 ) \quad \ \ H _ { 2 } = \int _ { 0 } ^ { t } H _ { 1 } ( x, y ) \ e ^ { - \delta d ( t - y ) } \sqrt { \frac { \delta y } { t - y } } J _ { 1 } ( \, 2 \, \sqrt { \delta \, y \, ( t - y ) } \ d y.$$

Besides, by setting

$$\sigma ^ { 2 } \ = \ s \, + \, a \, + \, \frac { \delta } { s + \delta d } \, + \, \frac { \varepsilon } { s + \, \beta \varepsilon },$$

and by denoting

$$( 2. 1 7 ) & & H = H _ { 1 } - H _ { 2 },$$

the following theorem holds:

Theorem 1. In the half-plane glyph[Rfractur] e s &gt; max ( -a, -βε, -δd ) the Laplace integral L t H converges absolutely for all x &gt; 0 , and it results:

$$\mathcal { L } _ { t } \, H \, \equiv \, \hat { H } \, = \int _ { 0 } ^ { \infty } e ^ { - s t } \, H \left ( x, t \right ) \, d t \, = \, \frac { 1 } { \sqrt { D } } \, \frac { e ^ { - \frac { | x | } { \sqrt { D } } \, \sigma } } { 2 \, \sigma }.$$

Moreover, function H(x,t) satisfies some properties typical of the fundamental solution of heat equation, such as:

- a) H x,t ( ) ∈ C ∞ , t &gt; 0 , x ∈ glyph[Rfractur] ,
- b) for fixed t &gt; 0 , H and its derivatives are vanishing esponentially fast as | x | tends to infinity.
- c) In addition, it results lim t → 0 H x,t ( ) = 0 , for any fixed η &gt; 0 , uniformly for all | x | ≥ η.

Proof. Since for all real z one has | J 1 ( z ) | ≤ 1 , the Fubini -Tonelli theorem assures that

$$L _ { t } & \left ( \frac { 1 } { 2 } \quad \int _ { 0 } ^ { t } \frac { e ^ { - \frac { x ^ { 2 } } { 4 D y } } - a y } { \sqrt { t - y } } \ \frac { \sqrt { \varepsilon } \ e ^ { - \beta \varepsilon } ( t - y ) } { \sqrt { \pi } \ D } J _ { 1 } ( 2 \sqrt { \varepsilon } \, y ( t - y ) \ \right ) d y = \\ & = \ \frac { \sqrt { \varepsilon } } { 2 \, \sqrt { \pi } D } \quad \int _ { 0 } ^ { \infty } e ^ { - ( s + a ) y - \frac { x ^ { 2 } } { 4 D y } } \, d y \ \int _ { 0 } ^ { \infty } e ^ { - ( s + \beta \varepsilon ) t } \ J _ { 1 } ( 2 \sqrt { \varepsilon } \, y t \, ) \ \frac { d t } { \sqrt { t } } \\ \text{and bein} \sigma$$

and being

$$( 2. 1 9 ) \quad \int _ { 0 } ^ { \infty } \, e ^ { - p t } \, \sqrt { \frac { c } { t } } \ J _ { 1 } ( 2 \sqrt { c t } \, ) \, d t = \, 1 \, - \, e ^ { - c / p } \quad ( \Re e \, p > \, 0 \, ),$$

$$\int _ { 0 } ^ { \infty } \frac { e ^ { - x ^ { 2 } / 4 t } } { \sqrt { \pi \, t } } \, e ^ { - ( s + \alpha ) t } \, d t = \frac { e ^ { - x \, \sqrt { s + \alpha } } } { \sqrt { s + \alpha } }$$

it results:

$$\hat { H } _ { 1 } ( x, s ) = \frac { 1 } { 2 \sqrt { \pi D } } \int _ { 0 } ^ { \infty } e ^ { - \frac { x ^ { 2 } } { 4 D y } - ( s + a + \frac { e } { s + \beta \varepsilon } ) y } \frac { d y } { \sqrt { y } } = \frac { 1 } { 2 \sqrt { D } } \frac { e ^ { - \frac { | x | } { \sqrt { D } } r } } { r }$$

where

$$r ^ { 2 } = s + a \, + \, \frac { \varepsilon } { s + \beta \varepsilon }.$$

Besides, since Fubini -Tonelli theorem and (2.19) one has:

$$\hat { H } _ { 2 } = \hat { H } _ { 1 } - \frac { 1 } { 2 \, \sqrt { \pi D } } \, \int _ { 0 } ^ { \infty } \, e ^ { - \frac { x ^ { 2 } } { 4 D y } - ( s + a + \frac { \varepsilon } { s + \beta \varepsilon } + \frac { \delta } { s + d \delta } ) \, y } \, \frac { d y } { \sqrt { y } }$$

from which, since (2.20),

$$\hat { H } \left ( x, s \right ) \, = \, \frac { 1 } { 2 \, \sqrt { D } } \ \frac { e ^ { - \frac { | x | } { \sqrt { D } } \sigma } } { \sigma }$$

is deduced.

Besides, by means of property of convolution for which f ∗ g = g ∗ f, since (2.14) and (2.15), property a) is evident. Moreover, properties b) and c) are proved following theorem 3.2.1 of [25]. In particular, as for property c), for | x | ≥ η and since | J 1 ( z ) | ≤ 1 , it results:

$$| H _ { 1 } | \, = \, \frac { e ^ { - \frac { \eta ^ { 2 } } { 4 D t } } } { 2 \sqrt { \pi D t } } \ e ^ { - a \, t } + \sqrt { \frac { \varepsilon \, t } { \pi D } } ;$$

$$| H _ { 2 } | \, \leq \sqrt { \frac { \delta \, t } { 4 \pi D } } + 2 \, \sqrt { \frac { \varepsilon \delta } { \pi D } } \, t$$

$$& ( 2. 2 6 ) & \quad | H _ { 2 } | \ \leq \sqrt { \frac { \partial \, t } { 4 \pi D } } + 2 \, \sqrt { \frac { \varepsilon \partial } { \pi D } } \ t \\ & \quad \text{from which property follows.}$$

Now, introducing the following functions:

$$\varphi, \text{introducing the following functions} \colon \\ \begin{cases} \quad \varphi ( x, t ) = \frac { e ^ { - \frac { x ^ { 2 } } { 4 D t } } } { 2 \sqrt { \pi D t } } \ e ^ { - a t } ; \\ \quad \psi _ { \varepsilon } ( y, t ) = \frac { \sqrt { \varepsilon } \, y } { \sqrt { t } - y } J _ { 1 } ( \, 2 \, \sqrt { \varepsilon } \, y \, ( t - y ) \, ; \\ \quad \psi _ { \delta } ( y, t ) = \frac { \sqrt { \delta } \, y } { \sqrt { t } - y } J _ { 1 } ( \, 2 \, \sqrt { \delta } \, y \, ( t - y ) \, ; \\ \quad \text{s} \colon \end{cases}$$

it results:

$$H _ { 1 } ( x, t ) & = \varphi ( x, t ) - \int _ { 0 } ^ { t } \varphi ( x, y ) \ \psi _ { \varepsilon } ( y, t ) \, d y$$

$$H _ { 2 } ( x, t ) = \int _ { 0 } ^ { t } H _ { 1 } ( x, y ) \ \psi _ { \delta } ( y, t ) \ d y.$$

Moreover, by denoting

$$g _ { 1 } ( x, t ) \ast g _ { 2 } ( x, t ) = \int _ { 0 } ^ { t } g _ { 1 } ( x, t - \tau ) g _ { 2 } ( x, \tau ) \, d \tau$$

the convolution with respect to t, for t &gt; 0 , as proved in [26] by means of formula (20),(21) and (24), it results:

$$( 2. 3 0 ) \quad \ \ ( \partial _ { t } + a - D \partial _ { x x } ) H _ { 1 } = - \varepsilon e ^ { - \varepsilon \beta t } * H _ { 1 } ( x, t ) = - \varepsilon K _ { \varepsilon } ( x, t )$$

where K ε is given by

$$K _ { \varepsilon } ( x, t ) \, = \, \frac { 1 } { 2 \, \sqrt { \pi \, D } } \, \int _ { 0 } ^ { t } \, e ^ { - \frac { x ^ { 2 } } { 4 D y } \, - a \, y \, - \beta \varepsilon ( t - y ) } \, \ J _ { 0 } \, ( 2 \, \sqrt { \varepsilon \, y \, ( t - y ) } \, ) \, \frac { d y } { \sqrt { y } }.$$

Hence, the following theorem can be proved:

Theorem 2. For t &gt; 0 , it results LH = 0 , i.e.

$$( 2. 3 2 ) \quad \ \ H _ { t } - D H _ { x x } + a H + \int _ { 0 } ^ { t } [ \varepsilon \, e ^ { - \varepsilon \beta ( t - \tau ) } \, + \delta e ^ { - \delta d ( t - \tau ) } ] H ( x, \tau ) \, d \tau = 0.$$

Proof. Let us consider that:

$$( \partial _ { t } + a - D \partial _ { x x } ) H _ { 2 } & = H _ { 1 } ( x, t ) \psi _ { \delta } ( t, t ) + \int _ { 0 } ^ { t } H _ { 1 } ( x, y ) \left [ \, \partial _ { t } \psi _ { \delta } ( y, t ) + a \psi _ { \delta } ( y, t ) \right ] \\ & \quad - D \, \int _ { 0 } ^ { t } \psi _ { \delta } ( y, t ) \, \partial _ { x x } H _ { 1 } ( x, y ) \ d y.$$

Accordingly, given relation (2.30), one has

$$( \partial _ { t } + a - D \partial _ { x x } ) H _ { 2 } & = \int _ { 0 } ^ { t } \left [ H _ { 1 } ( x, y ) \ \partial _ { t } \psi _ { \delta } ( y, t ) - \psi _ { \delta } ( y, t ) \ \partial _ { y } H _ { 1 } ( x, y ) \right ] d y + \\ ( 2. 3 4 ) & + H _ { 1 } ( x, t ) \psi _ { \delta } ( t, t ) - \varepsilon \int _ { 0 } ^ { t } K _ { \varepsilon } ( x, y ) \psi _ { \delta } ( y, t ) d y. \\ - \ \. \quad \. \quad.$$

Besides, considering that:

$$( 2. 3 5 ) \quad \int _ { 0 } ^ { t } \psi _ { \delta } ( y, t ) \partial _ { y } H _ { 1 } ( x, y ) \, d y = H _ { 1 } ( x, t ) \psi _ { \delta } ( t, t ) - \int _ { 0 } ^ { t } H _ { 1 } ( x, y ) \, \partial _ { y } \psi _ { \delta } ( y, t ) \, d y,$$

one has:

$$( \partial _ { t } + a - D \partial _ { x x } ) H _ { 2 } = \int _ { 0 } ^ { t } H _ { 1 } ( x, y ) \, \left [ \partial _ { t } \psi _ { \delta } ( y, t ) + \partial _ { y } \psi _ { \delta } ( y, t ) \right ] d y$$

$$- \varepsilon \int _ { \infty } ^ { t } K _ { \varepsilon } ( x, y ) \psi _ { \delta } ( y, t ) d y$$

$$( \partial _ { t } + a - D \partial _ { x x } ) H _ { 2 } = \int _ { 0 } ^ { t } H _ { 1 } ( x, \\ \\ ( 2. 3 6 ) & & - \varepsilon \int _ { 0 } ^ { t } K _ { \varepsilon } ( x,$$

where it results:

$$\partial _ { t } \psi _ { \delta } ( y, t ) + \partial _ { y } \psi _ { \delta } ( y, t ) = \delta e ^ { - \delta d \left ( t - y \right ) } J _ { 0 } ( \, 2 \, \sqrt { \delta \, y \left ( t - y \right ) } \,.$$

So that, denoting by

$$K _ { \delta } ( x, t ) \equiv \int _ { 0 } ^ { t } \, e ^ { - \delta d \, ( t - y ) } \, H _ { 1 } ( x, y ) \, J _ { 0 } \, ( \, 2 \, \sqrt { \delta \, y ( t - y ) } \, ) \, d y,$$

one has:

$$( \partial _ { t } + a - D \partial _ { x x } ) H _ { 2 } = \delta K _ { \delta } - \varepsilon \int _ { 0 } ^ { t } K _ { \varepsilon } ( x, y ) \psi _ { \delta } ( y, t ) d y.$$

Consequently, since for equation (2.30) one has K x,t ε ( ) = H x,t 1 ( ) ∗ e -εβ t , by means of Fubini -Tonelli theorem and (2.28), it is proved that:

$$( 2. 3 9 ) & & ( \partial _ { t } + a - D \partial _ { x x } ) H _ { 2 } = \delta K _ { \delta } - \varepsilon \, e ^ { - \varepsilon \beta t } * H _ { 2 } ( x, t ).$$

On the other hand, the convolution e -δd t ∗ H x,t ( ) is given by

$$e ^ { - \delta d t } * H ( x, t ) = e ^ { - \delta d t } * H _ { 1 } ( x, t ) - \int _ { 0 } ^ { t } H _ { 1 } ( r, y ) d y \int _ { y } ^ { t } e ^ { - \delta d ( t - \tau ) } \psi _ { \delta } ( y, \tau ) \, d \tau$$

with

$$( 2. 4 0 ) \quad \int _ { y } ^ { t } e ^ { - \delta d ( t - \tau ) } \psi _ { \delta } ( y, \tau ) d \tau = e ^ { - \delta d ( t - y ) } \int _ { y } ^ { t } \sqrt { \frac { \delta \, y } { \tau - y } } \ J _ { 1 } \left ( 2 \, \sqrt { \delta \, y ( \tau - y ) } \, \right ) \, d \tau =$$

$$= \ e ^ { - \, \delta d \, ( t - y ) } \, \left \lceil \, 1 \, - \, J _ { 0 } \, ( \, 2 \, \sqrt { \delta \, y ( t - y ) } \, ) \, \right \rceil.$$

As a consequence, it results:

$$e ^ { - \delta d t } * H = K _ { \delta }.$$

Therefore, given relations (2.30), (2.39), (2.41), theorem holds.

## 3 Solution related to the (FHR) problem

To provide the solution by means of the integral expression (2.13), some convolutions need to be determined.

In order to evaluate ∫ t 0 dτ ∫ glyph[Rfractur] H ξ, τ ( ) dξ, let us start to observe that

$$\int _ { \Re } d \xi \int _ { 0 } ^ { t } H _ { 2 } ( \xi, \tau ) d \tau = \int _ { \Re } d \xi \int _ { 0 } ^ { t } H _ { 1 } ( \xi, y ) d y \int _ { y } ^ { t } \psi _ { \delta } ( y, \tau ) d \tau$$

with

(3.42)

$$- \delta \, d \int _ { y } ^ { t } e ^ { - \delta \, d ( t - y ) } J _ { 0 } ( 2 \sqrt { \delta y ( t - y ) } d \tau.$$

Consequently, for (2.38), one has:

$$\int _ { \Re } d \xi \int _ { 0 } ^ { t } H _ { 2 } ( \xi, \tau ) d \tau & = \int _ { \Re } d \xi \int _ { 0 } ^ { t } H _ { 1 } ( \xi, \tau ) d \tau - \int _ { \Re } K _ { \delta } ( \xi, t ) d \xi + \\ - \delta d \int _ { \Re } d \xi \int _ { 0 } ^ { t } K _ { \delta } ( \xi, \tau ) d \tau \\ \varepsilon \quad. \quad. \quad. \quad. \quad. \quad \varepsilon \quad. \quad. \quad.$$

So that, according to (2.17), it results:

$$( 3. 4 4 ) \quad \int _ { 0 } ^ { t } \, d \tau \, \int _ { \Re } \, H ( \xi, \tau ) \ d \xi = \int _ { \Re } K _ { \delta } ( \xi, t ) d \xi \, + \delta d \, \int _ { \Re } d \xi \int _ { 0 } ^ { t } K _ { \delta } ( \xi, \tau ) d \tau.$$

Now, let us evaluate

$$\int _ { y } ^ { t } \psi _ { \delta } ( y, \tau ) d \tau = 1 - e ^ { - \delta d ( t - y ) } J _ { 0 } ( 2 \sqrt { \delta y ( t - y ) } +$$

$$e ^ { - \varepsilon \beta t } * H _ { 2 } = \int _ { 0 } ^ { t } \, H _ { 1 } \left ( x, y \right ) d y \, \int _ { y } ^ { t } \, e ^ { - \beta \varepsilon \left ( t - \tau \right ) } \, \psi _ { \delta } \left ( y, \tau \right ) d \tau.$$

Considering (2 27) , after an integration by parts, one obtains: . 3

$$e ^ { - \varepsilon \beta t } * H _ { 2 } & = e ^ { - \beta \varepsilon t } * H _ { 1 } ( x, t ) - \int _ { 0 } ^ { t } e ^ { - \delta d ( t - y ) } H _ { 1 } ( x, y ) J _ { 0 } ( 2 \sqrt { \delta y ( t - y ) } ) d y + \\ ( 3. 4 5 ) & = \int _ { \real } ^ { \real } \Big { w }$$

$$\theta _ { 0 } ^ { \prime } \left ( \varepsilon \beta - \delta d \right ) \int _ { 0 } ^ { t } d \tau \int _ { 0 } ^ { \tau } H _ { 1 } \left ( x, y \right ) e ^ { - \delta d \left ( \tau - y \right ) } \, e ^ { - \varepsilon \beta \left ( t - \tau \right ) } \, J _ { 0 } \left ( \, 2 \, \sqrt { \delta \, y ( \tau - y ) } \, \right ) d y,$$

and, for (2.38), it results:

$$( 3. 4 6 ) \quad \ e ^ { - \varepsilon \beta t } * H _ { 2 } = e ^ { - \varepsilon \beta t } * H _ { 1 } - K _ { \delta } + ( \varepsilon \beta - \delta d ) e ^ { - \beta \varepsilon t } * K _ { \delta }.$$

Moreover, since (2.17) and (3.46), one deduces:

$$e ^ { - \varepsilon \beta t } * H = K _ { \delta } + ( \delta d - \varepsilon \beta ) e ^ { - \beta \varepsilon t } * K _ { \delta }.$$

Now, let us denote by

$$f _ { 1 } ( x, t ) \ast f _ { 2 } ( x, t ) = \int _ { \Re } f _ { 1 } ( \xi, t ) f _ { 2 } ( x - \xi, t ) \, d \xi$$

the convolution with respect to the space x, and let

$$H \otimes F \, = \, \int _ { 0 } ^ { t } \, d \tau \, \int _ { \Re } \, H ( x - \xi, t - \tau ) \ F \left [ \xi, \tau, u ( \xi, \tau ) \right ] d \xi.$$

Considering (2.41) and (3.47) one has:

$$\begin{cases} \ H \otimes e ^ { - \delta d t } = \int _ { \Re } K _ { \delta } ( \xi, t ) \ d \xi, \\ \ H \otimes e ^ { - \beta \varepsilon t } = \int _ { \Re } \left [ K _ { \delta } + ( \delta d - \varepsilon \beta ) e ^ { - \beta \varepsilon t } * K _ { \delta } \right ] d \xi. \\ \quad \ \. \end{cases}$$

Moreover, it results:

$$\begin{cases} \ H \otimes ( y _ { 0 } ( x ) \, e ^ { - \delta d t } ) = y _ { 0 } \ast K _ { \delta } \\ \ H \otimes ( w _ { 0 } ( x ) \, e ^ { - \beta \varepsilon t } ) = w _ { 0 } \ast [ K _ { \delta } + ( \delta d - \varepsilon \beta ) e ^ { - \beta \varepsilon t } \ast K _ { \delta } \, ] \end{cases}$$

where

$$w _ { 0 } \ast ( \delta d - \varepsilon \beta ) e ^ { - \beta \varepsilon \, t } \ast K _ { \delta } = ( \delta d - \varepsilon \beta ) \, w _ { 0 } \, \otimes ( e ^ { - \beta \varepsilon \, t } K _ { \delta } \, ).$$

Consequently, given (1.9) and (2.13), for (3.44),(3.50),(3.51) and (3.52), it results:

$$u ( x, t ) \, = \, u _ { 0 } ( x ) \ast H \, + ( y _ { 0 } ( x ) - \, w _ { 0 } ( x ) ) \ast K _ { \delta } \, + \, \varphi ( u ) \otimes H \, + \\ + ( \varepsilon \beta - \delta d ) w _ { 0 } ( x ) \otimes e ^ { - \beta \varepsilon t } K _ { \delta } + \frac { c } { \beta } ( \delta d - \varepsilon \beta ) e ^ { - \beta \varepsilon t } \otimes K _ { \delta } + \\ + \left ( \frac { h } { d } - \frac { c } { \beta } \right ) \delta d \otimes K _ { \delta } \\ \right ) \, + \, \left ( \frac { h } { d } - \frac { c } { \beta } \right ) \delta d \otimes K _ { \delta } \right ) \right ) \, + \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cd: \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdclots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \collections$$

and this formula, together with relations (1.8), allow us to determine also v x, t ( ) and y x, t ( ) in terms of the data.

## 4 Explicit solutions

Several methods have been developed to find exact solutions related to partial differential equations [27-31]. In this case, by referring to [7], an extra term is added in order to achieve some solutions. Accordingly, let us consider

$$( 4. 5 4 ) & \begin{cases} \begin{array} { c } \partial u \\ \overline { \partial t } \end{array} = \begin{array} { c } D \, \frac { \partial ^ { 2 } u } { \partial x ^ { 2 } } \, - \, w \, + y \, + f ( u ) \\ \\ \overline { \partial w } \end{array} = \begin{array} { c } \varepsilon ( - \beta w + c + u ) + k \, u ^ { 2 } \\ \\ \overline { \partial y } \end{array} = \begin{array} { c } \delta ( - u + h - d y ) \end{array} \\ \text{where $k \neq 0$ and let us assume $\varepsilon \beta = \delta d$ and $f ( u ) = 2 u ( a - u ) ( u - 1 ) } \end{cases}$$

glyph[negationslash]

where k = 0 and let us assume εβ = δ d and f ( u ) = 2 u a ( -u ) ( u -1) . Under these conditions, problem (1.10) turns into:

$$( 4. 5 5 ) \quad \begin{cases} \ u _ { t } - D u _ { x x } + 2 a u + ( \varepsilon + \delta ) \int _ { 0 } ^ { t } e ^ { - \varepsilon \beta ( t - \tau ) } \, u ( x, \tau ) \, d \tau \, = \, F ( x, t, u ) \\ \ u ( x, 0 ) \, = u _ { 0 } ( x ) \quad x \, \in \, \Re, \end{cases}$$

where, by denoting ϕ 1 = 2 u 2 ( a +1 -u ) + k ∫ t 0 e -ε β ( t -τ ) u 2 ( x, τ ) dτ, it results:

$$F = \varphi _ { 1 } ( u ) - w _ { 0 } e ^ { - \varepsilon \beta t } + y _ { 0 } e ^ { - \delta d t } - \frac { c } { \beta } ( 1 - e ^ { - \varepsilon \beta t } ) + \frac { h } { d } ( 1 - e ^ { - \delta d t } )$$

and

$$( 4. 5 7 ) \quad \begin{cases} \ w = w _ { 0 } e ^ { - \varepsilon \beta t } + \frac { c } { \beta } ( 1 - e ^ { - \varepsilon \beta t } ) + \int _ { 0 } ^ { t } e ^ { - \varepsilon \beta ( t - \tau ) } \left [ \varepsilon \, u ( x, \tau ) + k u ^ { 2 } ( x, t ) \right ] d \tau \\ \ v = y _ { 0 } \, e ^ { - \delta \, d \, t } \, + \, \frac { h } { d } \, ( 1 - e ^ { - \delta \, d \, t } ) \, - \delta \int _ { 0 } ^ { t } e ^ { - \delta \, d \, ( t - \tau ) } \, u ( x, \tau ) \, d \tau. \\ \psi \quad, \cdots \, \hat { \ }, \cdots \, \hat { \ }, \cdots \, \hat { \ }, \cdots \, \hat { \ }, \cdots \, \end{cases}$$

In order to find explicit solutions, let us introduce

$$z = x - C \, t,$$

obtaining, from system (4.54), the following equation:

$$& D C \, u _ { z z z } + ( C ^ { 2 } - \varepsilon \beta \, D ) u _ { z z } - 6 C u ^ { 2 } u _ { z } + 4 C ( a + 1 ) u \, u _ { z } + 2 \varepsilon \beta u ^ { 3 } + k \, u ^ { 2 } \\ & - C \, ( 2 \, a \, + \varepsilon \beta ) \, u _ { z } - 2 \, \varepsilon \beta ( a + 1 ) u ^ { 2 } + 2 \, \varepsilon \beta \, a \, u + ( \varepsilon + \delta ) \, u + \varepsilon \, c - \delta h = 0.$$

Now, let us consider

$$( 4. 5 9 ) & & f ( z ) = \sqrt { y } \, \tanh { ( \sqrt { y } \, ( z - z _ { 0 } ) ) },$$

that is a solution of Riccati type equation:

$$f _ { z } + f ^ { 2 } - y = 0,$$

and let us assume

$$u ( z ) & = A \ f ( z ) + b.$$

Since

$$u _ { z } = - A \, f ^ { 2 } ( z ) + A \, y \\ u _ { z z } = 2 \, A f ^ { 3 } - 2 \, A \, f \, y \\ u _ { z z z } = - 6 A f ^ { 4 } ( z ) + 8 \, A \, y f ^ { 2 } ( z ) - 2 \, A \, y ^ { 2 } \\ ( 4. 6 1 ) \quad \ u \, u _ { z } = - A ^ { 2 } \, f ^ { 3 } ( z ) - A b f ^ { 2 } + A ^ { 2 } y \, f + A b y \\ u ^ { 2 } \, u _ { z } = - A ^ { 3 } f ^ { 4 } - 2 A ^ { 2 } \, b \, f ^ { 3 } + ( A ^ { 3 } y - A b ^ { 2 } ) f ^ { 2 } + 2 A ^ { 2 } b \, y f + A b ^ { 2 } y, \\ \cdot \quad \cdot \quad \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \end{$$

in order to satisfy equation (4.58), constant b needs to assume the following expression:

$$b & = \frac { \varepsilon \beta \, A } { 2 C } - \frac { \varepsilon \, + \, \delta } { 2 \, k }$$

and, moreover, it has to be

$$D y$$

$$A ^ { 2 } = D \\ a + 1 = 3 b + \frac { C \, A } { 2 \, D } \\ D \, y = 3 \, b ^ { 2 } + \left ( \frac { C } { A } - 3 - \frac { k } { \varepsilon \beta } \right ) b - \frac { C } { 2 A } - \frac { \varepsilon + \delta } { 2 \, \varepsilon \, \beta } + 1$$

and

$$\varepsilon \, c - \delta \, h = 2 \, C D A y ^ { 2 } + \varepsilon \beta \, C A y - 2 \varepsilon \beta \, b ^ { 3 } + 2 \, \varepsilon \beta \, ( a + 1 ) \, b ^ { 2 } + \\ - 2 \, \varepsilon \beta \, a \, b + 6 \, C \, A \, y \, b ^ { 2 } - 4 C \, y \, ( a + 1 ) A b + 2 \, a \, C \, A \, y \, - \, ( \varepsilon + \delta ) \, b - k \, b ^ { 2 }.$$

So, if for instance, it is assumed εβ = 1 , z 0 = 0, ε + δ = 0 04 . , k = 0 01 and . C = 1 with D &gt; 0 019 . , for A = √ D, by introducing

$$g ( D ) = \sqrt { 3 9 0 0 \sqrt { D } - 1 5 0 1 \, D + 1 5 0 \, D \, \sqrt { D } - 5 0 0 },$$

equation (4.60) gives

$$u ( z ) = \frac { \sqrt { 2 } } { 2 0 \, D ^ { 1 / 4 } } \ g ( D ) \, \tanh \left ( \frac { \sqrt { 2 } \ z \, g ( D ) } { 2 0 \, D ^ { 3 / 4 } } \right ) + \frac { \sqrt { D } } { 2 } - 2.$$

In Fig.1, solutions u(z) expressed by means of formula (4.64) are illustrated for different values of D, by showing that the amplitude increases as 0 &lt; D &lt; 1 increases.

Remark When fast variable u simulates the membrane potential of a nerve cell, while slow variable w and superslow variable y determine the corresponding ion

Figure 1: Solution u(z) when εβ = 1 , z 0 = 0, ε + δ = 0 04 . , k = 0 01 . , C = 1 , A = √ D and for D = 0 02 . , D = 0 05 . , D = 0 5 . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00997v4-figure-20e3893b-6325-47d0-8cae-4fedffc78b10.jpg)

number densities, model (4.54) with its solutions can be of interest in applications to understand how impulses are propagated from one neuron to another. Moreover, as similarly underlined in [4], the knowledge of exact solutions can help in testing different applications of models in neuroscience.

## Acknowledgements

The present work has been developed with the economic support of MIUR (Italian Ministry of University and Research) performing the activities of the project ARS01 -00861 'Integrated collaborative systems for smart factory - ICOSAF'. The paper has been performed under the auspices of G.N.F.M. of INdAM.

The authors are grateful to anonymous referees for their helpful comments.

## References

- [1] R. Bertram, T. Manish J. Butte,T. Kiemel and A. Sherman, Topological and phenomenological classification of bursting oscillations, Bull. Math.Biol, Vol. 57, No. 3, pp. 413, (1995)
- [2] J.Wojcik, A. Shilnikov, Voltage Interval Mappings for an Elliptic Bursting Model in Nonlinear Dynamics New Directions Theoretical Aspects Gonz´leza Aguilar H; Ugalde E. (Eds.) 219 pp Springer,(2015)
- [3] Zemlyanukhin A. I., Bochkarev A. V., Analytical Properties and Solutions of the FitzHugh - Rinzel Model, Rus. J. Nonlin. Dyn., Vol. 15, no. 1, pp. 3-12, (2019)

| [4]   | N. A. Kudryashov, On Integrability of the FitzHugh - Rinzel Model,Russian Journal of Nonlinear Dynamics, 2019, vol. 15, no. 1, pp. 13-19.                                                                                                               |
|-------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [5]   | Izhikevich E.M., Dynamical Systems in Neuroscience: The Geometry of Ex- citability and Bursting,p.397. The MIT press. England (2007)                                                                                                                    |
| [6]   | M.De Angelis, P. Renno, Asymptotic effects of boundary perturbations in excitable systems, Discrete and continuous dynamical systems series B, 19, no 7 (2014)                                                                                          |
| [7]   | N.K.Kudryashov, K.R. Rybka, A. Sboev, Analytical properties of the perturbed FitzHugh-Nagumo model, Applied Mathematics Letters, 76 , 142-147 ,(2018)                                                                                                   |
| [8]   | M. De Angelis, On a model of superconductivity and biology, Advances and Applications in Mathematical Sciences Volume 7, Issue 1, Pages 41-50, (2010)                                                                                                   |
| [9]   | Rionero, S. A rigorous reduction of the L 2 -stability of the solutions to a nonlin- ear binary reaction-diffusion system of PDE's to the stability of the solutions to a linear binary system of ODE's, J. Math. Anal. Appl. 319 no. 2, 377-397 (2006) |
| [10]  | Rionero, S, Torcicollo,I. On the dynamics of a nonlinear reaction-diffusion duopoly mode,International Journal of Non-Linear Mechanics,Volume 99, Pages 105-111,(2018)                                                                                  |
| [11]  | M. De Angelis, Asymptotic Estimates Related to an Integro Differential Equa- tion, Nonlinear Dynamics and Systems Theory, 13 (3) 217-228 (2013)                                                                                                         |
| [12]  | G. Gambino,M. C. Lombardo,G. Rubino, M. Sammartino, Pattern selection in the 2D FitzHugh-Nagumo model, Ric. di Mat. https://doi.org/10.1007/s11587 1-15 (2018)                                                                                          |
| [13]  | Keener, J. P. Sneyd,J. Mathematical Physiology . Springer-Verlag, N.Y, 470 pp, (1998)                                                                                                                                                                   |
| [14]  | A. Quarteroni, A. Manzoni and C. Vergara, The cardiovascular sys- tem:Mathematical modelling, numerical algorithms and clinical applications, Acta Numerica, vol 26,pp. 365-590,(2017)                                                                  |
| [15]  | F. De Angelis, D. Cancellara, L. Grassia, A. D'Amore, The influence of load- ing rates on hardening effects in elasto/viscoplastic strain-hardening materials, Mechanics of Time-Dependent Materials, 22 (4) 533-551, (2018)                            |
| [16]  | Carillo, S.; Chipot, M.; Valente, V.; Vergara Caffarelli, G. On weak regularity requirements of the relaxation modulus in viscoelasticity. Commun. Appl. Ind. Math. 10 , no. 1, 78-87,(2019)                                                            |

| [17]   | F. De Angelis, D. Cancellara, Dynamic analysis and vulnerability reduction of asymmetric structures: Fixed base vs base isolated system, Composite Struc- tures, 219 203-220, (2019).                                                                          |
|--------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [18]   | F. De Angelis, Extended formulations of evolutive laws and constitutive rela- tions in non-smooth plasticity and viscoplasticity, Composite Structures, 193 35-41, (2018).                                                                                     |
| [19]   | H. Simo, P. Woafo,Bursting oscillations in electromechanical systems, Mechan- ics Research Communications 38 537 - 541, (2011)                                                                                                                                 |
| [20]   | E.Juzekaeva, A. Nasretdinov, S. Battistoni, T. Berzina, S. Iannotta, R. Khazipov, V. Erokhin, M. Mukhtarov, Coupling Cortical Neurons through Elec- tronic Memristive Synapse, Adv. Mater. Technol. 4, 1800350 (6) (2019)                                      |
| [21]   | F. Corinto, V. Lanza, A. Ascoli, and Marco Gilli, Synchronization in Net- works of FitzHugh-Nagumo Neurons with Memristor Synapses, in 20th Euro- pean Conference on Circuit Theory and Design (ECCTD) IEEE. (2011)                                            |
| [22]   | M De Angelis, A priori estimates for excitable models,Meccanica, Volume 48, Issue 10, pp 2491-2496 (2013)                                                                                                                                                      |
| [23]   | Murray, J.D. Mathematical Biology I, . Springer-Verlag, N.Y, 767 pp, (2003)                                                                                                                                                                                    |
| [24]   | Wenxian Xie, Jianwen Xu, Li Cai, Yanfei Jin, Dynamics and Geometric Desin- gularization of the Multiple Time Scale FitzHugh Nagumo Rinzel Model with Fold Singularity , Communications in Nonlinear Science and Numerical Simu- lation, 63, p. 322-338. (2018) |
| [25]   | J. R. Cannon, The one-dimensional heat equation, Addison-Wesley Publishing Company 483 pp, (1984)                                                                                                                                                              |
| [26]   | De Angelis, M. Renno, P, Existence, uniqueness and a priori estimates for a non linear integro - differential equation, Ricerche di Mat. 57 95-109 (2008)                                                                                                      |
| [27]   | M. De Angelis, On the transition from parabolicity to hyperbolicity for a nonlinear equation under Neumann boundary conditions, Meccanica, Volume 53, Issue 15, pp 3651-3659, (2018)                                                                           |
| [28]   | Li H., Guoa Y.: New exact solutions to the Fitzhugh Nagumo equation, Applied Mathematics and Computation 180 , 2, 524-528 (2006)                                                                                                                               |
| [29]   | M. De Angelis, A wave equation perturbed by viscous terms: fast and slow times diffusion effects in a Neumann problem,Ricerche di Matematica, Volume 68, Issue 1, pp 237-252, (2019)                                                                           |

- [30] B.Prinaria, F.Demontis, Sitai Li, T.P.Horikis, Inverse scattering transform and soliton solutions for square matrix nonlinear Schr¨dinger equations with nono zero boundary conditions,Physica D: Nonlinear Phenomena,Volume 368, Pages 22-49, (2018)
- [31] N.K.Kudryashov, Asymptotic and Exact Solutions of the FitzHugh-Nagumo Model, Regul. Chaotic Dyn., vol 23, No 2, 152-160, (2018)
Fabio De Angelis Monica De Angelis
On solutions related to FitzHugh-Rinzel type model
A ternary autonomous dynamical system of FitzHugh-Rinzel type is analyzed.
The system, at start, is reduced to a nonlinear integro differential equation.
The fundamental solution $ H(x,t) $ is explicitly determined and the initial
value problem is analyzed in the whole space. The solution is expressed by
means of an integral equation involving $ H(x,t) $. Moreover, adding an extra
control term, explicit solutions are achieved.