
Artificial intelligence (AI) and generative AI (GenAI) have sparked confusion and concern regarding their impact on education. Beyond the assessment integrity risks that currently draw the most attention, technologies such as ChatGPT, Copilot, and Gemini have also been identified as tools that can support learning. Project work, especially when there is no single correct solution, provides a great opportunity for integration, fostering technology knowledge and higher learning standards. However, no AI-integration framework for project-based work is available, resulting in a limited understanding of how AI integration can occur or be maximized. To address this, a collaborative effort of 16 educators from 9 Australian universities has led to the development of a generic AI implementation framework, built upon the CDIO approach. With a focus on engineering education, this framework can be adapted to other project-based learning contexts, where educators can pick and choose the relevant implementation items as needed. This framework is called the Project-work Artificial Intelligence Integration Framework (PAIIF), and its development and structure are outlined here. Initial implementations have shown the effectiveness of promoting reflection and guidance on where and how AI integration can occur.
Citation: Sasha Nikolic, Zach Quince, Anna Lidfors Lindqvist, Peter Neal, Sarah Grundy, May Lim, Faham Tahmasebinia, Shannon Rios, Josh Burridge, Kathy Petkoff, Ashfaque Ahmed Chowdhury, Wendy S.L. Lee, Rita Prestigiacomo, Hamish Fernando, Peter Lok, Mark Symes. Project-work Artificial Intelligence Integration Framework (PAIIF): Developing a CDIO-based framework for educational integration[J]. STEM Education, 2025, 5(2): 310-332. doi: 10.3934/steme.2025016
[1] | Hardik Joshi, Brajesh Kumar Jha, Mehmet Yavuz . Modelling and analysis of fractional-order vaccination model for control of COVID-19 outbreak using real data. Mathematical Biosciences and Engineering, 2023, 20(1): 213-240. doi: 10.3934/mbe.2023010 |
[2] | Shina D. Oloniiju, Olumuyiwa Otegbeye, Absalom E. Ezugwu . Investigating the impact of vaccination and non-pharmaceutical measures in curbing COVID-19 spread: A South Africa perspective. Mathematical Biosciences and Engineering, 2022, 19(1): 1058-1077. doi: 10.3934/mbe.2022049 |
[3] | Emmanuel Addai, Mercy Ngungu, Musibau Abayomi Omoloye, Edmore Marinda . Modelling the impact of vaccination and environmental transmission on the dynamics of monkeypox virus under Caputo operator. Mathematical Biosciences and Engineering, 2023, 20(6): 10174-10199. doi: 10.3934/mbe.2023446 |
[4] | Francesco Salvarani, Gabriel Turinici . Optimal individual strategies for influenza vaccines with imperfect efficacy and durability of protection. Mathematical Biosciences and Engineering, 2018, 15(3): 629-652. doi: 10.3934/mbe.2018028 |
[5] | Simphiwe M. Simelane, Justin B. Munyakazi, Phumlani G. Dlamini, Oluwaseun F. Egbelowo . Projections of human papillomavirus vaccination and its impact on cervical cancer using the Caputo fractional derivative. Mathematical Biosciences and Engineering, 2023, 20(7): 11605-11626. doi: 10.3934/mbe.2023515 |
[6] | Yujie Sheng, Jing-An Cui, Songbai Guo . The modeling and analysis of the COVID-19 pandemic with vaccination and isolation: a case study of Italy. Mathematical Biosciences and Engineering, 2023, 20(3): 5966-5992. doi: 10.3934/mbe.2023258 |
[7] | Muntaser Safan . Mathematical analysis of an SIR respiratory infection model with sex and gender disparity: special reference to influenza A. Mathematical Biosciences and Engineering, 2019, 16(4): 2613-2649. doi: 10.3934/mbe.2019131 |
[8] | Jutarat Kongson, Chatthai Thaiprayoon, Apichat Neamvonk, Jehad Alzabut, Weerawat Sudsutad . Investigation of fractal-fractional HIV infection by evaluating the drug therapy effect in the Atangana-Baleanu sense. Mathematical Biosciences and Engineering, 2022, 19(11): 10762-10808. doi: 10.3934/mbe.2022504 |
[9] | Xiaojing Wang, Yu Liang, Jiahui Li, Maoxing Liu . Modeling COVID-19 transmission dynamics incorporating media coverage and vaccination. Mathematical Biosciences and Engineering, 2023, 20(6): 10392-10403. doi: 10.3934/mbe.2023456 |
[10] | Ahmed Alshehri, Saif Ullah . A numerical study of COVID-19 epidemic model with vaccination and diffusion. Mathematical Biosciences and Engineering, 2023, 20(3): 4643-4672. doi: 10.3934/mbe.2023215 |
Artificial intelligence (AI) and generative AI (GenAI) have sparked confusion and concern regarding their impact on education. Beyond the assessment integrity risks that currently draw the most attention, technologies such as ChatGPT, Copilot, and Gemini have also been identified as tools that can support learning. Project work, especially when there is no single correct solution, provides a great opportunity for integration, fostering technology knowledge and higher learning standards. However, no AI-integration framework for project-based work is available, resulting in a limited understanding of how AI integration can occur or be maximized. To address this, a collaborative effort of 16 educators from 9 Australian universities has led to the development of a generic AI implementation framework, built upon the CDIO approach. With a focus on engineering education, this framework can be adapted to other project-based learning contexts, where educators can pick and choose the relevant implementation items as needed. This framework is called the Project-work Artificial Intelligence Integration Framework (PAIIF), and its development and structure are outlined here. Initial implementations have shown the effectiveness of promoting reflection and guidance on where and how AI integration can occur.
Mathematical forecasts are yet to gain its improvement in fighting the COVID-19 pandemic. The meteoric spread, the course of actions, and the methods concerned in the containment of the infection demand the soonest and quickest understanding in perceiving solutions in line along with the biological, physiological, environmental, and continual aspects of life with much better predictions by mathematical computerized modeling. Almost all the countries are still suffering from the infection and its destructiveness is ascending gradually. The disease has spread, and is spreading amongst different countries or regions mainly through air-travel mode [1]. The infection doesn't happen on its own and the present occurrence is possibly because of our ignorance or casual thinking and way of living [2]. As a matter of fact, there are many looming queries about COVID-19 spread, such as "number of individuals that will be infected in coming days", "time of appearance of rate of inflexion point of infection", "infections during the pinnacle period and number of individuals that are infected at that time", "effectiveness of present action to control the spread", etc. Usually, in analysis of epidemiology, virus causing diseases grow exponentially maintaining a determined rate of reproduction [3].
One of the most economical and efficient means to prevent and control the infection spread is mass vaccination. The world is now getting leaky vaccinated to reduce the infection as the vaccines obstruct signs and symptoms of lethal infections but do not put a stop to contamination or onward spread [4,5,6,7,8,9,10,11]. It is seen in India that the government is facing a huge challenge to attain preferable vaccination coverage in order for group immunity to be in effect [12]. Hysteresis loops can surface in social replica of dynamical behaviour of vaccination [13]. The impact of varying vaccine efficacies has been investigated in response to control the spread of disease [14,15,16]. The vaccine is effective (infection can be prevented) for η=1 (here we denote η as vaccine efficacy), whilst the value η=0 means the vaccine is ineffective to prevent the virus infection. In the bio-mathematics regime, fractional calculus is now one of the vital studies to elaborate precise cellular behavior, and many human disease's mathematical models have been described by various types of fractional-order derivative operators [17,18,19,20].
FDEs (short for Fractional Differential Equation) are used to include the memory of the phenomenon in the system's dynamics. As per literature [21], there are multiple kinds of fractional derivative operators, amongst which we present the study of the model by fractional differential equation by CF (short for Caputo-Fabrizio derivative) and ABC (short for Atangana-Baleanu-Caputo derivative) sense. ODEs (short for Ordinary Differential Equation) are particular cases of FDEs. FDEs provide descriptions of the courses of endemics, epidemics, pandemics, and some other natural processes in a better way than the ODEs [22]. In reference to deterministic models of the contagiousness of any disease, ODE models are well-structured yet sometimes researchers use models using FDE [23,24]. Some recent research papers focusing on FDE for various real-life problems which are modelled mathematically to understand the problems and their solutions has been studied [25,26,27,28,29,30,31,32,33,34,35,36].
As per the author's knowledge, the remark and application for vaccine efficacy, and the study and experimentation has not yet looked into comparison between the mentioned fractional derivative operators for the SEVR COVID-19 model with the memory effect. The quantitative analysis makes the two operators propitious choices in modeling the dynamics of the disease mathematically. In the upcoming segments, the model is formed with fractional differential equation by CF and ABC sense. This technique is used to find the course of action of the infection whilst vaccination program is on the go. Using Next Generation Matrix Method the expression for the basic reproduction number, Rn, is found and by dint of that local stability are investigated for the equilibrium points. The existence and uniqueness for the models are done using the Banach fixed-point theorem. Considering the course of events for beginning of the transmission to clinical recuperation and seriousness of illness, the appropriateness of utilizing the two operators, joined by relative similitudes, are examined just as their disparities. Lastly, the article is wrapped-up by scrutinizing the seriousness of the spread of infection in terms of CF and ABC, and depict their upcoming involvement for the growth of infection control policy at worldwide level to limit the spread amongst humans.
The follow-up sections deal with the essence of the manuscript. In Section 2, some preliminary requisite knowings of CF and ABC non-singular kernel fractional derivatives is done to get acquainted with customary notations. This section also has model formulation with a box diagram to have a glimpse on the transmission dynamics. The parameters are broadly described with proper justification. The equilibrium points are evaluated and also Rn is calculated. The stability analysis for the infection-free equilibrium has been done in regard to Rn and Jacobian matrix. In Section 3, the solutions of the CF and ABC systems are shown to exist and unique using Banach fixed-point theory and Lipschitz condition. Section 4 deals with discussion and numerical simulation of both the systems in terms of graphs, order of fractional derivative and Rn. At last, in Section 5, we draw the conclusions of the entire manuscript.
Here, the preliminary and primary definitions in terms of properties of the mentioned two fractional derivative operators with exponential and Mittag-Leffler kernels are presented.
Definition 1. [37] Consider f∈H1(c1,c2), c2>c1, α∈]0,1[ then the CF fractional derivative is defined as:
CFc1Dαt[f(t)]=N(α)1−α∫tc1f′(x)e[−αt−x1−α]dx, |
where N(α) is a normalization function satisfying N(0)=N(1)=1 if f∉H1(c1,c2) then
Dαt(f(t))=αN(α)1−α∫tc1(f(t)−f(x))e[−αt−x1−α]dx. |
Definition 2. [37,38] Let us assume that 0<α<1, then the integral of the fractional order α for a function f(t) is defined as:
Iαt(f(t))=2(1−α)(2−α)N(α)f(t)+2α(2−α)N(α)∫t0f(y)dy, |
t≥0.
Definition 3. [39] Consider f∈H1(c1,c2), c2>c1, α∈]0,1[ then the ABC is defined as:
ABCc1Dαt[f(t)]=N(α)1−α∫tc1f′(x)Mα[−α(t−x)α1−α]dx, |
where Mα represent the Mittag-Leffler function and N(α) is a normalization function satisfying N(0)=N(1)=1 and N(α)=1−α+αΓ(α).
Definition 4. [39] Let us assume that 0<α<1, then the integral of the ABC, fractional order α for a function f(t) is defined as:
ABCc1Iαt[f(t)]=1−αN(α)f(t)+αN(α)Γ(α)×∫tc1(t−y)α−1f(y)dy, |
t≥0.
It is important to identify and then study the dynamical behaviour of unwanted results of COVID-19 disease on the populace growth and to predict what may occur. Keeping track of the infection has been a steadily tiring job and vaccination has been a usual approach for lessening the infection, day-to-day vaccination is currently supplied in all nations. It has been believed that the immune system will generate antibody against the virus after the vaccine doses. Even so, it may be not in a fully protective level. The memory effect is an important element in several biological phenomenons, which is included in the models. Hence, it is practicable to use fractional alternation, as fractional derivatives essentially involve the effect of memory. One of the best ways to tackle such issue is to put the fractional derivatives of order α∈(0,1) for the state variables. The next step is to have some modification in the biological parameters such that for a FDE system the both sides of an equation must have dimension 1(time)α. This can be done by raising each of the parameters by α [40]. Motivated by these details, the SEV model presented Liu et al. [42] and Abuasad et al. [41], the updated SEVR model with vaccination by CF and ABC non-singular kernel fractional derivative operators are as follows:
CF0DαtS(t)=λα−γαSE−(dα+να)S, | (2.1) |
CF0DαtE(t)=γαSE+(1−η)γαVE−(dα+βα)E, | (2.2) |
CF0DαtV(t)=ναS−dαV−(1−η)γαVE, | (2.3) |
CF0DαtR(t)=βαE−dαR. | (2.4) |
and
ABC0DαtS(t)=λα−γαSE−(dα+να)S, | (2.5) |
ABC0DαtE(t)=γαSE+(1−η)γαVE−(dα+βα)E, | (2.6) |
ABC0DαtV(t)=ναS−dαV−(1−η)γαVE, | (2.7) |
ABC0DαtR(t)=βαE−dαR. | (2.8) |
The initial conditions for both the systems are S(0)=S0>0, E(0)=E0≥0, V(0)=V0≥0, and R(0)=R0≥0. The virus spread method is considered amongst humans. In the models (2.1)–(2.4) and (2.5)–(2.8), S(t) denotes the fractions of susceptible individuals, E(t) denotes the fractions of virus effected individuals, V(t) denotes fractions of vaccinated individuals, and R(t) denotes fractions of individuals who have recovered from the disease, at any time t. The explanation for the parameters (non-state variables) is as following:
i) The parameter values involved in the aforementioned systems are considered to be non-negative, and λ>0.
ii) λ: rate at which individuals are considered to be born susceptible.
iii) d: death rate of all the state variables.
iv) ν: rate at which susceptibles are vaccinated.
v) γ: rate at which susceptibles are effected by the virus.
vi) (1−η)γ: rate at which vaccinated individuals are affected, which can be possible due to findings present in [43].
v) β: recovery rate of affected individuals but are open to infection again if proper medical and government guidelines are not followed.
The vaccine efficacy can be given a numerical value [44], and as per our knowledge the vaccine efficacy has not yet been modelled to any mathematical COVID-19 model. So for the solutions of system (2.1)–(2.4) and (2.5)–(2.8) we take two values of η in order to compare the solution trajectories. Next, one of the values for death rate of state variables, d, is taken as 0.002 to have a view at the infection in terms of the basic reproduction number. The values of fractional order α are taken arbitrarily. This section is followed by two other subsections dealing with finding the equilibrium points and the basic reproduction number along with stability analysis of the equilibrium points.
Here, only two equilibrium points of the models given by the system (2.1)–(2.4) and (2.5)–(2.8) are presented with biological clarification. The points are as following:
i) X0(λαdα+να,0,λαναdα(dα+να),0), which is named as "Infection-free equilibrium (IFE) point". This is the case where affected individuals are zero and individuals are in infection free state. The vaccination for this point is fruitful for both short and long run as the infection can be contained.
ii) X∗(S∗,E∗,V∗,R∗), which is named as "Co-existing equilibrium (CoE) point". The point X∗ satisfies the below mentioned system of the equations.
λα−γαS∗E∗−(dα+να)S∗=0,γαS∗E∗+(1−η)γαV∗E∗−(dα+βα)E∗=0,ναS∗−dαV∗−(1−η)γαV∗E∗=0,βαE∗−dαR∗=0. |
In this case, the affected individuals co-exist with susceptible individuals and this can be an alarming issue for the infection spread. The vaccination cannot be applied to the affected compartment for this situation as it will wreck havoc to the immune system, and also rate of positive infection will rise leading to indefinite death toll.
It is important information in the mathematical regime of infectious disease models as it helps in finding whether a pandemic will happen or not. The number Rn resembles the expected cases of infections caused by an infectious individual inside susceptible populace. The next-generation matrix method [45] is used to find Rn of the model presented in this article.
The Jacobian matrices F and V with the rate of appearance of new infections and rate of transfer of individuals into and out of compartments, respectively, at X0 are given by
F=(00000γαλα{dα+(1−η)να}dα(dα+να)00να0000βα00),V=(dα+ναγαλαdα+να000dα+βα000(1−η)γαλαναdα(dα+να)dα0000dα). |
Proceeding further, we get
FV−1=diag(0,Rn,0,0). |
The eigenvalues for the diagonal matrix FV−1 are the diagonal elements itself, and amongst them the dominant eigenvalue is Rn=γαλα{dα+(1−η)να}dα(dα+βα)(dα+να), which is the reproduction number for the proposed models.
The infection will self-annihilate if Rn<1, otherwise for Rn>1 the infection will conquer and begin to a pandemic or an epidemic until and unless the medical as well as government guidelines are not implemented amongst the mass.
The discussion will be for the case Rn<1, then the system (2.1)–(2.4) and (2.5)–(2.8) has IFE point X0 and it is globally asymptotically stable in an invariant region [41,42]
Γ={(S,E,V,R):S>0,E≥0,V≥0,R≥0,0<S+E+V+R≤1}. |
Now, to check local stablity the Jacobian of the models (2.1)–(2.4) or (2.5)–(2.8) at X0 is evaluated and it is given by the matrix
J(X0)=F(X0)−V(X0)=(−(dα+να)−γαλαdα+να000(dα+βα)(Rn−1)00να−(1−η)γαλαναdα(dα+να)−dα00βα0−dα). |
The characteristic equation for the matrix J(X0) in variable ξ is given by
4∏i=1(aii−ξ)=0, |
where a11=−(dα+να), a22=(dα+βα)(Rn−1), a33=−dα, and a44=−dα are the eigenvalues. Therefore, the IFE point X0 of the models (2.1)–(2.4) and (2.5)–(2.8) is locally asymptotically stable for a22<0, implies that Rn<1. So, it can be stated in a theorem as following:
Theorem 1. For Rn<1, the IFE point X0 is locally asymptotically stable for the COVID-19 models (2.1)–(2.4) and (2.5)–(2.8).
The follow-up section is about two of the most essential parts of any system of equations viz, existence and uniqueness of solutions.
Here, in this section, necessary theorems are provided and discussed thoroughly to support the existence and uniqueness of solutions of models (2.1)–(2.4) and (2.5)–(2.8) together with the initial conditions, by applying the Banach fixed-point theory.
The findings for the system of equations of COVID-19 CF model using the fractional integral are as follows [38,46]:
S(t)=CF0Iαt{λα−γαSE−(dα+να)S}+S0,E(t)=CF0Iαt{γαSE+(1−η)γαVE−(dα+βα)E}+E0,V(t)=CF0Iαt{ναS−dαV−(1−η)γαVE}+V0,R(t)=CF0Iαt{βαE−dαR}+R0. |
The symbols, after arranged in the view of the problem, gives the following results:
S(t)=S0+2(1−α)(2−α)N(α)[λα−γαS(t)E(t)−(dα+να)S(t)]+2α(2−α)N(α)∫t0[λα−γαS(y)E(y)−(dα+να)S(y)]dy,E(t)=E0+2(1−α)(2−α)N(α)[γαS(t)E(t)+(1−η)γαV(t)E(t)−(dα+βα)E(t)]+2α(2−α)N(α)∫t0[γαS(y)E(y)+(1−η)γαV(y)E(y)−(dα+βα)E(y)]dy,V(t)=V0+2(1−α)(2−α)N(α)[ναS(t)−dαV(t)−(1−η)γαV(t)E(t)]+2α(2−α)N(α)∫t0[ναS(y)−dαV(y)−(1−η)γαV(y)E(y)]dy,R(t)=R0+2(1−α)(2−α)N(α)[βαE(t)−dαR(t)]+2α(2−α)N(α)∫t0[βαE(y)−dαR(y)]dy. |
For the sake of simplicity, we consider
H1(S,t)=[λα−γαS(t)E(t)−(dα+να)S(t)], | (3.1) |
H2(E,t)=[γαS(t)E(t)+(1−η)γαV(t)E(t)−(dα+βα)E(t)], | (3.2) |
H3(V,t)=[ναS(t)−dαV(t)−(1−η)γαV(t)E(t)], | (3.3) |
H4(R,t)=[βαE(t)−dαR(t)]. | (3.4) |
Theorem 2. The Lipschitz condition and contraction are justified by the kernel
i) H1 if 0≤(γαω2+dα+να)<1.
ii) H2 if 0≤(γαω1+(1−η)γαω3+dα+βα)<1.
iii) H3 if 0≤(dα+(1−η)γαω2)<1.
iv) H4 if 0≤dα<1.
Proof. Let us assume that the state variables present in the COVID-19 models are bounded, and therefore it can be said that ∥S(t)∥≤ω1, ∥E(t)∥≤ω2, ∥V(t)∥≤ω3, and ∥R(t)∥≤ω4 [46].
Let us consider S(t) and S′(t) be the two functions, and therefore
∥H1(S,t)−H1(S′,t)∥=∥−γα(S(t)−S′(t))E(t)−(dα+να)(S(t)−S′(t))∥≤∥γαE(t)(S(t)−S′(t))∥+∥(dα+να)(S(t)−S′(t))∥≤[γα∥E(t)∥+(dα+να)]∥S(t)−S′(t)∥≤[γαω2+dα+να]∥S(t)−S′(t)∥. |
So, it can be mentioned now that
∥H1(S,t)−H1(S′,t)∥≤ϕ1∥S(t)−S′(t)∥, |
where ϕ1=(γαω2+dα+να). Thus, the Lipschitz condition is achieved for kernel H1 and 0≤(γαω2+dα+να)<1 [46] gives H1, and also contraction is justified.
Like wise, contraction and the Lipschitz condition can be done and written for the kernels H2, H3, and H4 as mentioned below;
∥H2(E,t)−H2(E′,t)∥≤ϕ2∥E(t)−E′(t)∥,∥H3(V,t)−H3(V′,t)∥≤ϕ3∥V(t)−V′(t)∥,∥H4(R,t)−H4(R′,t)∥≤ϕ4∥R(t)−R′(t)∥, |
where ϕ2=(γαω1+(1−η)γαω3+dα+βα), ϕ3=(dα+(1−η)γαω2), and ϕ4=dα.
Taking the kernels H1, H2, H3, and H4 into consideration the following equations are obtained, which are as follows:
S(t)=S0+2(2−α)N(α)[(1−α)H1(S,t)+α∫t0H1(S,z)dz], | (3.5) |
E(t)=E0+2(2−α)N(α)[(1−α)H2(E,t)+α∫t0H2(E,z)dz], | (3.6) |
V(t)=V0+2(2−α)N(α)[(1−α)H3(V,t)+α∫t0H3(V,z)dz], | (3.7) |
R(t)=R0+2(2−α)N(α)[(1−α)H4(R,t)+α∫t0H4(R,z)dz]. | (3.8) |
The following iterative formulae are now seen closely, which are as follows:
Sn(t)=S0+2(2−α)N(α)[(1−α)H1(Sn−1,t)+α∫t0H1(Sn−1,z)dz], | (3.9) |
En(t)=E0+2(2−α)N(α)[(1−α)H2(En−1,t)+α∫t0H2(En−1,z)dz], | (3.10) |
Vn(t)=V0+2(2−α)N(α)[(1−α)H3(Vn−1,t)+α∫t0H3(Vn−1,z)dz], | (3.11) |
Rn(t)=R0+2(2−α)N(α)[(1−α)H4(Rn−1,t)+α∫t0H4(Rn−1,z)dz], | (3.12) |
with the initial conditions of the COVID-19 model given in the system of Eqs (2.1)–(2.4). The upcoming expressions for the difference of sequential terms are written as:
ρn(t)=Sn(t)−Sn−1(t)=2(1−α)(2−α)N(α)[(H1(Sn−1,t)−H1(Sn−2,t))]+2α(2−α)N(α)∫t0(H1(Sn−1,z)−H1(Sn−2,z))dz,σn(t)=En(t)−En−1(t)=2(1−α)(2−α)N(α)[(H2(En−1,t)−H2(En−2,t))]+2α(2−α)N(α)∫t0(H2(En−1,z)−H2(En−2,z))dz,δn(t)=Vn(t)−Vn−1(t)=2(1−α)(2−α)N(α)[(H3(Vn−1,t)−H3(Vn−2,t))]+2α(2−α)N(α)∫t0(H3(Vn−1,z)−H3(Vn−2,z))dz,πn(t)=Rn(t)−Rn−1(t)=2(1−α)(2−α)N(α)[(H4(Rn−1,t)−H4(Rn−2,t))]+2α(2−α)N(α)∫t0(H4(Rn−1,z)−H4(Rn−2,z))dz. |
It is evident that
Sn(t)=n∑m=1ρm(t),En(t)=n∑m=1σm(t),Vn(t)=n∑m=1δm(t),Rn(t)=n∑m=1πm(t). |
Involving norm to ρn(t) then using triangle inequality, we get
∥ρn(t)∥=∥Sn(t)−Sn−1(t)∥≤2(1−α)(2−α)N(α)∥(H1(Sn−1,t)−H1(Sn−2,t))∥+2α(2−α)N(α)∥∫t0(H1(Sn−1,z)−H1(Sn−2,z))dz∥. |
It has been shown that the the Lipschitz condition is satisfied by the kernel H1, therefore
∥Sn(t)−Sn−1(t)∥≤2ϕ1(2−α)N(α)[(1−α)∥Sn−1−Sn−2∥+α∫t0∥Sn−1−Sn−2∥dz]. |
Thus, it can be deduced that
∥ρn(t)∥≤2ϕ1(2−α)N(α)[(1−α)∥ρn−1(t)∥+α∫t0∥ρn−1(z)∥dz]. | (3.13) |
In the same way, other results are followed as;
∥σn(t)∥≤2ϕ2(2−α)N(α)[(1−α)∥σn−1(t)∥+α∫t0∥σn−1(z)∥dz], | (3.14) |
∥δn(t)∥≤2ϕ3(2−α)N(α)[(1−α)∥δn−1(t)∥+α∫t0∥δn−1(z)∥dz], | (3.15) |
∥πn(t)∥≤2ϕ4(2−α)N(α)[(1−α)∥πn−1(t)∥+α∫t0∥πn−1(z)∥dz]. | (3.16) |
By taking the inequalities (3.13)–(3.16) into consideration the following theorem can be stated, which is:
Theorem 3. The system (2.1)–(2.4) has a solution if there exist a real number τ0, satisfying
2ϕi(2−α)N(α)[1+α(τ0−1)]<1, |
for all i=1,2,3,4.
Proof. Inequations (3.13)–(3.16) are considered along with the fact that the state variables are bounded and the kernels Hi (i=1,2,3,4) satisfy the Lipschitz condition. The following mentioned relations employ the iterative method [46]:
∥ρn(t)∥≤∥Sn(0)∥{2ϕ1(2−α)N(α)[1+α(τ−1)]}n, | (3.17) |
∥σn(t)∥≤∥En(0)∥{2ϕ2(2−α)N(α)[1+α(τ−1)]}n, | (3.18) |
∥δn(t)∥≤∥Vn(0)∥{2ϕ3(2−α)N(α)[1+α(τ−1)]}n, | (3.19) |
∥πn(t)∥≤∥Rn(0)∥{2ϕ4(2−α)N(α)[1+α(τ−1)]}n. | (3.20) |
Thus, solutions of the system (2.1)–(2.4) exist, which are given by the functions (or inequalities) (3.17)–(3.20), and are also continuous. Let us now consider that
S(t)−S0=Sn(t)−Sn(t),E(t)−E0=En(t)−En(t),V(t)−V0=Vn(t)−Vn(t),R(t)−R0=Rn(t)−Rn(t). |
Taking Sn(t) into consideration, we now achieve
∥Sn(t)∥=2(2−α)N(α)∥{(1−α)(H1(S,t)−H1(Sn−1,t))}+α∫t0(H1(S,z)−H1(Sn−1,z))dz∥≤ϕ1∥S−Sn−1∥{2(1−α)(2−α)N(α)[1+α(τ−1)]}. |
If this process is iterated then at τ=τ0, we get
∥Sn(t)∥≤{2ϕ1(2−α)N(α)[1+α(τ0−1)]}n+1ω1. |
Now, it is clear that limn→∞∥Sn(t)∥=0. In the same manner it can be put that En(t)→0, Vn(t)→0, and Rn(t)→0 whenever n→∞.
To clarify the uniqueness for the solutions of the system (2.1)–(2.4), we now assume that S′(t), E′(t), V′(t), and R′(t) are a distinct set of solutions for the said system, then
S(t)−S′(t)=2(1−α)(2−α)N(α)(H1(S,t)−H1(S′,t))+2α(2−α)N(α)∫t0(H1(S,z)−H1(S′,z))dz. |
Let us consider that kernel H1 justify the Lipschitz condition and taking the norm on the above equation, we get the inequality given below;
∥S(t)−S′(t)∥{1−2ϕ1(2−α)N(α)[1+α(τ−1)]}≤0. | (3.21) |
Similarly,
∥E(t)−E′(t)∥{1−2ϕ2(2−α)N(α)[1+α(τ−1)]}≤0, | (3.22) |
∥V(t)−V′(t)∥{1−2ϕ3(2−α)N(α)[1+α(τ−1)]}≤0, | (3.23) |
∥R(t)−R′(t)∥{1−2ϕ4(2−α)N(α)[1+α(τ−1)]}≤0. | (3.24) |
Theorem 4. For the system (2.1)–(2.4), a unique solution exists if the following inequalities holds good
{1−2ϕi(2−α)N(α)[1+α(τ−1)]}>0, |
for all i = 1, 2, 3, 4.
Proof. If we consider the assumption made, then from the inequalities given by (3.21)–(3.24), we get
∥S(t)−S′(t)∥=0,∥E(t)−E′(t)∥=0,∥V(t)−V′(t)∥=0,∥R(t)−R′(t)∥=0. |
This implies that S(t)=S′(t), E(t)=E′(t), V(t)=V′(t), and R(t)=R′(t). And this justifies the proof of COVID-19 model (2.1)–(2.4).
The system of equations of COVID-19 ABC model using the fractional integral [38,46], the findings involve the following:
S(t)=ABC0Iαt{λα−γαSE−(dα+να)S}+S0,E(t)=ABC0Iαt{γαSE+(1−η)γαVE−(dα+βα)E}+E0,V(t)=ABC0Iαt{ναS−dαV−(1−η)γαVE}+V0,R(t)=ABC0Iαt{βαE−dαR}+R0. |
The symbols, after arranged in view of the problem, give the following results;
S(t)=S0+(1−α)N(α)[λα−γαS(t)E(t)−(dα+να)S(t)]+αN(α)Γ(α)∫t0(t−y)α−1[λα−γαS(y)E(y)−(dα+να)S(y)]dy,E(t)=E0+(1−α)N(α)[γαS(t)E(t)+(1−η)γαV(t)E(t)−(dα+βα)E(t)]+αN(α)Γ(α)∫t0(t−y)α−1[γαS(y)E(y)+(1−η)γαV(y)E(y)−(dα+βα)E(y)]dy,V(t)=V0+(1−α)N(α)[ναS(t)−dαV(t)−(1−η)γαV(t)E(t)]+αN(α)Γ(α)∫t0(t−y)α−1[ναS(y)−dαV(y)−(1−η)γαV(y)E(y)]dy,R(t)=R0+(1−α)N(α)[βαE(t)−dαR(t)]+αN(α)Γ(α)∫t0[βαE(y)−dαR(y)]dy. |
For the sake of simplicity, we consider
ˉH1(S,t)=[λα−γαS(t)E(t)−(dα+να)S(t)], | (3.25) |
ˉH2(E,t)=[γαS(t)E(t)+(1−η)γαV(t)E(t)−(dα+βα)E(t)], | (3.26) |
ˉH3(V,t)=[ναS(t)−dαV(t)−(1−η)γαV(t)E(t)], | (3.27) |
ˉH4(R,t)=[βαE(t)−dαR(t)]. | (3.28) |
Theorem 5. The Lipschitz condition and contraction are justified by the kernel
i) ˉH1 if 0≤(γαˉω2+dα+να)<1.
ii) ˉH2 if 0≤(γαˉω1+(1−η)γαˉω3+dα+βα)<1.
iii) ˉH3 if 0≤(dα+(1−η)γαˉω2)<1.
iv) ˉH4 if 0≤dα<1.
Proof. Let us assume that the state variables present in the COVID-19 models are bounded, and therefore it can be said that ∥S(t)∥≤ˉω1, ∥E(t)∥≤ˉω2, ∥V(t)∥≤ˉω3, and ∥R(t)∥≤ˉω4 [46].
Let us consider S(t) and S′(t) be the two functions, and therefore
∥ˉH1(S,t)−ˉH1(S′,t)∥=∥−γα(S(t)−S′(t))E(t)−(dα+να)(S(t)−S′(t))∥≤∥γαE(t)(S(t)−S′(t))∥+∥(dα+να)(S(t)−S′(t))∥≤[γα∥E(t)∥+(dα+να)]∥S(t)−S′(t)∥≤[γαˉω2+dα+να]∥S(t)−S′(t)∥. |
So, it can be mentioned now that
∥ˉH1(S,t)−ˉH1(S′,t)∥≤ˉϕ1∥S(t)−S′(t)∥, |
where ˉϕ1=(γαˉω2+dα+να). Thus, the Lipschitz condition is achieved for kernel ˉH1 and 0≤(γαˉω2+dα+να)<1 [46] gives ˉH1, and also contraction is justified.
Like wise, the Lipschitz condition and contraction can be performed and written for the kernels ˉH2, ˉH3, and ˉH4 as mentioned below;
∥ˉH2(E,t)−ˉH2(E′,t)∥≤ˉϕ2∥E(t)−E′(t)∥,∥ˉH3(V,t)−ˉH3(V′,t)∥≤ˉϕ3∥V(t)−V′(t)∥,∥ˉH4(R,t)−ˉH4(R′,t)∥≤ϕ4∥R(t)−R′(t)∥, |
where ˉϕ2=(γαˉω1+(1−η)γαˉω3+dα+βα), ˉϕ3=(dα+(1−η)γαˉω2), and ϕ4=dα.
Taking the kernels ˉH1, ˉH2, ˉH3, and ˉH4 into consideration the following equations are obtained:
S(t)=S0+1N(α)Γ(α)[Γ(α)(1−α)ˉH1(S,t)+α∫t0(t−z)α−1ˉH1(S,z)dz], | (3.29) |
E(t)=E0+1N(α)Γ(α)[Γ(α)(1−α)ˉH2(E,t)+α∫t0(t−z)α−1ˉH2(E,z)dz], | (3.30) |
V(t)=V0+1N(α)Γ(α)[Γ(α)(1−α)ˉH3(V,t)+α∫t0(t−z)α−1ˉH3(V,z)dz], | (3.31) |
R(t)=R0+1N(α)Γ(α)[Γ(α)(1−α)ˉH4(R,t)+α∫t0(t−z)α−1ˉH4(R,z)dz]. | (3.32) |
The following iterative formulae with the initial conditions of the COVID-19 model given in the system of Eqs (2.5)–(2.8) are now seen closely, which are as follows:
S(t)=S0+1N(α)Γ(α)[Γ(α)(1−α)ˉH1(Sn−1,t)+α∫t0(t−z)α−1ˉH1(Sn−1,z)dz], | (3.33) |
E(t)=E0+1N(α)Γ(α)[Γ(α)(1−α)ˉH2(En−1,t)+α∫t0(t−z)α−1ˉH2(En−1,z)dz], | (3.34) |
V(t)=V0+1N(α)Γ(α)[Γ(α)(1−α)ˉH3(Vn−1,t)+α∫t0(t−z)α−1ˉH3(Vn−1,z)dz], | (3.35) |
R(t)=R0+1N(α)Γ(α)[Γ(α)(1−α)ˉH4(Rn−1,t)+α∫t0(t−z)α−1ˉH4(Rn−1,z)dz]. | (3.36) |
The upcoming expressions for the difference of sequential terms are written as follows:
ˉρn(t)=Sn(t)−Sn−1(t)=(1−α)N(1−α)(ˉH1(Sn−1,t)−ˉH1(Sn−2,t))+αN(α)Γ(α)∫t0(t−z)α−1(ˉH1(Sn−1,z)−ˉH1(Sn−2,z))dz,ˉσn(t)=En(t)−En−1(t)=(1−α)N(1−α)(ˉH2(En−1,t)−ˉH2(En−2,t))+αN(α)Γ(α)∫t0(t−z)α−1(ˉH2(En−1,z)−ˉH2(En−2,z))dz, |
ˉδn(t)=Vn(t)−Vn−1(t)=(1−α)N(1−α)(ˉH3(Vn−1,t)−ˉH3(Vn−2,t))+αN(α)Γ(α)∫t0(t−z)α−1(ˉH3(Vn−1,z)−ˉH3(Vn−2,z))dz,ˉπn(t)=Rn(t)−Rn−1(t)=(1−α)N(1−α)(ˉH4(Rn−1,t)−ˉH4(Rn−2,t))+αN(α)Γ(α)∫t0(t−z)α−1(ˉH4(Rn−1,z)−ˉH4(Rn−2,z))dz. |
It is evident that
Sn(t)=n∑m=1ˉρm(t),En(t)=n∑m=1ˉσm(t),Vn(t)=n∑m=1ˉδm(t),Rn(t)=n∑m=1ˉπm(t). |
Involving norm to ˉρn(t) then using triangle inequality, we get
∥ˉρn(t)∥≤(1−α)N(α)∥(ˉH1(Sn−1,t)−ˉH1(Sn−2,t))∥+αN(α)Γ(α)∥∫t0(t−z)α−1(H1(Sn−1,z)−H1(Sn−2,z))dz∥. |
It has been shown that the Lipschitz condition is satisfied by the kernel ˉH1, therefore
∥Sn(t)−Sn−1(t)∥≤(1−α)ˉϕ1N(α)∥Sn−1−Sn−2∥+αˉϕ1N(α)Γ(α)[∫t0(t−z)α−1∥Sn−1−Sn−2∥dz]. |
Thus, it can be deduced that
∥ˉρn(t)∥≤ˉϕ1N(α)Γ(α)[Γ(α)∥ˉρn−1(t)∥+α∫t0(t−z)α−1∥ˉρn−1(z)∥dz]. | (3.37) |
In the same way, other results are followed as:
∥ˉσn(t)∥≤ˉϕ2N(α)Γ(α)[Γ(α)∥ˉσn−1(t)∥+α∫t0(t−z)α−1∥ˉσn−1(z)∥dz], | (3.38) |
∥ˉδn(t)∥≤ˉϕ3N(α)Γ(α)[Γ(α)∥ˉδn−1(t)∥+α∫t0(t−z)α−1∥ˉδn−1(z)∥dz], | (3.39) |
∥ˉπn(t)∥≤ˉϕ4N(α)Γ(α)[Γ(α)∥ˉπn−1(t)∥+α∫t0(t−z)α−1∥ˉπn−1(z)∥dz]. | (3.40) |
Taking the inequalities (3.37)–(3.40) into consideration, the following theorem can be put forward:
Theorem 6. The system (2.5)–(2.8) has a solution if there exist a real number τ0, satisfying
ˉϕiN(α)Γ(α)[Γ(α)(1−α)+τα0]<1, |
for all i=1,2,3,4.
Proof. The Inequations (3.37)–(3.40) are considered along with the fact that the state variables are bounded and the kernels ˉHi (i=1,2,3,4) satisfy the Lipschitz condition. The following mentioned relations employ the iterative method [46]:
∥ˉρn(t)∥≤∥Sn(0)∥{ˉϕ1N(α)Γ(α)[Γ(α)(1−α)+τα0]}n, | (3.41) |
∥ˉσn(t)∥≤∥En(0)∥{ˉϕ2N(α)Γ(α)[Γ(α)(1−α)+τα0]}n, | (3.42) |
∥ˉδn(t)∥≤∥Vn(0)∥{ˉϕ3N(α)Γ(α)[Γ(α)(1−α)+τα0]}n, | (3.43) |
∥ˉπn(t)∥≤∥Rn(0)∥{ˉϕ4N(α)Γ(α)[Γ(α)(1−α)+τα0]}n. | (3.44) |
Thus, solutions of the system (2.5)–(2.8) exist, which are given by the functions (or inequalities) (3.41)–(3.44), and are also continuous. Let us now consider that
S(t)−S0=Sn(t)−ˉSn(t),E(t)−E0=En(t)−ˉEn(t),V(t)−V0=Vn(t)−ˉVn(t),R(t)−R0=Rn(t)−ˉRn(t). |
Taking ˉSn(t) into consideration, we now achieve
∥ˉSn(t)∥=1N(α)Γ(α)∥{(1−α)(ˉH1(S,t)−ˉH1(Sn−1,t))}+α∫t0(t−z)α−1(ˉH1(S,t)−ˉH1(Sn−1,t))dz∥≤ˉϕ1∥Sn−Sn−1∥{1N(α)Γ(α)[Γ(α)(1−α)+τα]}. |
If this process is iterated then at τ=τ0, we get
∥ˉSn(t)∥≤{1N(α)Γ(α)[Γ(α)(1−α)+τα0]}ˉω1. |
Now, it is clear that limn→∞∥ˉSn(t)∥=0. In the same manner it can be put that ˉEn(t)→0, ˉVn(t)→0, and ˉRn(t)→0 whenever n→∞.
To clarify the uniqueness for the solutions of the system (2.5)–(2.8), we now assume that S′(t), E′(t), V′(t), and R′(t) are a distinct set of solutions for the said system, then
S(t)−S′(t)=(1−α)N(α)(ˉH1(S,t)−ˉH1(S′,t))+αN(α)Γ(α)∫t0(t−z)α−1(ˉH1(S,z)−ˉH1(S′,z))dz. |
Let us consider that kernel ˉH1 justify the Lipschitz condition and taking norm on above equation, we get the following inequality given below as:
∥S(t)−S′(t)∥{1−ˉϕ1N(α)Γ(α)[Γ(α)(1−α)+τα]}≤0. | (3.45) |
Similarly,
∥E(t)−E′(t)∥{1−ˉϕ2N(α)Γ(α)[Γ(α)(1−α)+τα]}≤0, | (3.46) |
∥V(t)−V′(t)∥{1−ˉϕ3N(α)Γ(α)[Γ(α)(1−α)+τα]}≤0, | (3.47) |
∥R(t)−R′(t)∥{1−ˉϕ4N(α)Γ(α)[Γ(α)(1−α)+τα]}≤0. | (3.48) |
Theorem 7. For the system (2.5)–(2.8), a unique solution exists if the following inequalities hold good
{1−ˉϕiN(α)Γ(α)[Γ(α)(1−α)+τα]}>0 |
for all i=1,2,3,4.
Proof. If we consider the assumption made, then from the inequalities given by (3.45)–(3.48),
∥S(t)−S′(t)∥=0,∥E(t)−E′(t)∥=0,∥V(t)−V′(t)∥=0,∥R(t)−R′(t)∥=0. |
This implies that S(t)=S′(t), E(t)=E′(t), V(t)=V′(t), and R(t)=R′(t). And this justifies the proof for the uniqueness of the solutions for COVID-19 model (2.5)–(2.8).
The numerical simulation is done using Mathematica for both CF and ABC approach models (2.1)–(2.4) and (2.5)–(2.8), respectively. Next, the comparison and discussion of the following results are done. The dynamics of the novel COVID-19 infection-vaccination model are given by the system (2.1)–(2.4) and (2.5)–(2.8), respectively, for various values of α, randomly selected, are plotted in Figures 2–5. The parametric plot amongst the state variables of the SEVR model are plotted in Figures 6 and 7. The initial values of the state variables are taken as: S0=0.1, E0=0.05, V0=0.5, and R0=0 [41,46].
Due to lack of clinical evidence and support, it is not sure the precise death rate of individuals due to COVID-19 infection, and so we manipulate the value of d to understand the death dynamics. The numerical simulations' comparison of Figures 2(a)–(d), 3(a)–(d), 4(a)–(d), 5(a)–(d) and Figures 2(e)–(h), 3(e)–(h), 4(e)–(h), 5(e)–(h) for both the CF- and ABC approach, respectively, show that for same non-natural values of α, parameter values given in Table 1 and same initial conditions the SEVR model shows different trajectories. It is observed that susceptible populace decrease at a slower rate with the ABC approach as compared to the CF approach for the values of α in Table 1, which can be clearly seen in Figure 2(e), (f) and 2(a), (b). But Figure 2(c), (d) and 2(g), (h) show that both CF- and ABC approach provides precisely same result. It is deduced from Figures 3(a), (b) vs. 3(e), (f) and 4(a), (b) vs. 4(e), (f) about comparisons that the ABC approach gives much variation in effected and vaccinated populaces than the CF approach, for d=0.002 and η=0.5, and 0.8. One more prominent aspect suggests that the effected and vaccinated populaces provide a sharp increase because of high spread of COVID-19 infection presented by Hellewell et al.[47] and Gupta et al. [48]. But Figures 3(c), (d) vs. 3(g), (h) and 4(c), (d) vs. 4(g), (h) explain that both the CF- and ABC approach yields almost same outcome.
ν | α | η | d | β | γ | λ |
0.01 | 0.8, 0.7, 0.6 | 0.5, 0.8 | 0.002, 0.05 | 0.2 | 0.7 | 0.01 |
[41] | estimated | estimated | estimated, [41] | [41] | [41] | [41] |
Moreover, on minute speculation of comparisons shown in Figure 3(a), (b), (e), (f) with Figure 5(a), (b), (e), (f) we see a difference between the peak of effected and recovered populaces, when d=0.002. Be that as it may, it is noted that in Figure 3(c), (d), (g), (h) with Figure 5(c), (d), (g), (h) the difference between the decreasing peak points is approximately 5 days, calculated at d=0.05. The gradient of curves for different values of α share similarities with the findings. This shows that affected individuals are recovering very quickly by the ABC approach, whereas the CF approach shows almost a delay in the transfer of affected individuals to recover compartments.
In all this, we have come across that for almost any efficacy of vaccine, and for arbitrary chosen α if the death rate d=0.002, then affected individuals exhibits virulent characteristics. This may result in tending towards the co-existing equilibrium point X∗. From the Table 2, it can be deduced mathematically that for any vaccine efficacy ranging between 50% to 80% effectiveness can bring stability to the susceptible individuals to fight against the virus. It would be interesting to see how some state variables of the SEVR model are responding to each other for the cases when Rn>1. Therefore, Figures 6 and 7 provide parametric plot between the state variables and it is noticeable that they experience a chaotic behaviour.
η | d | α | Rn | X0 Remark |
0.8 | 5.852 | Unstable | ||
0.002 | 0.7 | 4.434 | Unstable | |
0.6 | 3.341 | Unstable | ||
0.5 | 0.8 | 0.504 | Stable | |
0.05 | 0.7 | 0.495 | Stable | |
0.6 | 0.484 | Stable | ||
0.8 | 3.594 | Unstable | ||
0.002 | 0.7 | 2.821 | Unstable | |
0.6 | 2.206 | Unstable | ||
0.8 | 0.8 | 0.467 | Stable | |
0.05 | 0.7 | 0.454 | Stable | |
0.6 | 0.436 | Stable |
The virtue of the models (2.1)–(2.4) and (2.5)–(2.8), when contrasted, found remarkable differences under the parameter values given in Table 1. The reason is memory properties of the kernel in the fractional derivatives' definitions. The CF and ABC derivative has exponential kernel and generalized Mittag-Leffler kernel, respectively. The ABC type shows a power-law memory, and also partial exponential decay memory [49]. It is clear that the model lies upon the fractional order notably, for various values of α it displays a clear difference and suggested models (2.1)–(2.4) and (2.5)–(2.8) explores new aspects at the fractional values of α.
Infection-related contact ascends with a growth in the size of population. Thus, COVID-19 infection spread during vaccination program for non-integer order model is considered, using the CF and ABC sense. Numerical simulation output for the proposed fractional order models are shown with distinct fractional order values α∈]0,1[ and briefly compared, discussed and investigated. The vaccine efficacy and the death rate for state variables plays a significant role in escaping the infection. Since COVID-19 virus biology is changing rapidly, so the precise estimation of death rate, d, can unfold the appropriate time for vaccine administration and how effective it should be. The graphical results demonstrate the ABC approach provides better suitability for mild cases. Whereas, the CF approach provides superior and more flexible results for critical cases. The CF approach provided a guarantee for use in lenient cases, whereas the ABC approach has a notable impact on the dynamics of COVID-19 infection and has given new insight for its use as a device to further research in dynamics of the spread of the disease in critical cases. These results show that CF and ABC approach implementation in real life situations are both plausible and doable as per the severity of illness for patient management. The findings also underline that even incomplete vaccination can have delay positive repercussions in minding the symptoms and spread of the Coronavirus COVID-19. It is found from the analysis that proposed fractional order mathematical models yield better results than that of the classical ones. The entire manuscript deals with the solution of a dynamical system of COVID-19 model for various fractional order α ranging in the interval [0,1].
For future research work the analysis of COVID-19 infection spread for different geographical areas can be achieved by examining the models with relevant parameter values as per data trends of the region, and also including quarantine and number of dosages of vaccines compartments. It is anticipated that this research will provide significance and will thus strengthen the research relevant to COVID-19 transmission dynamics, so that progressive disease control policies are formulated to provide patients with better medical care to all in need.
The researchers would like to thank the Deanship of Scientific Research, Qassim University for funding the publication of this project. We are also thankful to the editor and anonymous reviewers for their fruitful comments and suggestions that help to improve this manuscript.
The authors declare there is no conflict of interest as well as no known competing financial interests.
[1] |
Nikolic, S., Daniel, S., Haque, R., Belkina, M., Hassan, G.M., Grundy, S., et al., ChatGPT versus Engineering Education Assessment: A Multidisciplinary and Multi-institutional Benchmarking and Analysis of this Generative Artificial Intelligence Tool to Investigate Assessment Integrity. European Journal of Engineering Education, 2023, 48(4): 559‒614. https://doi.org/10.1080/03043797.2023.2213169 doi: 10.1080/03043797.2023.2213169
![]() |
[2] | Mollick, E., Co-intelligence: Living and Working with AI. London. WH Allen. 2024. |
[3] |
Nikolic, S., Sandison, C., Haque, R., Daniel, S., Grundy, S., Belkina, M., et al., ChatGPT, Copilot, Gemini, SciSpace and Wolfram versus Higher Education Assessments: An Updated Multi-Institutional Study of the Academic Integrity Impacts of Generative Artificial Intelligence (GenAI) on Assessment, Teaching and Learning in Engineering. Australasian Journal of Engineering Education, 2024, 29(2): 126‒153. https://doi.org/10.1080/22054952.2024.2372154 doi: 10.1080/22054952.2024.2372154
![]() |
[4] |
Bearman, M., Tai, J., Dawson, P., Boud, D. and Ajjawi, R., Developing evaluative judgement for a time of generative artificial intelligence. Assessment & Evaluation in Higher Education, 2024, 49(6): 893‒905. https://doi.org/10.1080/02602938.2024.2335321 doi: 10.1080/02602938.2024.2335321
![]() |
[5] |
Kizilcec, R.F., Huber, E., Papanastasiou, E.C., Cram, A., Makridis, C.A., Smolansky, A., et al., Perceived impact of generative AI on assessments: Comparing educator and student perspectives in Australia, Cyprus, and the United States. Computers and Education: Artificial Intelligence, 2024, 7: 100269. https://doi.org/10.1016/j.caeai.2024.100269 doi: 10.1016/j.caeai.2024.100269
![]() |
[6] | Quince, Z., Petkoff, K., Michael, R.N., Daniel, S. and Nikolic, S., The current ethical considerations of using GenAI in engineering education and practice: A systematic literature review. 35th Annual Conference of the Australasian Association for Engineering Education, 2024, Christchurch, New Zealand. |
[7] | Hysaj, A., Farouqa, G., Khan, S.A. and Hiasat, L., A Tale of Academic Writing Using AI Tools: Lessons Learned from Multicultural Undergraduate Students. Social Computing and Social Media, Cham. 2024. https://doi.otg/10.1007/978-3-031-61305-0_3 |
[8] | Fatahi, B., Nguyen, L.D., Khabbaz, H. and Hadgraft, R., Virtual Teammates: Transforming Engineering Learning through Generative AI Integration. 35th Australasian Association of Engineering Education Conference, Christchurch, New Zealand. 2024. |
[9] |
Honig, C., Rios, S. and Desu, A., Generative AI in engineering education: understanding acceptance and use of new GPT teaching tools within a UTAUT framework. Australasian Journal of Engineering Education, 2025, 1‒13. https://doi.org/10.1080/22054952.2025.2467500 doi: 10.1080/22054952.2025.2467500
![]() |
[10] |
Ajjawi, R., Tai, J., Dollinger, M., Dawson, P., Boud, D. and Bearman, M., From authentic assessment to authenticity in assessment: broadening perspectives. Assessment & Evaluation in Higher Education, 2023, 49(4): 499‒510. https://doi.org/10.1080/02602938.2023.2271193 doi: 10.1080/02602938.2023.2271193
![]() |
[11] |
Miao, G., Ranaraja, I., Grundy, S., Brown, N., Belkina, M. and Goldfinch, T., Project-based learning in Australian & New Zealand universities: current practice and challenges. Australasian Journal of Engineering Education, 2024, 1‒13. https://doi.org/10.1080/22054952.2024.2358576 doi: 10.1080/22054952.2024.2358576
![]() |
[12] | Mills, J.E. and Treagust, D.F., Engineering education—Is problem-based or project-based learning the answer. Australasian Journal of Engineering Education, 2003, 3(2): 2‒16. |
[13] |
Baig, M.I. and Yadegaridehkordi, E., ChatGPT in the higher education: A systematic literature review and research challenges. International journal of educational research, 2024,127: 102411. https://doi.org/10.1016/j.ijer.2024.102411 doi: 10.1016/j.ijer.2024.102411
![]() |
[14] | Nikolic, S. and Beckman, K., Supporting Engineering Project-Based Learning through the Use of ChatGPT and Generative AI: A Case Study. In H. Crompton & D. Burke (Eds.), Artificial Intelligence Applications in Higher Education: Theories, Ethics, and Case Studies for Universities, 2025,215‒232. Routledge. https://doi.org/10.4324/9781003440178 |
[15] |
Grilli, L. and Pedota, M., Creativity and artificial intelligence: A multilevel perspective. Creativity and Innovation Management, 2024, 33(2): 234‒247. https://doi.org/10.1111/caim.12580 doi: 10.1111/caim.12580
![]() |
[16] | Crawley, E.F., Malmqvist, J., Lucas, W.A. and Brodeur, D.R., The CDIO syllabus v2. 0. An updated statement of goals for engineering education. Proceedings of the 7th International CDIO Conference, 2011. |
[17] |
Lee, D., Arnold, M., Srivastava, A., Plastow, K., Strelan, P., Ploeckl, F., et al., The impact of generative AI on higher education learning and teaching: A study of educators' perspectives. Computers and Education: Artificial Intelligence, 2024, 6: 100221. https://doi.org/10.1016/j.caeai.2024.100221 doi: 10.1016/j.caeai.2024.100221
![]() |
[18] |
Lim, W.M., Gunasekara, A., Pallant, J.L., Pallant, J.I. and Pechenkina, E., Generative AI and the future of education: Ragnarök or reformation? A paradoxical perspective from management educators. The International Journal of Management Education, 2023, 21(2): 100790. https://doi.org/10.1016/j.ijme.2023.100790 doi: 10.1016/j.ijme.2023.100790
![]() |
[19] |
Nikolic, S., Wentworth, I., Sheridan, L., Moss, S., Duursma, E., Jones, R.A., et al., A systematic literature review of attitudes, intentions and behaviours of teaching academics pertaining to AI and generative AI (GenAI) in higher education: An analysis of GenAI adoption using the UTAUT framework. Australasian Journal of Educational Technology, 2024, 40(6): 56‒75. https://doi.org/10.14742/ajet.9643 doi: 10.14742/ajet.9643
![]() |
[20] |
Ahmed, Z., Shanto, S.S. and Jony, A.I., Potentiality of generative AI tools in higher education: Evaluating ChatGPT's viability as a teaching assistant for introductory programming courses. STEM Education, 2024, 4(3): 165‒182. https://doi.org/10.3934/steme.2024011 doi: 10.3934/steme.2024011
![]() |
[21] |
Kim, D., Majdara, A. and Olson, W., A Pilot Study Inquiring into the Impact of ChatGPT on Lab Report Writing in Introductory Engineering Labs. International Journal of Technology in Education, 2024, 7(2): 259‒289. https://doi.org/10.46328/ijte.691 doi: 10.46328/ijte.691
![]() |
[22] | Nikolic, S., Heath, A., Vu, B.A., Daniel, S., Alimardani, A., Sandison, C., et al., Prompt Potential: A Pilot Assessment of Using Generative Artificial Intelligence (ChatGPT-4) as a Tutor for Engineering and Maths. 52nd Annual Conference of the European Society for Engineering Education (SEFI), 2024, Lausanne, Switzerland. |
[23] |
Liu, J., Wang, C., Liu, Z., Gao, M., Xu, Y., Chen, J., et al., A bibliometric analysis of generative AI in education: current status and development. Asia Pacific Journal of Education, 2024, 44(1): 156‒175. https://doi.org/10.1080/02188791.2024.2305170 doi: 10.1080/02188791.2024.2305170
![]() |
[24] |
Bobrytska, V.I., Krasylnykova, H.V., Beseda, N.А., Krasylnykov, S.R. and Skyrda, T.S., Artificial intelligence (AI) in Ukrainian Higher Education: A Comprehensive Study of Stakeholder Attitudes, Expectations and Concerns. International Journal of Learning, Teaching and Educational Research, 2024, 23(1): 400‒426. https://doi.org/10.26803/ijlter.23.1.20 doi: 10.26803/ijlter.23.1.20
![]() |
[25] |
Luo, J., A critical review of GenAI policies in higher education assessment: a call to reconsider the "originality" of students' work. Assessment & Evaluation in Higher Education, 2024, 49(5): 651‒664. https://doi.org/10.1080/02602938.2024.2309963 doi: 10.1080/02602938.2024.2309963
![]() |
[26] |
Southworth, J., Migliaccio, K., Glover, J., Reed, D., McCarty, C., Brendemuhl, J., et al, Developing a model for AI Across the curriculum: Transforming the higher education landscape via innovation in AI literacy. Computers and Education: Artificial Intelligence, 2023, 4: 100127. https://doi.org/10.1016/j.caeai.2023.100127 doi: 10.1016/j.caeai.2023.100127
![]() |
[27] |
Shailendra, S., Kadel, R., and Sharma, A., Framework for Adoption of Generative Artificial Intelligence (GenAI) in Education. IEEE Transactions on Education, 2024, 67(5): 777‒785. https://doi.org/10.1109/TE.2024.3432101 doi: 10.1109/TE.2024.3432101
![]() |
[28] |
Shanto, S.S., Ahmed, Z. and Jony, A.I., PAIGE: A generative AI-based framework for promoting assignment integrity in higher education. STEM Education, 2023, 3(4), 288‒305. https://doi.org/10.3934/steme.2023018 doi: 10.3934/steme.2023018
![]() |
[29] | Ambikairajah, E., Sirojan, T., Thiruvaran, T. and Sethu, V., ChatGPT in the Classroom: A Shift in Engineering Design Education. 2024 IEEE Global Engineering Education Conference (EDUCON), 2024, Kos, Greece. |
[30] |
Salinas-Navarro, D.E., Vilalta-Perdomo, E., Michel-Villarreal, R. and Montesinos, L., Designing experiential learning activities with generative artificial intelligence tools for authentic assessment. Interactive Technology and Smart Education, 2024, 21(4): 708‒734. https://doi.org/10.1108/ITSE-12-2023-0236 doi: 10.1108/ITSE-12-2023-0236
![]() |
[31] |
Prieto, S.A., Mengiste, E.T. and García de Soto, B., Investigating the Use of ChatGPT for the Scheduling of Construction Projects. Buildings, 2023, 13(4): 1‒16. https://doi.org/10.3390/buildings13040857 doi: 10.3390/buildings13040857
![]() |
[32] |
Fosso Wamba, S., Queiroz, M.M., Chiappetta Jabbour, C.J. and Shi, C., Are both generative AI and ChatGPT game changers for 21st-Century operations and supply chain excellence? International journal of production economics, 2023,265: 109015. https://doi.org/10.1016/j.ijpe.2023.109015 doi: 10.1016/j.ijpe.2023.109015
![]() |
[33] | Manresa, A., Sammour, A., Mas-Machuca, M., Chen, W. and Botchie, D., Humanizing GenAI at work: bridging the gap between technological innovation and employee engagement. Journal of Managerial Psychology, 2024, ahead-of-print. https://doi.org/10.1108/JMP-05-2024-0356 |
[34] | Humlum, A. and Vestergaard, E., The Adoption of ChatGPT. University of Chicago, Becker Friedman Institute for Economics Working Paper, 2024(2024-50). |
[35] |
Guo, P., Saab, N., Post, L.S. and Admiraal, W., A review of project-based learning in higher education: Student outcomes and measures. International journal of educational research, 2020,102: 101586. https://doi.org/https://doi.org/10.1016/j.ijer.2020.101586 doi: 10.1016/j.ijer.2020.101586
![]() |
[36] |
Gomez-del Rio, T. and Rodriguez, J., Design and assessment of a project-based learning in a laboratory for integrating knowledge and improving engineering design skills. Education for Chemical Engineers, 2022, 40: 17‒28. https://doi.org/https://doi.org/10.1016/j.ece.2022.04.002 doi: 10.1016/j.ece.2022.04.002
![]() |
[37] | Gregory, S., O'Connell, J., Butler, D., McDonald, M., Kerr, T., Schutt, S., et al., New applications, new global audiences: Educators repurposing and reusing 3D virtual and immersive learning resources. 2015 ASCILITE Annual Conference, 2015, Perth, Australia. |
[38] |
Nikolic, S., Suesse, T.F., Grundy, S., Haque, R., Lyden, S., Lal, S., et al., Assessment integrity and validity in the teaching laboratory: adapting to GenAI by developing an understanding of the verifiable learning objectives behind laboratory assessment selection. European Journal of Engineering Education, 2025, 1‒28. https://doi.org/10.1080/03043797.2025.2456944 doi: 10.1080/03043797.2025.2456944
![]() |
[39] |
Lee, M.J.W., Nikolic, S., Vial, P.J., Ritz, C., Li, W. and Goldfinch, T., Enhancing project-based learning through student and industry engagement in a video-augmented 3-D virtual trade fair. IEEE Transactions on Education, 2016, 59(4): 290‒298. https://doi.org/10.1109/TE.2016.2546230 doi: 10.1109/TE.2016.2546230
![]() |
[40] |
Witarsa, and Muhammad, S., Critical thinking as a necessity for social science students capacity development: How it can be strengthened through project based learning at university. Frontiers in Education, 2023, 7(1). https://doi.org/10.3389/feduc.2022.983292 doi: 10.3389/feduc.2022.983292
![]() |
[41] |
Markula, A. and Aksela, M., The key characteristics of project-based learning: how teachers implement projects in K-12 science education. Disciplinary and Interdisciplinary Science Education Research, 2022, 4(1): 2. https://doi.org/10.1186/s43031-021-00042-x doi: 10.1186/s43031-021-00042-x
![]() |
[42] |
Beneroso, D. and Robinson, J., Online project-based learning in engineering design: Supporting the acquisition of design skills. Education for Chemical Engineers, 2022, 38: 38‒47. https://doi.org/https://doi.org/10.1016/j.ece.2021.09.002 doi: 10.1016/j.ece.2021.09.002
![]() |
[43] |
Edström, K. and Kolmos, A., PBL and CDIO: complementary models for engineering education development. European Journal of Engineering Education, 2014, 39(5): 539‒555. https://doi.org/10.1080/03043797.2014.895703 doi: 10.1080/03043797.2014.895703
![]() |
[44] | Crawley, E.F., The CDIO Syllabus: A statement of goals for undergraduate engineering education. Massachusetts Institute of Technology Cambridge. 2001. |
[45] | Ramírez de Dampierre, M., Gaya-López, M.C. and Lara-Bercial, P.J., Evaluation of the Implementation of Project-Based-Learning in Engineering Programs: A Review of the Literature. Education Sciences, 2024, 14(10): 1107. https://www.mdpi.com/2227-7102/14/10/1107 |
[46] |
Tanveer, B. and Usman, M., An Empirical Study on the Use of CDIO in Software Engineering Education. IEEE Transactions on Education, 2022, 65(4): 684‒694. https://doi.org/10.1109/TE.2022.3163911 doi: 10.1109/TE.2022.3163911
![]() |
[47] |
O'Connor, S., Power, J. and Blom, N., A systematic review of CDIO knowledge library publications (2010–2020): An Overview of trends and recommendations for future research. Australasian Journal of Engineering Education, 2023, 28(2): 166‒180. https://doi.org/10.1080/22054952.2023.2220265 doi: 10.1080/22054952.2023.2220265
![]() |
[48] | Malmqvist, J., Lundqvist, U., Rosén, A. and Edström, K., The CDIO syllabus 3.0-an updated statement of goals. 18th International CDIO Conference, June 13-15 2022, Reykjavik. 2022. |
[49] |
Sablatzky, T., The Delphi method. Hypothesis: Research Journal for Health Information Professionals, 2022, 34(1). https://doi.org/10.18060/26224 doi: 10.18060/26224
![]() |
[50] | Linstone, H.A. and Turoff, M., The delphi method, Addison-Wesley Reading, MA, 1975. |
[51] |
Kokotsaki, D., Menzies, V. and Wiggins, A., Project-based learning: A review of the literature. Improving Schools, 2016, 19(3): 267‒277. https://doi.org/10.1177/1365480216659733 doi: 10.1177/1365480216659733
![]() |
[52] |
Francisco, M.G., Canciglieri Junior, O. and Sant'Anna, Â.M.O., Design for six sigma integrated product development reference model through systematic review. International Journal of Lean Six Sigma, 2020, 11(4): 767‒795. https://doi.org/10.1108/IJLSS-05-2019-0052 doi: 10.1108/IJLSS-05-2019-0052
![]() |
[53] | Dowling, D., Hadgraft, R., Carew, A., McCarthy, T., Hargreaves, D., Baillie, C., et al., Engineering Your Future, Melbourne. Wiley. 2019. |
[54] |
Bukar, U.A., Sayeed, M.S., Razak, S.F.A., Yogarayan, S., and Sneesl, R., Decision-Making Framework for the Utilization of Generative Artificial Intelligence in Education: A Case Study of ChatGPT. IEEE Access, 2024. https://doi.org/10.1109/ACCESS.2024.3425172 doi: 10.1109/ACCESS.2024.3425172
![]() |
[55] |
Guo, W., Li, W. and Tisdell, C.C., Effective pedagogy of guiding undergraduate engineering students solving first-order ordinary differential equations. Mathematics, 2021, 9(14): 1623. https://doi.org/10.3390/math9141623. doi: 10.3390/math9141623
![]() |
1. | Shuangfei Li, Yingxian Zhu, Yunxian Dai, Yiping Lin, Stability Switching Curves and Hopf Bifurcation of a Fractional Predator–Prey System with Two Nonidentical Delays, 2022, 14, 2073-8994, 643, 10.3390/sym14040643 | |
2. | Houssine Zine, El Mehdi Lotfi, Delfim F. M. Torres, Noura Yousfi, Taylor’s Formula for Generalized Weighted Fractional Derivatives with Nonsingular Kernels, 2022, 11, 2075-1680, 231, 10.3390/axioms11050231 | |
3. | Moirangthem Bidyaluxmi Devi, Arpita Devi, Praveen Kumar Gupta, Dharmendra Tripathi, Response of vaccination on community transmission of COVID-19: a dynamical approach, 2022, 231, 1951-6355, 3749, 10.1140/epjs/s11734-022-00652-0 | |
4. | H Yépez-Martínez, J F Gómez-Aguilar, Mustafa Inc, New modied Atangana-Baleanu fractional derivative applied to solve nonlinear fractional dierential equations, 2023, 98, 0031-8949, 035202, 10.1088/1402-4896/acb591 | |
5. | A. K. Misra, Jyoti Maurya, Mohammad Sajid, Modeling the effect of time delay in the increment of number of hospital beds to control an infectious disease, 2022, 19, 1551-0018, 11628, 10.3934/mbe.2022541 | |
6. | Djalal Boucenna, Amar Chidouh, Delfim F. M. Torres, Existence Results for a Multipoint Fractional Boundary Value Problem in the Fractional Derivative Banach Space, 2022, 11, 2075-1680, 295, 10.3390/axioms11060295 | |
7. | Mohammed Subhi Hadi, Bülent Bilgehan, Fractional COVID-19 Modeling and Analysis on Successive Optimal Control Policies, 2022, 6, 2504-3110, 533, 10.3390/fractalfract6100533 | |
8. | Mohammad Sajid, Biplab Dhar, Ahmed S. Almohaimeed, Differential order analysis and sensitivity analysis of a CoVID-19 infection system with memory effect, 2022, 7, 2473-6988, 20594, 10.3934/math.20221129 | |
9. | Asish Adak, Arpita Devi, Praveen Kumar Gupta, 2024, 3045, 0094-243X, 030009, 10.1063/5.0201315 | |
10. | Masuda Akter, Mohammad Sharif Ullah, Mutum Zico Meetei, Abdullah A. Zaagan, Ali M. Mahnashi, An innovative fractional-order evolutionary game theoretical study of personal protection, quarantine, and isolation policies for combating epidemic diseases, 2024, 14, 2045-2322, 10.1038/s41598-024-61211-2 | |
11. | Masakazu Onitsuka, Iz-iddine EL-Fassi, GENERALIZED CAPUTO-FABRIZIO FRACTIONAL DIFFERENTIAL EQUATION, 2024, 14, 2156-907X, 964, 10.11948/20230221 | |
12. | Ritu Agarwal, Pooja Airan, Mohammad Sajid, Numerical and graphical simulation of the non-linear fractional dynamical system of bone mineralization, 2024, 21, 1551-0018, 5138, 10.3934/mbe.2024227 | |
13. | P K Santra, G S Mahapatra, Sanjoy Basu, Stability analysis of fractional epidemic model for two infected classes incorporating hospitalization impact, 2024, 99, 0031-8949, 065237, 10.1088/1402-4896/ad4692 | |
14. | Biplab Dhar, Manisha Malik, Mathematical analysis of a new dimensionalized fractal-fractional operator with application to colorectal cancer treatment model, 2025, 1598-5865, 10.1007/s12190-025-02464-6 |
η | d | α | Rn | X0 Remark |
0.8 | 5.852 | Unstable | ||
0.002 | 0.7 | 4.434 | Unstable | |
0.6 | 3.341 | Unstable | ||
0.5 | 0.8 | 0.504 | Stable | |
0.05 | 0.7 | 0.495 | Stable | |
0.6 | 0.484 | Stable | ||
0.8 | 3.594 | Unstable | ||
0.002 | 0.7 | 2.821 | Unstable | |
0.6 | 2.206 | Unstable | ||
0.8 | 0.8 | 0.467 | Stable | |
0.05 | 0.7 | 0.454 | Stable | |
0.6 | 0.436 | Stable |
ν | α | η | d | β | γ | λ |
0.01 | 0.8, 0.7, 0.6 | 0.5, 0.8 | 0.002, 0.05 | 0.2 | 0.7 | 0.01 |
[41] | estimated | estimated | estimated, [41] | [41] | [41] | [41] |
η | d | α | Rn | X0 Remark |
0.8 | 5.852 | Unstable | ||
0.002 | 0.7 | 4.434 | Unstable | |
0.6 | 3.341 | Unstable | ||
0.5 | 0.8 | 0.504 | Stable | |
0.05 | 0.7 | 0.495 | Stable | |
0.6 | 0.484 | Stable | ||
0.8 | 3.594 | Unstable | ||
0.002 | 0.7 | 2.821 | Unstable | |
0.6 | 2.206 | Unstable | ||
0.8 | 0.8 | 0.467 | Stable | |
0.05 | 0.7 | 0.454 | Stable | |
0.6 | 0.436 | Stable |