A Novel Numerical Approach for a Nonlinear Fractional Dynamical Model of Interpersonal and Romantic Relationships

Abstract

In this paper, we propose a new numerical algorithm, namely q-homotopy analysis Sumudu transform method (q-HASTM), to obtain the approximate solution for the nonlinear fractional dynamical model of interpersonal and romantic relationships. The suggested algorithm examines the dynamics of love affairs between couples. The q-HASTM is a creative combination of Sumudu transform technique, q-homotopy analysis method and homotopy polynomials that makes the calculation very easy. To compare the results obtained by using q-HASTM, we solve the same nonlinear problem by Adomian’s decomposition method (ADM). The convergence of the q-HASTM series solution for the model is adapted and controlled by auxiliary parameter $\hslash$ and asymptotic parameter n. The numerical results are demonstrated graphically and in tabular form. The result obtained by employing the proposed scheme reveals that the approach is very accurate, effective, flexible, simple to apply and computationally very nice.

1. Introduction

The theory of entropy was originally related with thermodynamics, but in recent years it has also been used in other areas of investigation such as information theory, psychodynamics, thermoeconomics, human relationships and many more. The second law of thermodynamics states that entropy increases with time. It indicates the instability of a system over a period of time if there is nothing to stabilize it. Similarly, in human relationships, we have daily interactions and these relationships also become disordered. We require input to stabilize any relationship, to iron out the wrinkles or differences, so that we don’t possess and store things forever [1,2,3,4]. The study of interpersonal and romantic relations has been gaining popularity during the past ten years. Interpersonal relations occur in many ways, for example marriage, blood relations, close associations, work and clubs. Since 1957, marriage has been investigated scientifically and we can observe some general interpretations which lead to models of marital interactions. The authors were inspired to study why some married couples get a divorce, when some others couples do not divorce. Moreover, among married couples, some are happy, while some are not happy with each other. In today’s scenario around the globe the number of divorce cases in interpersonal and romantic relations for marriage is increasing day to day. A survey in the United States revealed that for a 40 year duration the chances a first marriage will end in divorce are approximately 50 to 67 percent. For the case of a second marriage for the same duration the chance it will end in divorce is 10 percent higher than for a first marriage. Worldwide, the United States has the highest number of divorce cases. In these areas, experiments are inconvenient and may be restricted due to ethical considerations. To study the dynamics of interpersonal and romantic relations in marriage mathematical models can play a key role and in recent years, many scientists and researchers have studied dynamical modelling of interpersonal relations [5,6,7,8].

Fractional order models are very important to study natural problems. It is well known that the nature of the trajectory of the fractional order derivatives is non-local, which describes that the fractional order derivative has memory effect features [9,10,11,12,13], and any dynamical or physical system related to fractional order differential operators has a memory effect, meaning that the future states depend on the present as well as the past states [14,15,16,17,18,19]. In view of the great importance of fractional approaches in real world problems we were motivated to study the nonlinear dynamical model of interpersonal and romantic relationships. There are various analytical and numerical techniques to study such problems. In 1992, Liao proposed an analytic technique known as the homotopy analysis method (HAM) [20,21] for managing nonlinear problems. More recently, the authors of [22,23] have suggested an extension of HAM known as q-homotopy analysis method (q-HAM) to discuss nonlinear mathematical models. Standard analytical schemes requisite more computer memory and computational time, thus to overcome these limitations the analytical methods need to be merged with standard integral transform operators to study the nonlinear equations appearing in science and engineering [24,25,26].The key purpose of the present work is to propose a novel numerical approach namely q-homotopy analysis Sumudu transform method (q-HASTM) to solve the nonlinear fractional dynamical model of interpersonal and romantic relations in marriages. The q-HASTM is an elegant amalgamation of the q-HAM, the standard Sumudu transform and the homotopy polynomials. The supremacy of this algorithm is its potential of merging two robust computational approaches for investigating nonlinear fractional differential equations. Moreover, the proposed technique contains an asymptotic parameter n by which we can insure to convergence of a series of problem solutions. The paper is organized as follows: In Section 2, fractional calculus and the Sumudu transform are discussed. In Section 3, a fractional order nonlinear dynamical model of interpersonal and romantic relations in marriages is discussed. In Section 4, the fundamental plan of q-HASTM is proposed. In Section 5, a solution of the fractional model is obtained by employing q-HASTM. In Section 6, the fundamental plan of ADM is discussed. In Section 7, a solution of the mathematical model is obtained by employing ADM. Results and discussion are presented in Section 8. Lastly, in Section 9 the concluding observations are highlighted.

2. Basic Definitions

In this section, we discuss the fundamental definitions of integrals and derivatives of non-integer order and the Sumudu transform.

Definition 1.

The left-sided Riemann-Liouville fractional integral operator of order $\mu >0,$ of a function $g(\xi )\in C_{\alpha },\alpha \ge -1$ is given in the following manner [9]:

$$J^{\mu }g(\xi )=\frac{1}{\Gamma (\mu )}{\displaystyle {\int }_0^{\xi }{(\xi -\eta )}^{\mu -1}}g(\eta )d\eta ,(\mu >0),$$

(1)

Definition 2.

If $g(\xi )$ is a function of a variable $\xi$, then the Caputo fractional derivative of the function is expressed as follows [10]:

$$D_{\xi }^{\mu }g(\xi )=J^{m-\mu }{D}^{m}g(\xi )=\frac{1}{\Gamma (m-\mu )}{\displaystyle {\int }_0^{\xi }{(\xi -\eta )}^{m-\mu -1}}{g}^{(m)}(\eta )d\eta ,$$

(2)

for $m-1<\mu \le m,m\in N,\xi >0.$

Now, we present the following relation for the fractional integral in the Riemann-Liouville sense and the fractional derivative operator in the Caputo sense:

$$J_{\xi }^{\mu }{D}_{\xi }^{\mu }g(\xi )=g(\xi )-{\displaystyle \sum _{l=0}^{m-1}{g}^{(l)}}(0+)\frac{{\xi }^l}{l!}.$$

(3)

Definition 3.

The Sumudu transform is a newly introduced integral operator, which was suggested and developed by Watugala [27].This newly integral transform is expressed and discussed over the set of functions:

$$\Theta =\{g(\xi ):\exists \mathrm{N},{\delta }_1,{\delta }_2>0,\left|g(\xi )\right|<{\mathrm{N}\mathrm{e}}^{\left|\xi \right|/{\delta }_{\mathrm{j}}},\mathrm{if}\xi \in {(-1)}^{\mathrm{j}}\times [0,\infty )\}$$

as follows:

$$G(\mathrm{u})=S[g(\xi )]={\displaystyle {\int }_0^{\infty }g(\mathrm{u}\xi {)\mathrm{e}}^{-\xi }\mathrm{d}\xi },\mathrm{u}\in (-{\delta }_1,{\delta }_2).$$

(4)

The pioneer work in connection with the formation of important and basic results of the newly integral transform was conducted by Asiru [28], Belgacem et al. [29,30] and Srivastava et al. [31].

Definition 4.

If $g(\xi )$ is a function of $\xi$, then the Sumudu transform of the Caputo fractional operator $D_{\xi }^{\mu }g(\xi )$ is written as follows [32]:

$$S[D_{\xi }^{\mu }g(\xi )]=u^{-\mu }S[g(\xi )]-{\displaystyle \sum _{l=0}^{m-1}{u}^{-\mu +l}}{g}^{(l)}(0+),(m-1<\mu \le m).$$

(5)

3. Model Description

We consider a nonlinear dynamical model of interpersonal relations with arbitrary order of the form [8]:

$$\begin{gathered} D_{\xi }^{\mu }x(\xi )=-a_{1}x+b_{1}y(1-\epsilon y^2)+c_1,0<\mu \le 1, \\ {D}_{\xi }^{\mu }y(\xi )=-a_{2}y+b_{2}x(1-\epsilon x^2)+c_2,0<\mu \le 1, \end{gathered}$$

(6)

with initial conditions:

$$x(0)=0,y(0)=0.$$

Here variable $x$ measures the love of individual 1 for his/her partner and $y$ measures the love of individual 2 for his/her partner. The positive and negative measurement indicates the feelings of the partners towards each other. In the model, $a_i>0,a_i,b_i,c_i(i=1,2)$ and $\epsilon$ are real constants. These parameters are described as the oblivion, reaction, and attraction constants. In this model, we suppose that in the absence of partners feelings of each other decay exponentially fast. The romantic style of individuals 1 and 2 specified by the parameters $a_i,b_i$ and $c_i(i=1,2).$ For example, $a_i$ represents the expanse to which the individual $i$ is inspired by one's own feelings. The parameter $b_i$ describes the expanse to which the individual $i$ is inspired by the partner and it is expected that individual’s partner are supportive in nature. It measures the inclination to seek or avoid the closeness in the love affair between a couple. Hence, $-a_i{x}_i$ indicates that love measures of $i$ when the partner is not present, decay exponentially. The term 1/a_i indicates the time required for love to decay.

4. Basic Planof q-HASTM

To discuss the key plan of q-HASTM, we take a nonlinear fractional differential equation of the form:

$$D_{\xi }^{\mu }x(\xi )+Rx(\xi )+Nx(\xi )=\psi (\xi ),m-1<\mu \le m,$$

(7)

Here $x$ is an unknown function of variable $\xi$ and $D_{\xi }^{\mu }=\frac{d^{\mu }}{d{\xi }^{\mu }}$ denotes the fractional operator of order $\mu$ proposed by Caputo, $m\in N,R$ indicates the bounded linear operator in $\xi$ and $N$ represents the general nonlinear differential operator in $\xi$ Employing the Sumudu transform on Equation (7), we obtain:

$$S[D_{\xi }^{\mu }x]+S[Rx+Nx]=S[\psi (\xi )].$$

(8)

Applying the differentiation aspect of the Sumudu transform, we have the following result:

$$u^{-\mu }S[x]-{\displaystyle \sum _{l=0}^{m-1}{u}^{-\mu +l}{x}^{(l)}(0)}+S[Rx+Nx]=S[\psi (\xi )].$$

(9)

On simplification, we have:

$$S[x]-u^{\mu }{\displaystyle \sum _{l=0}^{m-1}{u}^{-\mu +l}{x}^{(l)}(0)}+u^{\mu }\left[S[Rx+Nx]-S[\psi (\xi )]\right]=0.$$

(10)

We define the nonlinear operator by keeping view of Equation (10):

$$\begin{gathered} \wp [\vartheta (\xi ;q)]=S[\vartheta (\xi ;q)]-u^{\mu }{\displaystyle \sum _{l=0}^{m-1}{u}^{-\mu +l}{\vartheta }^{(l)}(\xi ;q)(0)} \\ +u^{\mu }\left[S[R\vartheta (\xi ;q)+N\vartheta (\xi ;q)]-S[\psi (\xi )]\right], \end{gathered}$$

(11)

In Equation (11), $q\in \left[1,\frac{1}{n}\right]$ is the embedding parameter and $\vartheta (\xi ;q)$ is known as a real function of $\xi$ and $q.$ Now we develop the homotopy which is given by Equation (12):

$$(1-nq)S[\vartheta (\xi ;q)-x_0]=\hslash N[x(\xi )],$$

(12)

where S indicates the Sumudu transform operator, $n\ge 1,\hslash \ne 0$ specify an auxiliary parameter, $x_0(\xi )$ show an initial approximation of $x(\xi ).$ If we set the embedding parameter $q=0$ and $q=\frac{1}{n},$ then we have:

$$\vartheta (\xi ;0)=x_0(\xi ),\vartheta (\xi ;\frac{1}{n})=x(\xi ),$$

(13)

Thus, when $q$ varies from 0 to $\frac{1}{n}$, then the solution $\vartheta (\xi ;q)$ takes place from the initial approximation $x_0(\xi )$ to the solution $x(\xi ).$ Expanding the function $\vartheta (\xi ;q)$ in the series form by exerting the Taylor series about $q,$ we obtain:

$$\vartheta (\xi ;q)=x_0(\xi )+{\displaystyle \sum _{k=1}^{\infty }{x}_k(\xi )q^k},$$

(14)

where:

$$x_k(\xi )=\frac{1}{k!}\frac{{\partial }^k}{\partial q^k}\left\{\vartheta (\xi ;q)\right\}|{}_{q=0}.$$

(15)

If the initial approximation $x_0(\xi ),$ asymptotic parameter $n$ and the auxiliary parameter $\hslash$ are chosen properly, then Equation (14) converges at $q=\frac{1}{n}$, then we get the following result:

$$x(\xi )=x_0(\xi )+{\displaystyle \sum _{k=1}^{\infty }{x}_k(\xi ){\left(\frac{1}{n}\right)}^k},$$

(16)

which must be one of the solutions of the nonlinear fractional differential Equation (7). By taking into consideration Equation (16), the governing equation can be found from the Equation (12).

We define the vectors as follows:

$${\overrightarrow{x}}_k=\{x_0(\xi ),x_1(\xi ),x_2(\xi ),x_3(\xi ),\dots ,x_k(\xi )\}.$$

(17)

Now differentiating the equation (12) k-times w.r.t. $q$, then divide by k! and finally setting $q=0,$ it yields the following result:

$$S[x_k(\xi )-{\chi }_k{x}_{k-1}(\xi )]=\hslash {\Im }_k({\overrightarrow{x}}_{k-1}).$$

(18)

Now employing the inverse of the Sumudu transform operator on Equation (18), we obtain:

$$x_k(\xi )={\chi }_k{x}_{k-1}(\xi )+\hslash S^{-1}[{\Im }_k({\overrightarrow{x}}_{k-1})],$$

(19)

where ${\chi }_k$ is defined and given as follows:

$${\chi }_k=\{\begin{array}{l}0,k\le 1\\ n,k>1,\end{array}$$

(20)

We define the value of ${\Im }_k({\overrightarrow{x}}_{k-1})$ in a new way as:

$$\begin{gathered} {\Im }_k({\overrightarrow{x}}_{k-1})=S[x_{k-1}(\xi )]-\left(1-\frac{{\chi }_k}{n}\right)u^{\mu }\left({\displaystyle \sum _{l=0}^{m-1}{u}^{-\mu +l}{x}^{(l)}(0)+S[\varphi (\xi )]}\right) \\ +u^{\mu }S[R{x}_{k-1}+P_{k-1}], \end{gathered}$$

(21)

In the Equation (21), $P_k$ represents the homotopy polynomial [33] and given in the following form:

$$P_k=\frac{1}{\Gamma (k)}{\left[\frac{{\partial }^k}{\partial q^k}N\vartheta (\xi ,q)\right]}_{q=0},$$

(22)

and:

$$\vartheta (\xi ,q)={\vartheta }_0+q{\vartheta }_1+q^2{\vartheta }_2+\cdots .$$

(23)

Using Equation (21) in Equation (19), we get:

$$\begin{gathered} x_k(\xi )=({\chi }_k+\hslash )x_{k-1}(\xi )-\hslash \left(1-\frac{{\chi }_k}{n}\right)S^{-1}\left(u^{\mu }{\displaystyle \sum _{l=0}^{m-1}{u}^{-\mu +l}{x}^{(l)}(\xi )+u^{\mu }S[\psi (\xi )]}\right) \\ +\hslash S^{-1}\left[u^{\mu }S[R{x}_{k-1}+P_{k-1}]\right]. \end{gathered}$$

(24)

The novelty in our proposed approach is that a new correction function (24) is introduced by employing homotopy polynomials. Therefore, from Equation (24), we can calculate the distinct iterates $x_k(\xi )$ for $k>1$ and the q-HASTM series solution is given as follows:

$$x(\xi )={\displaystyle \sum _{k=0}^{\infty }{x}_k(\xi )}{\left(\frac{1}{n}\right)}^k$$

(25)

5. q-HASTM Solution for Nonlinear Fractional Dynamical Model of Interpersonal and Romantic Relationships

In this section, we find the solution of nonlinear fractional dynamical model of interpersonal and romantic relationship written in the following form:

$$\{\begin{array}{l}{D}_{\xi }^{\mu }x=-a_{1}x+b_{1}y(1-\epsilon y^2)+c_1,0<\mu \le 1\\ D_{\xi }^{\mu }y=-a_{2}y+b_{2}x(1-\epsilon x^2)+c_2,0<\mu \le 1,\end{array}$$

(26)

withinitial conditions:

$$\{\begin{array}{l}x(0)=0\\ y(0)=0.\end{array}$$

(27)

Employing the Sumudu transform operator on Equation (26) and using the initial conditions (27), we get the following results:

$$\begin{gathered} S[x(\xi )]-u^{\mu }{c}_1-u^{\mu }S[-a_{1}x+b_{1}y(1-\epsilon y^2)]=0, \\ S[y(\xi )]-u^{\mu }{c}_2-u^{\mu }S[-a_{2}y+b_{2}x(1-\epsilon x^2)]=0. \end{gathered}$$

(28)

Now, we write the nonlinear operator as follows:

$${\wp }_1[{\vartheta }_1(\xi ;q)]=S[{\vartheta }_1(\xi ;q)]-u^{\mu }{c}_1-u^{\mu }S[-a_1{\vartheta }_1(\xi ;q)+b_1{\vartheta }_2(\xi ;q)(1-\epsilon {\vartheta }_2^2(\xi ;q))],$$

(29)

and:

$${\wp }_2[{\vartheta }_2(\xi ;q)]=S[{\vartheta }_2(\xi ;q))]-u^{\mu }{c}_2-u^{\mu }S[-a_2{\vartheta }_2(\xi ;q)+b_2{\vartheta }_1(\xi ;q)(1-\epsilon {\vartheta }_1^2(\xi ;q))],$$

(30)

thus, we have:

$$\begin{gathered} {\Re }_{1k}({\overrightarrow{x}}_{k-1})=S[x_{k-1}(\xi )]-\left(1-\frac{{\chi }_k}{n}\right)u^{\mu }{c}_1 \\ -u^{\mu }S\left[-a_1{x}_{k-1}(\xi )+b_1{y}_{k-1}(\xi )-b_1\epsilon {\displaystyle \sum _{r=0}^{k-1}{y}_{k-1-r}\left({\displaystyle \sum _{j=0}^r{y}_j{y}_{r-j}}\right)}\right], \\ {\Re }_{2k}({\overrightarrow{y}}_{k-1})=S[y_{k-1}(\xi )]-\left(1-\frac{{\chi }_k}{n}\right)u^{\mu }{c}_2 \\ -u^{\mu }S\left[-a_2{y}_{k-1}(\xi )+b_2{x}_{k-1}(\xi )-b_2\epsilon {\displaystyle \sum _{r=0}^{k-1}{x}_{k-1-r}\left({\displaystyle \sum _{j=0}^r{x}_j{x}_{r-j}}\right)}\right]. \end{gathered}$$

(31)

The m-th–order deformation equations are given as follows:

$$\begin{gathered} S[x_k(\xi )-{\chi }_k{x}_{k-1}(\xi )]=\hslash {\Re }_{1k}({\overrightarrow{x}}_{k-1}), \\ \mathrm{and}S[y_k(\xi )-{\chi }_k{y}_{k-1}(\xi )]=\hslash {\Re }_{2k}({\overrightarrow{y}}_{k-1}). \end{gathered}$$

(32)

Next, on employing the inverse Sumudu operator, this gives:

$$\begin{gathered} x_k(\xi )={\chi }_k{x}_{k-1}(\xi )+\hslash S^{-1}[{\Re }_{1k}({\overrightarrow{x}}_{k-1})], \\ \mathrm{and}{y}_k(\xi )={\chi }_k{y}_{k-1}(\xi )+\hslash S^{-1}[{\Re }_{2k}({\overrightarrow{y}}_{k-1})]. \end{gathered}$$

(33)

Now using the initial approximation $x_0=0,y_0=0$ and recursive relation (33), we obtain the following iterations of the q-HASTM series solution:

$$x_1(\xi )=-c_1\hslash \frac{{\xi }^{\mu }}{\Gamma (1+\mu )},y_1(\xi )=-c_2\hslash \frac{{\xi }^{\mu }}{\Gamma (1+\mu )},$$

(34)

Following the same way the remaining components $x_k,y_k,k\ge 2$ of the q-HASTM solution can be easilyobtained, hence we obtain the complete solution. Finally, we arrive at the following series solution:

$$\begin{gathered} x(\xi )=\underset{M\to \infty }{\lim }{\displaystyle \sum _{k=0}^M{x}_k(\xi )}{\left(\frac{1}{n}\right)}^k \\ \mathrm{and}y(\xi )=\underset{M\to \infty }{\lim }{\displaystyle \sum _{k=0}^M{y}_k(\xi )}{\left(\frac{1}{n}\right)}^k. \end{gathered}$$

(35)

6. Basic Idea of ADM

To demonstrate the ADM solution procedure [34,35], we use a fractional nonlinear differential equation with the initial condition of the form:

$$D_{\xi }^{\mu }x(\xi )+Rx(\xi )+Nx(\xi )=\psi (\xi ),m-1<\mu \le m,$$

(36)

here $x$ is an unknown function of variable $\xi$ and $D_{\xi }^{\mu }=\frac{d^{\mu }}{d{\xi }^{\mu }}$ denotes the fractional operator of order $\mu$ proposed by Caputo, $m\in N,R$ indicates the bounded linear operator in $\xi$ and $N$ represents the general nonlinear differential operator in $\xi$.

Employing the operator $J_{\xi }^{\mu }$ on both sides of Equation (36) and making use of result (3), we get:

$$x(\xi )={\displaystyle \sum _{l=0}^{m-1}{\left(\frac{d^{l}x}{d{\xi }^l}\right)}_{\xi =0}}\frac{{\xi }^l}{l!}+J_{\xi }^{\mu }\psi (\xi )-J_{\xi }^{\mu }[Rx(\xi )+Nx(\xi )].$$

(37)

Further, we decompose the function $x(\xi )$ into sum of an infinite number of components expressed by the decomposition series as follows:

$$x={\displaystyle \sum _{n=0}^{\infty }{x}_n,}$$

(38)

and the nonlinear term can be decomposed as:

$$Nx={\displaystyle \sum _{n=0}^{\infty }{A}_n},$$

(39)

where $A_n$ denotes the Adomian polynomials, given as follows:

$$A_n=\frac{1}{n!}{\left[\frac{d^n}{d{\lambda }^n}N\left({\displaystyle \sum _{i=0}^n{\lambda }^i{x}_i}\right)\right]}_{\lambda =0},n=0,1,2,\dots$$

(40)

The components $x_0,x_1,x_2,\dots$ are determined recursively by substituting (38) and (39) into (37) leads to:

$${\displaystyle \sum _{n=0}^{\infty }{x}_n}={\displaystyle \sum _{l=0}^{m-1}{\left(\frac{d^{l}x}{d{\xi }^k}\right)}_{\xi =0}}\frac{{\xi }^l}{l!}+J_{\xi }^{\mu }\psi (\xi )-J_{\xi }^{\mu }\left[R\left({\displaystyle \sum _{n=0}^{\infty }{x}_n}\right)+{\displaystyle \sum _{n=0}^{\infty }{A}_n}\right].$$

(41)

This can be written as:

$$\begin{gathered} x_0+x_1+x_2+\cdots ={\displaystyle \sum _{l=0}^{m-1}{\left(\frac{d^{l}x}{d{\xi }^l}\right)}_{\xi =0}}\frac{{\xi }^l}{l!}+J_{\xi }^{\mu }\psi (\xi ) \\ -J_{\xi }^{\mu }\left[R(x_0+x_1+x_2+\cdots )+(A_0+A_1+A_2+\cdots )\right]. \end{gathered}$$

(42)

ADM uses the formal recursive relations as:

$$x_0={\displaystyle \sum _{l=0}^{m-1}{\left(\frac{d^{l}x}{d{\xi }^l}\right)}_{\xi =0}}\frac{{\xi }^l}{l!}+J_{\xi }^{\mu }\psi (\xi ),x_{n+1}=-J_{\xi }^{\mu }[R(x_n)+A_n],n\ge 0.$$

(43)

7. ADM Solution for Nonlinear Fractional Dynamical Model of Interpersonal and Romantic Relationship

To solve the nonlinear fractional dynamical model of interpersonal and romantic relations (26), we apply the operator $J_{\xi }^{\mu }$ on both sides of Equation (26) and use result (3) to obtain:

$$\begin{gathered} x={\displaystyle \sum _{l=0}^{1-1}\frac{{\xi }^l}{l!}{\left[D_{\xi }^{l}x\right]}_{\xi =0}}+J_{\xi }^{\mu }{c}_1+J_{\xi }^{\mu }\left[-a_{1}x+b_{1}y(1-\epsilon y^2)\right], \\ y={\displaystyle \sum _{l=0}^{1-1}\frac{{\xi }^l}{l!}{\left[D_{\xi }^{l}y\right]}_{\xi =0}}++J_{\xi }^{\mu }{c}_2+J_{\xi }^{\mu }\left[-a_{2}y+b_{1}x(1-\epsilon x^2)\right]. \end{gathered}$$

(44)

This yields the following recursive relations using Equation (43):

$$x_0={\displaystyle \sum _{l=0}^0\frac{{\xi }^l}{l!}{\left[D_{\xi }^{l}x\right]}_{\xi =0}+J_{\xi }^{\mu }{c}_1},y_0={\displaystyle \sum _{l=0}^0\frac{{\xi }^l}{l!}{\left[D_{\xi }^{l}y\right]}_{\xi =0}+J_{\xi }^{\mu }{c}_2}.$$

(45)

$$\begin{gathered} x_{n+1}=J_{\xi }^{\mu }[-a_1{x}_n+b_1{y}_n-\epsilon b_1{A}_n],n=0,1,2,\dots \\ {y}_{n+1}=J_{\xi }^{\mu }[-a_2{y}_n+b_2{x}_n-\epsilon b_2{A}^{\prime }{}_n], \end{gathered}$$

(46)

where:

$${\displaystyle \sum _{n=0}^{\infty }{A}_n}=y^3,{\displaystyle \sum _{n=0}^{\infty }{A}^{\prime }{}_n}=x^3.$$

(47)

The first few components of the Adomian polynomials are given as follows:

$$\begin{gathered} A_0=y_0^3,A_1=3y_0^2{y}_1,A_2=3\left[y_0^2{y}_2+y_0{y}_1^2\right], \\ {A}_0^{\prime }=x_0^3,A_1^{\prime }=3x_0^2{x}_1,A_2^{\prime }=3\left[x_0^2{x}_2+x_0{x}_1^2\right], \end{gathered}$$

(48)

The components of the ADM solution can be simply obtained by making use of the above recursive relation:

$$\begin{gathered} x_0=c_1\frac{1}{\Gamma (1+\mu )}{\xi }^{\mu }, \\ {y}_0=c_2\frac{1}{\Gamma (1+\mu )}{\xi }^{\mu }, \\ {x}_1=\frac{-a_1{c}_1}{\Gamma (1+2\mu )}{\xi }^{2\mu }+\frac{b_1{c}_2}{\Gamma (1+2\mu )}{\xi }^{2\mu }-\frac{\epsilon b_1{c}_2^3}{{\left(\Gamma (1+\mu )\right)}^3}\frac{\Gamma (1+3\mu )}{\Gamma (1+4\mu )}{\xi }^{4\mu }, \\ {y}_1=\frac{-a_2{c}_2}{\Gamma (1+2\mu )}{\xi }^{2\mu }+\frac{b_2{c}_1}{\Gamma (1+2\mu )}{\xi }^{2\mu }-\frac{\epsilon b_2{c}_1^3}{{\left(\Gamma (1+\mu )\right)}^3}\frac{\Gamma (1+3\mu )}{\Gamma (1+4\mu )}{\xi }^{4\mu }, \end{gathered}$$

(49)

and so on. In this way the remaining components of the ADM solution can be determined. Hence, the ADM solution of Equation (26) is given as follows:

$$\begin{gathered} x(\xi )=c_1\frac{1}{\Gamma (1+\mu )}{\xi }^{\mu }+\frac{-a_1{c}_1}{\Gamma (1+2\mu )}{\xi }^{2\mu }+\frac{b_1{c}_2}{\Gamma (1+2\mu )}{\xi }^{2\mu } \\ -\frac{\epsilon b_1{c}_2^3}{{\left(\Gamma (1+\mu )\right)}^3}\frac{\Gamma (1+3\mu )}{\Gamma (1+4\mu )}{\xi }^{4\mu }+\dots \\ y(\xi )=c_2\frac{1}{\Gamma (1+\mu )}{\xi }^{\mu }+\frac{-a_2{c}_2}{\Gamma (1+2\mu )}{\xi }^{2\mu }+\frac{b_2{c}_1}{\Gamma (1+2\mu )}{\xi }^{2\mu } \\ -\frac{\epsilon b_2{c}_1^3}{{\left(\Gamma (1+\mu )\right)}^3}\frac{\Gamma (1+3\mu )}{\Gamma (1+4\mu )}{\xi }^{4\mu }+\dots \end{gathered}$$

(50)

which can be recovered from the q-HASTM solution by setting $\hslash =-1$ and $n=1.$ However, most of the results obtained by using other analytical schemes such as HPM, VIM, DTM, and ADM converge to the corresponding numerical results in a very small region. Instead of these methods, in the case of q-HASTM, we can easily settle and restrict the region of convergence of series solution obtained by q-HASTM by setting appropriate values of $\hslash$ and $n.$

8. Numerical Simulations

In this section, we perform the numerical simulations for $x(\xi )$ and at the distinct fractional Brownian motions $\mu =0.95$ and $\mu =0.90$ and also for the standard motion $\mu =1.$ The special solution of fractional dynamical model of the problem is obtained by employing q-HASTM and ADM with different parameters as $a_1=0.05,b_1=0.04,c_1=0.2,a_2=0.07,b_2=0.06,c_2=0.3$ and $\epsilon =0.01.$ The results are shown through

Table 1. Comparative study between q-HASTM and ADM for $x$ when $\mu$ = 1, $a_1=0.05,b_1=0.04,c_1=0.2,a_2=0.07,b_2=0.06,c_2=0.3$ and $\epsilon =0.01.$

$\mathit{\xi }$	ADM	q-HASTM(for $\mathit{n}\mathbf{=}\mathbf{1}$ and $\mathit{\hslash }\mathbf{=}\mathbf{-}\mathbf{1}$)	q-HASTM(for $\mathit{n}\mathbf{=}\mathbf{2}$ and $\mathit{\hslash }\mathbf{=}\mathbf{-}\mathbf{1.99}$)
0	0	0	0
2	0.4033806333	0.4033788000	0.4033790454
4	0.8110188698	0.8109674666	0.8109706878
6	1.219148908	1.218802800	1.218817897
8	1.624049986	1.622732800	1.622778592
10	2.022138493	2.018416666	2.018525638
12	2.410163488	2.401324800	2.401546847
14	2.785587314	2.766738800	2.767144981
16	3.147219587	3.109751466	3.110437747
18	3.496167079	3.425266800	3.426357807
20	3.837175265	3.708000000	3.709652760

and

Table 2. Comparative study between q-HASTM and ADM for $y$ when $\mu$ = 1,$a_1=0.05,b_1=0.04,c_1=0.2,a_2=0.07,b_2=0.06,c_2=0.3$ and $\epsilon =0.01.$

$\mathit{\xi }$	ADM	q-HASTM(for $\mathit{n}\mathbf{=}\mathbf{1}$ and $\mathit{\hslash }\mathbf{=}\mathbf{-}\mathbf{1}$)	q-HASTM(for $\mathit{n}\mathbf{=}\mathbf{2}$ and $\mathit{\hslash }\mathbf{=}\mathbf{-}\mathbf{1.99}$)
0	0	0	0
2	0.5829289274	0.5829274000	0.5829287007
4	1.134891434	1.134838400	1.134860301
6	1.659550026	1.659119400	1.659232190
8	2.159330141	2.157414400	2.157773939
10	2.635727526	2.629625000	2.630507313
12	3.089599578	3.073910400	3.075746263
14	3.521399539	3.486687400	3.490096920
16	3.931326626	3.862630400	3.868457617
18	4.319398257	4.194671400	4.204018861
20	4.685515415	4.474000000	4.488263354

and Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11. From Table 1 and Table 2, it can be noticed that the values of the approximate solution at distinct grid points obtained by the q-HASTM and ADM are in a very good agreement. Figure 1 and Figure 2 depict the behavior of $x(\xi )$ and $y(\xi )$ for the different values of $\mu .$ From Figure 1 we can observe that when we decrease the values of $\mu$ in the fractional dynamic model then the love of individual 1 for his/her partner decreases. From Figure 2 we can observe that when we decrease the value of $\mu$ in the fractional dynamic model then the love of individual 2 for his/her partner decreases. Figure 3, Figure 4 and Figure 5 reveals the behavior of $x(\xi )$ and $y(\xi )$ for the different value of $\mu$ at $\hslash =-1$ and $n=1.$ From Figure 3, Figure 4 and Figure 5, we see that when we decrease the values of $\mu$ in the fractional dynamic model then romantic relation between the couple decreases. Figure 6, Figure 7 and Figure 8 demonstrate the $\hslash$-curves for $x(\xi )$ at distinct values of $\mu$ and $n.$ Figure 9, Figure 10 and Figure 11 show the $\hslash$-curves of $y(\xi )$ at distinct values of $\mu$ and $n.$ The horizontal line segment in $\hslash$-curves shows the range of convergence of q-HASTM solution. From the Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, it can be noticed that the range of convergence is directly proportional to the value of asymptotic parameter n. Hence, the results obtained by using the proposed technique converge very fast as compared to other existing analytical sachems.

9. Conclusions

In this article, we have proposed a novel numerical method named q-HASTM to solve nonlinear dynamical models of interpersonal and romantic relationships for marriages with a fractional approach and compared the results with ADM. The numerical and graphical results reveal the successfully application of q-HASTM for solving the fractional dynamical model. The results derived with the aid of both the techniques are in an excellent agreement. The displacement shows a new nature for the time fractional derivative compared to the integer order derivative. The convergence of the q-HASTM solution can be adjusted and controlled with the aid of the auxiliary parameter $\hslash$ and asymptotic parameter n. The results presented in the form of graphs and $\hslash$-curves indicate that the suggested scheme is very efficient and accurate. Hence, it can be concluded that the proposed approach is highly logical and can be employed to investigate a wide range of nonlinear mathematical models of fractional order appearing in real world problems. Moreover, q-HASTM opens new doors in the fields of mathematical modeling and fractional calculus.