Collective creativity in science classrooms: Scale adaptation and an investigation in terms of demographic variables

Menşure Alkış Küçükaydın; Çiğdem Akkanat Avşar; Menşure Alkış Küçükaydın; Çiğdem Akkanat Avşar

doi:10.3934/steme.2025011

STEM Education

2025, Volume 5, Issue 2: 207-228. doi: 10.3934/steme.2025011

Previous Article Next Article

Research article Topical Sections

Collective creativity in science classrooms: Scale adaptation and an investigation in terms of demographic variables

Menşure Alkış Küçükaydın ^{1,
,
,},
Çiğdem Akkanat Avşar ²

1.
Department of Basic Education, Necmettin Erbakan University, Eregli Faculty of Education, 42310 Ereğli, Konya, Türkiye
2.
Ministry of National Education, Ankara, Türkiye; cakkanat@gmail.com

Academic Editor: Giuseppe Carbone

Received: 08 October 2024 Revised: 21 January 2025 Accepted: 22 January 2025 Published: 21 February 2025

Academics and traditional schooling normally focus on achievement in the classroom, and so students' inherent creativity rarely gets explored or expressed. In this study, the psychometric properties of the Turkish version of the Science Classroom Creativity Scale (SCC Scale) were evaluated. The scale has nine factors and a total of 49 items to determine science-specific creativity. A total of 422 students in grades 5–10 in the Central Anatolia Region of Turkey participated in the study. The results showed that the model fit values of the SCC scale are as follows: χ2 /df = 2.07, RMSEA = 0.05 [90% CI: 0.065; 0.078], S-RMR = 0.05, AGFI = 0.91, NFI = 0.90, IFI = 0.95, GFI = 0.91, CFI = 0.92, TLI = 0.91. Also reliability analyses showed that Cronbach's alpha coefficient ranged between 0.90 and 0.96, while McDonald's omega coefficient ranged between 0.90 and 0.98. The validity evidence also indicated that the scale was appropriate for the Turkish sample. In addition, students' creativity was analyzed in terms of demographic variables. Accordingly, the science-specific creativity level varied according to grade level, but there was no statistical difference according to gender. The results were discussed in terms of student creativity in science classrooms.

Keywords:

Citation: Menşure Alkış Küçükaydın, Çiğdem Akkanat Avşar. Collective creativity in science classrooms: Scale adaptation and an investigation in terms of demographic variables[J]. STEM Education, 2025, 5(2): 207-228. doi: 10.3934/steme.2025011

Related Papers:

[1]	Hasib Khan, Jehad Alzabut, Anwar Shah, Sina Etemad, Shahram Rezapour, Choonkil Park . A study on the fractal-fractional tobacco smoking model. AIMS Mathematics, 2022, 7(8): 13887-13909. doi: 10.3934/math.2022767
[2]	Weerawat Sudsutad, Jutarat Kongson, Chatthai Thaiprayoon, Nantapat Jarasthitikulchai, Marisa Kaewsuwan . A generalized Gronwall inequality via $ \psi $-Hilfer proportional fractional operators and its applications to nonlocal Cauchy-type system. AIMS Mathematics, 2024, 9(9): 24443-24479. doi: 10.3934/math.20241191
[3]	Mdi Begum Jeelani, Abeer S. Alnahdi, Mohammed A. Almalahi, Mohammed S. Abdo, Hanan A. Wahash, M. A. Abdelkawy . Study of the Atangana-Baleanu-Caputo type fractional system with a generalized Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(2): 2001-2018. doi: 10.3934/math.2022115
[4]	Khaled M. Saad, Manal Alqhtani . Numerical simulation of the fractal-fractional reaction diffusion equations with general nonlinear. AIMS Mathematics, 2021, 6(4): 3788-3804. doi: 10.3934/math.2021225
[5]	Rahat Zarin, Amir Khan, Pushpendra Kumar, Usa Wannasingha Humphries . Fractional-order dynamics of Chagas-HIV epidemic model with different fractional operators. AIMS Mathematics, 2022, 7(10): 18897-18924. doi: 10.3934/math.20221041
[6]	Aziz Khan, Thabet Abdeljawad, Manar A. Alqudah . Neural networking study of worms in a wireless sensor model in the sense of fractal fractional. AIMS Mathematics, 2023, 8(11): 26406-26424. doi: 10.3934/math.20231348
[7]	Muhammad Farman, Ali Akgül, Sameh Askar, Thongchai Botmart, Aqeel Ahmad, Hijaz Ahmad . Modeling and analysis of fractional order Zika model. AIMS Mathematics, 2022, 7(3): 3912-3938. doi: 10.3934/math.2022216
[8]	Yuehong Zhang, Zhiying Li, Wangdong Jiang, Wei Liu . The stability of anti-periodic solutions for fractional-order inertial BAM neural networks with time-delays. AIMS Mathematics, 2023, 8(3): 6176-6190. doi: 10.3934/math.2023312
[9]	Amir Ali, Abid Ullah Khan, Obaid Algahtani, Sayed Saifullah . Semi-analytical and numerical computation of fractal-fractional sine-Gordon equation with non-singular kernels. AIMS Mathematics, 2022, 7(8): 14975-14990. doi: 10.3934/math.2022820
[10]	Mohamed A. Barakat, Abd-Allah Hyder, Doaa Rizk . New fractional results for Langevin equations through extensive fractional operators. AIMS Mathematics, 2023, 8(3): 6119-6135. doi: 10.3934/math.2023309

Abstract

1. Introduction

The biological sciences have advanced at a breakneck pace during the last few decades, and it's reasonable to predict that this trend will continue, aided by massive technological advancements in a long time. Mathematics has always both contributed to and benefited society. There is huge progress in the natural sciences because of mathematics, and it can do the same for biological science ^[1]. Biology generates intriguing problems, and mathematics develops models to help us understand them, and biology then tests the mathematical models. The handling of complex mathematical systems has been made easier thanks to recent improvements in computer algebra systems. This has allowed scientists to concentrate on understanding mathematical biology rather than the mechanics of solving problems ^[2].

Cancer has always been considered a complex system as it causes different diseases around two hundred with different characteristics. Therefore, many scientists are still trying to investigate the interactions between tumor cells and immune cells by employing different strategies to better understand the dynamics of cancer ^[3,4]. Research is devoted to understanding the interaction between the immune system and tumor cells ^[5]. Ordinary differential equations mathematical models have been proven valuable to understand the dynamics of the tumor-immune system how host immune cells and cancerous cells evolve interact; See e.g. ^[6,7,8]. However, fractional-order differential equations have more characteristics than classical derivatives in mathematical modeling.

Accordingly, the subject of fractional calculus has gained popularity and importance due to its demonstrated applications in system biology ^[9], and different fields of sciences ^[10]. Derivatives and integrals of any non-negative order are allowed in fractional calculus. The benefit of fractional derivatives (integral as well) is that they are not a local (or point) attribute (or quantity) ^[11]. We noticed that in epidemiology, mathematical models are frequently employed to comprehend the complexity of infectious diseases. The stability theory of differential equations is used to study the dysentery modeling approach with controls ^[12]. The most often employed operators, such as Caputo-Fabrizio, and Atangana-Baleanu, involve a local derivative with an exponential function, power-law, and a Mittag-Leffler function, respectively. On the other hand, real-world processes display intricate fractal behavior, including volatility, aquifers, porous media, biomedicine, and Darcy's law. In this context, it is necessary to have operators involving nonlocal differentiation in the kernel, which Atangana provided in the form of fractional fractal operators ^[13]. In the field of mathematical biology, epidemiological models design has been examined using classical and certain fractional-order operators such as conformable derivatives, Caputo-Fabrizio, Atangana-Baleanu, beta derivatives, and a few more. Dengue fever, tuberculosis, measles, Ebola, and other diseases are among the epidemiological designs supported by fractional operators. Nevertheless, fresh real-life applications studied on fractal-fractional operators have shown that the operator's effectiveness is best suited to curves derived under Caputo type fractal operators using actual banking data ^[14]. Fractal-fractional derivative and integration analysis develop complicated behaviors of diarrhoeal sickness. These operators have created the best curve of the actual data of the diseased population, and the model of diarrhea proves that these novel operators have a unique solution ^[15]. Some application of fractional order model with the local and nonlocal, nonsingular kernel is also studied in ^{[16,17,18,19,20,21]} having memory results for a dynamical system.

Fractional derivative in the Caputo sense is studied in ^[29,30,31]. A vitro model of HER 2 + breast cancer cells dynamics resulting from many dosages and timing of paclitaxel and trastuzumab combination regimens was proposed in ^[30]. Study the pattern and the trend of spread of this disease and prescribe a mathematical model which governs COVID-19 pandemic using Caputo type derivative ^[32]. Fractal-fractional operator in the Caputo sense, the condition for Ulam's type of stability of the solution to the models and its related work is in ^{[33,34,35,36]}. Distributed order time-fractional constitution model is put forward to analyze the unsteady natural convection boundary layer flow and heat transfer, in which the magnetic field effect is considered and its related work is in ^[37,38,39].

This paper, in Sections 1 and 2, consists of an introduction and basic definitions for analysis. Section 3 is for the Ulam-Hyres stability and uniqueness of the proposed scheme of the fractional order cancer model with the fractal fractional operator by using generalized Mittag-Leffler Kernel. A numerical algorithm for simulation and results is developed in Section 4. Description and conclusion of results are described in Sections 5 and 6.

2. Preliminaries

We consider some basic refinement of fractional calculus ^{[22,23,24,25]} in this section which are helpful for analysis and simulation of the problem.

Definition 2.1. Let $0\leq \xi$ , $\eta \leq 1$ , then $U(t)$ in the Riemann-Liouville for fractal-fractional with power-law kernel and the fractal-fractional integral have been given as ^[26]:

$\begin{eqnarray} {}_{}^{FFP} D_{0, t}^{\xi, \eta} (U(t)) = \frac{1}{\Gamma(m- \xi)} \frac{d}{dt^{\xi}} \int^{t}_{0}(t-\Psi )^{m - \xi -1} U(\Psi) d \Psi, \end{eqnarray}$

$\begin{eqnarray} \frac{d}{d\Psi^{\eta}} U(\Psi) = lim_{t\rightarrow \Psi} \frac{U(t) - U(\Psi)}{t^{\eta} - \Psi^{\eta}} \end{eqnarray}$

and

$\begin{eqnarray} {}_{}^{FFP} I_{0, t}^{\xi, \eta} (U(t)) = \frac{1}{\Gamma(\xi)} \int^{t}_{0} (t-\Psi )^{\upsilon -1} {\Psi}^{1-\eta} U(\Psi) d \Psi. \end{eqnarray}$

Definition 2.2. Let $0\leq \xi$ , $\eta \leq 1.$ Then $U(t)$ in the Riemann-Liouville for fractal fractional operator having exponentially decaying kernel and fractal-fractional integral have been presented by ^[26,27]:

$\begin{eqnarray} {}_{}^{FFE} D_{0, t}^{\xi, \eta} (U(t)) = \frac{M(\xi)}{\Gamma(m- \xi)} \frac{d}{dt^{\eta}} \int^{t}_{0} exp [- \frac{\xi}{1- \xi} (t-\Psi )^{n - \xi -1} U(\Psi) d \Psi, \end{eqnarray}$

and

$\begin{eqnarray} {}_{}^{FFE} D_{0, t}^{\xi, \eta} (U(t)) = \frac{\eta(1-\xi) t^{\eta -1} U(t)}{M(\xi)} + \frac{\xi \eta}{M(\xi)} \int^{t}_{0} {\Psi}^{\xi -1} U(\Psi) d \Psi. \end{eqnarray}$

Definition 2.3. Let $0\leq \xi$ , $\eta \leq 1$ then $U(t).$ Then the Riemann-Liouville for fractal fractional operator with generalized Mittag-Leffler kernel and fractal-fractional integral have been presented by ^[26,27]:

$\begin{eqnarray} {}_{}^{FFM} D_{0, t}^{\xi, \eta} (U(t)) = \frac{AB(\xi)}{1- \xi} \frac{d}{dt^{\eta}} \int^{t}_{0} E_{\xi} [- \frac{\xi}{1- \xi} (t-\Psi )^{\xi} U(\Psi) d \Psi, \end{eqnarray}$

and

$\begin{eqnarray} {}_{}^{FFM} D_{0, t}^{\xi, \eta} (U(t)) = \frac{\eta(1-\xi) t^{\eta -1} U(t)}{AB(\xi)} + \frac{\xi \eta}{AB(\xi)} \int^{t}_{0} {\Psi}^{\xi -1} (t-\Psi) U(\Psi) d \Psi. \end{eqnarray}$

3. Cancer model with fractal-fractional operator

Recently, many researchers studies to examine the mathematical model of cancer-immune ^[27,28]. By this motivation, IL-12 cytokine is used to increase it ability by increasing the number of CD4+T, CD8+T lymphocytes, and the help of anti-PD-L1 inhibitor with effect of CD8+T cells are added in Castiglione-Piccoli model. So, ^[27,28] explains the new time-dependent ordinary differential system. Equation (1) tells about components of immune system against cancer cells, which is made by adding IL-12 and anti-PD-L1 to components of immune system. The main objective of adding new variables in (1) is to fight cancer cells for the immune system more effectively. The system of nonlinear fractional order with Fractal Fractional operator is given as

$\begin{eqnarray} \begin{cases} ^{FFM}D^{\xi, \eta}_{0, t}(H(t)) = a_{0}+b_{0}DH(1-\frac{H}{f_{0}})+\lambda_{42}\frac{I_{2}}{K_{2}+I_{2}}H(1-\frac{H}{f_{0}})+\lambda_{412}\frac{I_{12}}{K_{12}+I_{12}}H(1-\frac{H}{f_{0}})-c_{0}H, \\ ^{FFM}D^{\xi, \eta}_{0, t}(M(t)) = a_{1}+b_{1}\frac{I_{2}}{K_{2}+I_{2}}(M+D)C(1-\frac{C}{f_{1}})+\lambda_{812}\frac{I_{12}}{K_{12}+I_{12}}C(1-\frac{C}{f_{1}})-c_{1}C, \\ ^{FFM}D^{\xi, \eta}_{0, t}(C(t)) = b_{2}M(1-\frac{M}{f_{2}})-d_{2}FMC, \\ ^{FFM}D^{\xi, \eta}_{0, t}(D(t)) = -d_{3}DC, \\ ^{FFM}D^{\xi, \eta}_{0, t}(I_{2}(t)) = b_{4}DH-e_{4}I_{2}C-c_{4}I_{2}, \\ ^{FFM}D^{\xi, \eta}_{0, t}(I_{12}(t)) = \lambda_{D_{I_{12}}}D-d_{I_{12}}I_{12}, \\ ^{FFM}D^{\xi, \eta}_{0, t}(Z(t)) = -\gamma {Z}. \end{cases} \end{eqnarray}$

(3.1)

where

$F = \frac{c_{pd}-1}{\pi}(\tan^{-1}((Z-1)k_{pd})+\frac{\pi}{2})$

$H_0 = H(0), C_0 = C(0), M_0 = M(0),$

$D_0 = D(0), I_{2_{0}} = I_{2}(0), I_{12_{0}} = I_{12}(0), Z_0 = Z(0)$

Here $H, C, M, D, I L - 2, I L-12$ and $Z$ symbolize $CD_4^+T$ and $CD_8^+T$ lymphocytes, cancer cells, dendritic cells, IL-2 and IL-12 cytokine, anti-PD-L1, respectively given in ^[27,28]

Theorem 1. The solution of the given cancer fractal-fractional model (3.1) along initial conditions is unique and bounded in $R_+^7$ .

Proof. We obtain

$\begin{eqnarray} ^{FFM}D^{\xi, \eta}_{0, t}(H(t))_{H = 0} & = & a_{0} \geq 0, \\ ^{FFM}D^{\xi, \eta}_{0, t}(M(t))_{M = 0} & = & a_{1}+\Big[b_{1}\frac{I_{2}}{K_{2}+I_{2}}D +\lambda_{812}\frac{I_{12}}{K_{12}+I_{12}}\Big]C(1-\frac{C}{f_{1}})-c_{1}C, \geq 0, \\ ^{FFM}D^{\xi, \eta}_{0, t}(C(t))_{C = 0} & = & b_{2}M(1-\frac{M}{f_{2}}) \geq 0, \\ ^{FFM}D^{\xi, \eta}_{0, t}(D(t))|_{D = 0} & = & 0 \geq 0, \\ ^{FFM}D^{\xi, \eta}_{0, t}(I_{2}(t))|_{I_{2} = 0} & = & b_{4}DH \geq 0, \\ ^{FFM}D^{\xi, \eta}_{0, t}(I_{12}(t))|_{I_{12} = 0} & = & \lambda_{D_{I_{12}}}D \geq 0, \\ ^{FFM}D^{\xi, \eta}_{0, t}(Z(t))|_{Z = 0} & = & 0 \geq 0. \end{eqnarray}$

(3.2)

If $(S(0), U(0), C(0), C_{a}(0), R(0), B(0) \in R_+^6$ , then according to Eq (3.2) thee solution cannot escape from the hyperplane. Also on each hyperplane bounding the non-negative orthant, the vector field points into $R_+^7$ , i.e., the domain $R_+^7$ is a positively invariant set.

Stability and existences of proposed operator

Here, we consider stability and existences with proposed operator,

$\begin{eqnarray} ^{FFM}D^{\xi, \eta}_{0, t}(H(t))& = &a_{0}+b_{0}DH(1-\frac{H}{f_{0}})+\lambda_{42}\frac{I_{2}}{K_{2}+I_{2}}H(1-\frac{H}{f_{0}})+\lambda_{412}\frac{I_{12}}{K_{12}+I_{12}}H(1-\frac{H}{f_{0}})-c_{0}H, \\ ^{FFM}D^{\xi, \eta}_{0, t}(M(t))& = &a_{1}+b_{1}\frac{I_{2}}{K_{2}+I_{2}}(M+D)C(1-\frac{C}{f_{1}})+\lambda_{812}\frac{I_{12}}{K_{12}+I_{12}}C(1-\frac{C}{f_{1}})-c_{1}C, \\ ^{FFM}D^{\xi, \eta}_{0, t}(C(t))& = &b_{2}M(1-\frac{M}{f_{2}})-d_{2}FMC, \\ ^{FFM}D^{\xi, \eta}_{0, t}(D(t))& = &-d_{3}DC, \\ ^{FFM}D^{\xi, \eta}_{0, t}(I_{2}(t))& = &b_{4}DH-e_{4}I_{2}C-c_{4}I_{2}, \\ ^{FFM}D^{\xi, \eta}_{0, t}(I_{12}(t))& = &\lambda_{D_{I_{12}}}D-d_{I_{12}}I_{12}, \\ ^{FFM}D^{\xi, \eta}_{0, t}(Z(t))& = &-\gamma {Z}. \end{eqnarray}$

(3.3)

For existence, we have

$\begin{eqnarray} \begin{cases} ^{ABR}_{0}D^{\xi, \eta}_{0, t}(H(t)) = \eta t^{\eta-1}Z_{1}(t, H, M, C, D, I_{2}, I_{12}, Z), \\ ^{ABR}_{0}D^{\xi, \eta}_{0, t}(M(t)) = \eta t^{\eta-1}Y_{1}(t, H, M, C, D, I_{2}, I_{12}, Z), \\ ^{ABR}_{0}D^{\xi, \eta}_{0, t}(C(t)) = \eta t^{\eta-1}X_{1}(t, H, M, C, D, I_{2}, I_{12}, Z), \\ ^{ABR}_{0}D^{\xi, \eta}_{0, t}(D(t)) = \eta t^{\eta-1}W_{1}(t, H, M, C, D, I_{2}, I_{12}, Z), \\ ^{ABR}_{0}D^{\xi, \eta}_{0, t}(I_{2}(t)) = \eta t^{\eta-1}V_{1}(t, H, M, C, D, I_{2}, I_{12}, Z), \\ ^{ABR}_{0}D^{\xi, \eta}_{0, t}(I_{12}(t)) = \eta t^{\eta-1}U_{1}(t, H, M, C, D, I_{2}, I_{12}, Z), \\ ^{ABR}_{0}D^{\xi, \eta}_{0, t}(Z(t)) = \eta t^{\eta-1}E_{1}(t, H, M, C, D, I_{2}, I_{12}, Z). \end{cases} \end{eqnarray}$

(3.4)

where

$\begin{eqnarray} \begin{cases} Z(t, H, M, C, D, I_{2}, I_{12}, Z) = a_{0}+b_{0}DH(1-\frac{H}{f_{0}})+\lambda_{42}\frac{I_{2}}{K_{2}+I_{2}}H(1-\frac{H}{f_{0}})+\lambda_{412}\frac{I_{12}}{K_{12}+I_{12}}H(1-\frac{H}{f_{0}})-c_{0}H, \\ Y(t, H, M, C, D, I_{2}, I_{12}, Z) = a_{1}+b_{1}\frac{I_{2}}{K_{2}+I_{2}}(M+D)C(1-\frac{C}{f_{1}})+\lambda_{812}\frac{I_{12}}{K_{12}+I_{12}}C(1-\frac{C}{f_{1}})-c_{1}C, \\ X(t, H, M, C, D, I_{2}, I_{12}, Z) = b_{2}M(1-\frac{M}{f_{2}})-d_{2}FMC, \\ W(t, H, M, C, D, I_{2}, I_{12}, Z) = -d_{3}DC, \\ V(t, H, M, C, D, I_{2}, I_{12}, Z) = b_{4}DH-e_{4}I_{2}C-c_{4}I_{2}, \\ U(t, H, M, C, D, I_{2}, I_{12}, Z) = \lambda_{D_{I_{12}}}D-d_{I_{12}}I_{12}, \\ E(t, H, M, C, D, I_{2}, I_{12}, Z) = -\gamma Z. \nonumber \end{cases} \end{eqnarray}$

We can write system (3.4) as:

$\left\{ \begin{array}{ll} ^{ABR}_{0}D^{\xi}_{t}\Omega(t) = \eta t^{\eta-1}\Lambda(t, \Omega(t)), \\ \Omega(0) = \Omega_{0}. \end{array} \right.$

By replacing $^{ABR}_{0}D^{\xi, \eta}_{0}$ by $^{ABC}_{0}D^{\xi, \eta}_{0}$ and applying fractional integral, we get

$\Omega(t) = \Omega(0)+\frac{\eta t^{\eta-1}(1-\xi)}{CD(\xi)}\Lambda(t, \Omega(t))+\frac{\xi\eta}{CD(\xi)\Gamma(\xi)}\int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(t, \Omega(t))d\omega,$

where

$\begin{eqnarray} \Omega(t) = \begin{cases} H(t)\\ M(t)\\ C(t)\\ D(t) \\ I_{2}(t)\\ I_{12}(t)\\ Z(t) \nonumber\\ \end{cases} \Omega(0) = \begin{cases} H(0)\\ M(0)\\ C(0)\\ D(0)\\ I_{2}(0)\\ I_{12}(0)\\ Z(0) \nonumber\\ \end{cases} \Lambda(t, \Omega(t)) = \begin{cases} Z(t, H, M, C, D, I_{2}, I_{12}, Z), \\ Y(t, H, M, C, D, I_{2}, I_{12}, Z), \\ X(t, H, M, C, D, I_{2}, I_{12}, Z), \\ W(t, H, M, C, D, I_{2}, I_{12}, Z), \\ V(t, H, M, C, D, I_{2}, I_{12}, Z), \\ U(t, H, M, C, D, I_{2}, I_{12}, Z)\\ E(t, H, M, C, D, I_{2}, I_{12}, Z)\nonumber\\ \end{cases} \end{eqnarray}$

For the existence theory, we define a Banach space $B = C\times C\times C\times C\times C\times C\times C$ , where $C = [0, T]$ under the norm

$\|\Omega\| = \max\limits_{t\epsilon[0, T]}|H(t)+M(t)+C(t)+D(t)+I_{2}(t)+I_{12}(t), Z(t)|.$

Define an operator $\wp:B \rightarrow B$ as:

$\begin{equation} \wp(\Omega)(t) = \Omega(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t, \Omega(t))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)} \int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(t, \Omega(t))d\omega. \end{equation}$

(3.5)

Now, we suppose that $\Lambda(t, (\Omega(t))$ is the non-linear function with growth and the Lipschitz conditions.

● For each $\Omega\epsilon B, \; \exists$ constants $A_{\Lambda}$ and $P_{\Lambda}$ such that

$\begin{equation} |\Lambda(t, \Omega(t))\leq A_{\Lambda}|\Omega(t)|+P_{\Lambda}. \end{equation}$

(3.6)

● For each $\Omega$ , $\overline{\Omega}\epsilon B, \exists\; a$ constant $K_{\Lambda} > 0$ such that

$\begin{equation} |\Lambda(t, \Omega(t))-\Lambda(t, \overline{\Omega}(t)|\leq K_{\Lambda}|\Omega(t)-\overline{\Omega}(t)|. \end{equation}$

(3.7)

Theorem 2. For the set of continuous function $\Lambda = [0, T]\times B\rightarrow R$ there exists at least single outcome for system (3.1) if the condition (3.6) holds.

Proof.

Let

$L = \ \{\Omega\epsilon B : \|\Omega\|\leq R, R > 0\}.$

Now for any $\Omega\epsilon B$ , we have

$\begin{eqnarray} |\wp(\Omega)| = \max\limits_{t\epsilon[0, T]}|\Omega(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t, \Omega(t))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)} \int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(t, \Omega(t))d\omega \end{eqnarray}$

$\begin{eqnarray} \leq\Omega(0)+\frac{\eta T^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}(A_{\Lambda}\|\Omega\|+P_{\Lambda})+\max\limits_{t\epsilon[0, T]}\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)} \int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}|\Lambda(t, \Omega(t))|d\omega \end{eqnarray}$

$\begin{eqnarray} \leq\Omega(0)+\frac{\eta M^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}(A_{\Lambda}\|\Omega\|+P_{\Lambda})+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}(A_{\Lambda}\|\Omega\|+P_{\Lambda})T^{\xi+\eta-1}L(\xi, \eta) \end{eqnarray}$

$\begin{eqnarray} \leq R. \end{eqnarray}$

Hence, $\wp$ is uniformly bounded, for equicontinuity of $\wp$ , let us take $t_{1} < t_{2}\leq T$ . So, we have

$\begin{eqnarray} |\wp(\Omega)(t_{2})-\wp(\Omega)(t_{1})| = |\frac{\eta t_{2}^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t_{2}, \Omega(t_{2}))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t_{2}}_{0}\omega^{\eta-1}(t_{2}-\omega)^{\eta-1}\Lambda(\omega, \Omega(\omega))d\omega \end{eqnarray}$

$\begin{eqnarray} -\frac{\eta t_{1}^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t_{1}, \Omega(t_{1}))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)} \int^{t_{1}}_{0}\omega^{\eta-1}(t_{1}-\omega)^{\eta-1}\Lambda(\omega, \Omega(\omega))d\omega| \end{eqnarray}$

$\begin{eqnarray} \leq|\frac{\eta t_{2}^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}(A_{\Lambda}|\Omega(t)|, P_{\Lambda})+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}(A_{\Lambda}|\Omega(t)|, P_{\Lambda})t_{2}^{\xi+\eta-1}L(\xi, \eta) \end{eqnarray}$

$\begin{eqnarray} -|\frac{\eta t_{1}^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}(A_{\Lambda}|\Omega(t)|, P_{\Lambda})+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}(A_{\Lambda}|\Omega(t)|, P_{\Lambda})t_{1}^{\xi+\eta-1}L(\xi, \eta), \end{eqnarray}$

when $t_{1}\rightarrow t_{2}$ then $|\wp(\Omega)(t_{2}-\wp(\Omega)(t_{1}|$ . Consequently, we can say that,

$\|\wp(\Omega)(t_{2}-\wp(\Omega)(t_{1}\|\rightarrow 0, t_{1}\rightarrow t_{2}.$

Hence, $\wp$ is equicontinuous. So, theorem is proved. Thus, system has at least one solution.

Theorem 3. Assume that (3.7) holds. If $\varrho < 1$ , where

$\varrho = (\frac{\eta T^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}T^{\xi+\eta-1}L(\xi, \eta))K_{\Lambda}.$

Then the considered model has a unique solution.

Proof. For $\Omega, \overline{\Omega} \epsilon B$ , we have

$\begin{eqnarray} |\wp(\Omega)-\wp(\overline{\Omega})| = \max\limits_{t\epsilon[0, T]}|\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}(\Lambda(t, \Omega(t))-\Lambda(t, \overline{\Omega(t)}))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)} \end{eqnarray}$

$\begin{eqnarray} \int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}d\omega(\Lambda(\omega, \Omega(\omega))-\Lambda(\omega, \overline{\Omega(\omega)}))| \end{eqnarray}$

$\begin{eqnarray} \leq[\frac{\eta T^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}T^{\xi+\eta-1}L(\xi, \eta)]\|\Omega-\overline{\Omega}\| \end{eqnarray}$

$\begin{eqnarray} \leq\varrho\|\Omega-\overline{\Omega}\|. \end{eqnarray}$

Hence, $\wp$ is a contraction. So, system has unique solution by using Banach contraction principle.

Definition 3.1. The proposed model is Ulam-Hyres stable if $\exists\; \wp_{\xi, \eta} \geq 0$ such that for any $\varepsilon > 0$ and for every $\Omega \in (C[0, T], R)$ satisfies the following inequality:

$\begin{eqnarray} |^{FFM}_{0}D^{\xi, \eta}_{t}\Omega(t)-\Lambda(t, \Omega(t))|\leq\varepsilon, t\epsilon[0, T], \end{eqnarray}$

(3.8)

and there exists a unique solution $\Upsilon \in (C[0, T], R)$ such that

$\begin{eqnarray} |\Omega(t)-\Upsilon(t)|\leq\wp_{\xi, \eta}\varepsilon, t\epsilon[0, T]. \end{eqnarray}$

(3.9)

Consider a small perturbation $\Psi \in C[0, T]$ such that $\Psi(0) = 0$ . Let

$\begin{eqnarray} \bullet|\psi(t)|\leq\varepsilon for \varepsilon > 0. \end{eqnarray}$

$\begin{eqnarray} \bullet ^{FFM}_{0}D^{\xi, \eta}_{t}\Omega(t) = \Lambda(t, \Omega(t))+\psi(t). \end{eqnarray}$

Theorem 4. The solution of the perturbed model

$\begin{eqnarray} ^{FFM}_{0}D^{\xi, \eta}_{t}\Omega(t) = \Lambda(t, \Omega(t))+\psi(t), \; \; \; \Omega(0) = \Omega_{0}, \end{eqnarray}$

fulfills the relation given below

$\begin{eqnarray} \label{eq9} |\wp(t)-(\Omega(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}\xi)}\Lambda(t, \Omega(t))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(\omega, \Omega(\omega))d\omega, )| \end{eqnarray}$

$\begin{equation} \leq\xi^{\ast}_{\xi, \eta}\varepsilon, \end{equation}$

(3.10)

where

$\xi^{\ast}_{\xi, \eta} = \frac{\eta T^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}T^{\xi+\eta-1}L(\xi, \eta).$

Proof. Under condition (3.7) along with Theorem 3.5, the solution of the proposed model is Ulam-Hyres stable if $\varrho < 1$ .

Let $\Upsilon\in B$ be a unique solution and $\Omega\in B$ be any solution of the proposed model, then

$\begin{eqnarray} |\Omega(t)-\Upsilon(t)| = |\Omega(t)-(\Upsilon(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t, \Upsilon(t))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)} \end{eqnarray}$

$\begin{eqnarray} \int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(\omega, \Upsilon(\omega))d\omega)| \end{eqnarray}$

$\begin{eqnarray} \leq|\Omega(t)-(\Omega(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t, \Omega(t))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(\omega, \Omega(\omega))d\omega)| \end{eqnarray}$

$\begin{eqnarray} +|\Omega(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t, \Omega(t))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(\omega, \Omega(\omega))d\omega)| \end{eqnarray}$

$\begin{eqnarray} -|\Upsilon(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}\Lambda(t, \Upsilon(t))+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\omega^{\eta-1}(t-\omega)^{\eta-1}\Lambda(\omega, \Upsilon(\omega))d\omega)| \end{eqnarray}$

$\begin{eqnarray} \leq\xi^{\ast}_{\xi, \eta}\varepsilon+(\frac{\eta T^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}T^{\xi+\eta-1}L(\xi, \eta))K_{\Lambda}|\Omega(t)-\Upsilon(t)| \end{eqnarray}$

$\begin{eqnarray} \leq\xi^{\ast}_{\xi, \eta}\varepsilon+\varrho|\Omega(t)-\Upsilon(t)|. \end{eqnarray}$

Consequently, one can write

$\begin{eqnarray} \|\Omega-\Upsilon\|\leq\xi^{\ast}_{\xi, \eta}\varepsilon+\varrho\|\Omega-\Upsilon\|. \end{eqnarray}$

We can write the above relation is

$\begin{eqnarray} \|\Omega-\Upsilon\|\leq\wp_{\xi, \eta}\varepsilon, \end{eqnarray}$

where $\wp_{\xi, \eta} = \frac{\xi^{\ast}_{\xi, \eta}}{1-\varrho}$ . Hence, proposed scheme is Ulam-Hyres stable.

4. Proposed scheme with fractional operator

In this section, The numerical scheme is establish for developed fractional order model, we get

$\begin{eqnarray} \begin{cases} H(t) = H(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{1}(t, H, M, C, D, I_{2}, I_{12}, Z)+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)\\ L_{1}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ M(t) = M(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{2}(t, H, M, C, D, I_{2}, I_{12}, Z)+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)\\ L_{2}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ C(t) = C(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{3}(t, H, M, C, D, I_{2}, I_{12}, Z)+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)\\ L_{3}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ D(t) = D(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{4}(t, H, M, C, D, I_{2}, I_{12}, Z)+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)\\ L_{4}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ I_{2}(t) = I_{2}(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{5}(t, H, M, C, D, I_{2}, I_{12}, Z)+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)\\ L_{5}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ I_{12}(t) = I_{12}(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{6}(t, H, M, C, D, I_{2}, I_{12}, Z)+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)\\ L_{6}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi\\ Z(t) = Z(0)+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{7}(t, H, M, C, D, I_{2}, I_{12}, Z)+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)\\ L_{7}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi. \end{cases} \end{eqnarray}$

(4.1)

Now, we derive the numerical scheme at $t = t_{x+1}$ , we have

$\begin{eqnarray} \begin{cases} H^{x+1} = H^{0}+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{1}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x}) +\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)^{\xi-1} \\L_{1}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ M^{x+1} = M^{0}+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{2}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x}) +\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)^{\xi-1}\\ L_{2}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ C^{x+1} = C^{0}+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{3}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)^{\xi-1}\\ L_{3}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ D^{x+1} = D^{0}+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{4}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)^{\xi-1}\\ L_{4}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ I_{2}^{x+1} = I_{2}^{0}+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{5}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)^{\xi-1}\\ L_{5}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ I_{12}^{x+1} = I_{12}^{0}+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{6}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x}) +\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)^{\xi-1}\\ L_{6}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi\\ Z^{x+1} = Z^{0}+\frac{\eta t^{\eta-1}(1-\xi)}{C_{1}D_{1}(\xi)}L_{7}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\int^{t}_{0}\phi^{\eta-1}(t-\phi)^{\xi-1}\\ L_{7}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi. \end{cases} \end{eqnarray}$

(4.2)

On system (4.2), we have

$\begin{eqnarray} \begin{cases} H^{x+1} = H^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{1}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})\\ +\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\sum^{x}_{f = 0}\int^{t_{f+1}}_{t_{f}}\phi^{\eta-1}(t_{x+1}-\phi)^{\xi-1}L_{1}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ M^{x+1} = M^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{2}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\\ \frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\sum^{x}_{f = 0}\int^{t_{f+1}}_{t_{f}}\phi^{\eta-1}(t_{x+1}-\phi)^{\xi-1}L_{2}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ C^{x+1} = C^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{3}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})\\ +\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\sum^{x}_{f = 0}\int^{t_{f+1}}_{t_{f}}\phi^{\eta-1}(t_{x+1}-\phi)^{\xi-1}L_{3}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ D^{x+1} = D^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{4}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\\ \frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\sum^{x}_{f = 0}\int^{t_{f+1}}_{t_{f}}\phi^{\eta-1}(t_{x+1}-\phi)^{\xi-1}L_{4}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ I_{2}^{x+1} = I_{2}^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{5}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})\\ +\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\sum^{x}_{f = 0}\int^{t_{f+1}}_{t_{f}}\phi^{\eta-1}(t_{x+1}-\phi)^{\xi-1}L_{5}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi, \\ I_{12}^{x+1} = I_{12}^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{6}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})\\ +\frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\sum^{x}_{f = 0}\int^{t_{f+1}}_{t_{f}}\phi^{\eta-1}(t_{x+1}-\phi)^{\xi-1}L_{6}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi\\ Z^{x+1} = Z^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{7}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\\ \frac{\xi\eta}{C_{1}D_{1}(\xi)\Gamma(\xi)}\sum^{x}_{f = 0}\int^{t_{f+1}}_{t_{f}}\phi^{\eta-1}(t_{x+1}-\phi)^{\xi-1}L_{7}(\phi, H, M, C, D, I_{2}, I_{12}, Z)d\phi. \end{cases} \end{eqnarray}$

(4.3)

Using Lagrangian piece-wise interpolation with in the finite interval $[t_f, t_{f+1}]$ , we get

$\begin{eqnarray} \begin{cases} J_{f}(\phi) = \frac{\phi-t_{f-1}}{t_{f}-t_{f-1}}t^{\eta-1}_{f}L_{1}(t^{f}, H^{f}, M^{f}, C^{f}, D^{f}, I_{2}^{f}, I_{12}^{f}, Z^{f})\\ -\frac{\phi-t_{f}}{t_{f}-t_{f-1}}t^{\eta-1}_{f-1}L_{1}(t^{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}), \\ O_{f}(\phi) = \frac{\phi-t_{f-1}}{t_{f}-t_{f-1}}t^{\eta-1}_{f}L_{2}(t^{f}, H^{f}, M^{f}, C^{f}, D^{f}, I_{2}^{f}, I_{12}^{f}, Z^{f})\\ -\frac{\phi-t_{f}}{t_{f}-t_{f-1}}t^{\eta-1}_{f-1}L_{2}(t^{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}), \\ G_{f}(\phi) = \frac{\phi-t_{f-1}}{t_{f}-t_{f-1}}t^{\eta-1}_{f}L_{3}(t^{f}, H^{f}, M^{f}, C^{f}, D^{f}, I_{2}^{f}, I_{12}^{f}, Z^{f})\\ -\frac{\phi-t_{f}}{t_{f}-t_{f-1}}t^{\eta-1}_{f-1}L_{3}(t^{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}), \\ N_{f}(\phi) = \frac{\phi-t_{f-1}}{t_{f}-t_{f-1}}t^{\eta-1}_{f}L_{4}(t^{f}, H^{f}, M^{f}, C^{f}, D^{f}, I_{2}^{f}, I_{12}^{f}, Z^{f})\\ -\frac{\phi-t_{f}}{t_{f}-t_{f-1}}t^{\eta-1}_{f-1}L_{4}(t^{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}), \\ F_{f}(\phi) = \frac{\phi-t_{f-1}}{t_{f}-t_{f-1}}t^{\eta-1}_{f}L_{5}(t^{f}, H^{f}, M^{f}, C^{f}, D^{f}, I_{2}^{f}, I_{12}^{f}, Z^{f})\\ -\frac{\phi-t_{f}}{t_{f}-t_{f-1}}t^{\eta-1}_{f-1}L_{5}(t^{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}), \\ Q_{f}(\phi) = \frac{\phi-t_{f-1}}{t_{f}-t_{f-1}}t^{\eta-1}_{f}L_{6}(t^{f}, H^{f}, M^{f}, C^{f}, D^{f}, I_{2}^{f}, I_{12}^{f}, Z^{f})\\ -\frac{\phi-t_{f}}{t_{f}-t_{f-1}}t^{\eta-1}_{f-1}L_{6}(t^{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}), \\ R_{f}(\phi) = \frac{\phi-t_{f-1}}{t_{f}-t_{f-1}}t^{\eta-1}_{f}L_{7}(t^{f}, H^{f}, M^{f}, C^{f}, D^{f}, I_{2}^{f}, I_{12}^{f}, Z^{f})\\ -\frac{\phi-t_{f}}{t_{f}-t_{f-1}}t^{\eta-1}_{f-1}L_{7}(t^{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}). \end{cases} \end{eqnarray}$

(4.4)

We have

$\begin{eqnarray} \begin{cases} H^{x+1} = H^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{1}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\eta(\Theta t)^{\xi}}{C_{1}D_{1}(\xi)\Gamma(\xi+2)} \Sigma^{x}_{f = 0}\\ \Big[t^{\eta-1}_{f}L_{1}(t_{f}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}) \times((x+1-f)^{\xi}(x-f+2+\xi)\\ -(x-2)^{\xi}(2+2\xi+x-2))-t^{\eta-1}_{f-1}L_{1}(t_{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})\\ \times((1+x-f)^{\xi+1}-(x-f)^{\xi}(1+\xi+x-f))\Big], \\ M^{x+1} = M^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{2}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\eta(\Theta t)^{\xi}}{C_{1}D_{1}(\xi)\Gamma(\xi+2)} \Sigma^{x}_{ff = 0} \\ \Big[t^{\eta-1}_{f}L_{2}(t_{f}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})\times((x+1-f)^{\xi}(x-f+2+\xi)\\ -(x-2)^{\xi}(2+2\xi+x-2))-t^{\eta-1}_{f-1}L_{2}(t_{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}) \\ \times((1+x-f)^{\xi+1}-(x-f)^{\xi}(1+\xi+x-f))\Big], \\ C^{x+1} = C^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{3}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\eta(\Theta t)^{\xi}}{C_{1}D_{1}(\xi)\Gamma(\xi+2)} \Sigma^{x}_{f = 0}\\ \Big[t^{\eta-1}_{f}L_{3}(t_{f}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})\times((x+1-f)^{\xi}(x-f+2+\xi)\\ -(x-2)^{\xi}(2+2\xi+x-2))-t^{\eta-1}_{f-1}L_{3}(t_{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})\\ \times((1+x-f)^{\xi+1}-(x-f)^{\xi}(1+\xi+x-f))\Big], \\ D^{x+1} = D^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{4}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\eta(\Theta t)^{\xi}}{C_{1}D_{1}(\xi)\Gamma(\xi+2)} \Sigma^{x}_{f = 0} \\ \Big[t^{\eta-1}_{f}L_{4}(t_{f}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}) \times((x+1-f)^{\xi}(x-f+2+\xi) \\ -(x-2)^{\xi}(2+2\xi+x-2))-t^{\eta-1}_{f-1}L_{4}(t_{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})\\ \times((1+x-f)^{\xi+1}-(x-f)^{\xi}(1+\xi+x-f))\Big], \\ I_{2}^{x+1} = I_{2}^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{5}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\eta(\Theta t)^{\xi}}{C_{1}D_{1}(\xi)\Gamma(\xi+2)}\Sigma^{x}_{f = 0}\\ \Big[t^{\eta-1}_{f}L_{5}(t_{f}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}) \times((x+1-f)^{\xi}(x-f+2+\xi)\\ -(x-2)^{\xi}(2+2\xi+x-2))-t^{\eta-1}_{f-1}L_{5}(t_{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}) \\ \times((1+x-f)^{\xi+1}-(x-f)^{\xi}(1+\xi+x-f))\Big], \\ I_{12}^{x+1} = I_{12}^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{6}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\eta(\Theta t)^{\xi}}{C_{1}D_{1}(\xi)\Gamma(\xi+2)} \Sigma^{x}_{f = 0}\\ \Big[t^{\eta-1}_{f}L_{6}(t_{f}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})times((x+1-f)^{\xi}(x-f+2+\xi)\\ -(x-2)^{\xi}(\xi2+2\xi+x-2))-t^{\eta-1}_{f-1}L_{6}(t_{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})\\ \times((1+x-f)^{\xi+1}-(x-f)^{\xi}(1+\xi+x-f))\Big], \\ Z^{x+1} = Z^{0}+\frac{\eta t^{\eta-1}_{x}(1-\xi)}{C_{1}D_{1}(\xi)}L_{7}(t_{x}, H^{x}, M^{x}, C^{x}, D^{x}, I_{2}^{x}, I_{12}^{x}, Z^{x})+\frac{\eta(\Theta t)^{\xi}}{C_{1}D_{1}(\xi)\Gamma(\xi+2)} \Sigma^{x}_{f = 0}\\ \Big[t^{\eta-1}_{f}L_{7}(t_{f}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1}) \times((x+1-f)^{\xi}(x-f+2+\xi)\\ -(x-2)^{\xi}(2+2\xi+x-2))-t^{\eta-1}_{f-1}L_{7}(t_{f-1}, H^{f-1}, M^{f-1}, C^{f-1}, D^{f-1}, I_{2}^{f-1}, I_{12}^{f-1}, Z^{f-1})\\ \times((1+x-f)^{\xi+1}-(x-f)^{\xi}(1+\xi+x-f))\Big]. \end{cases} \end{eqnarray}$

(4.5)

5. Numerical results and discussion

The mathematical analysis of the tumor-immune system consisting of a fractional differential equation has been presented. In order to conduct a credible study, several $\alpha$ values are evaluated. By using non-integer parameter values, surprising reactions are achieved from the compartments of the proposed fractional-order model. Solution for $H(t), M(t), C(t)$ , $I_{2}(t)$ and $I_{12}(t)$ start decreasing by decreasing the fractional values while $D(t)$ and $Z(t)$ start increasing by decreasing fractional values respectively and can be seen easily from – which converges to steady state with dimension 0.7. Solution for $H(t), M(t), D(t)$ , $I_{2}(t)$ , $I_{12}(t)$ and $Z(t)$ start increasing by decreasing the fractional values while $C(t)$ start decreasing by decreasing fractional values respectively and can be seen easily from Figures 8–14 which converges to steady state with dimension 0.9. It is easily observed from Figure 3 and Figure 9 that cancer cell start decreasing by decreasing the fractional with both dimensions. But Fractal-fractional technique with dimension 0.9 gives better results for all the compartments which converges rapidly to steady state as compare to dimension 0.7. It is easily observed from simulation that fractional order model check the relationship between cancer cells and immune system which further increase the immune system of CD4+T and CD8+T lymphocytes production and decrease the cancer cell by adding IL-12 cytokine and anti-PD-L1 inhibitor. The variable anti-PD-L1 inhibitor was introduced to catch cells escaping from the immune system and destroy these cells which can be seen in Figures 7 and 14 with both dimensions. Also the variable IL-12 cytokine is used to increase the number of CD4+T and CD8+T lymphocytes production which can be seen in Figures 1, 2 and 8, 9 with dimension 0.7 and 0.9 respectively. This kind of study helps to analyze the effect of cancer cell and its anti-PD-L1 inhibitor in human body. Numerical simulations are important in sense for developing algorithm of treatment cell and impact of vaccination to control the disease. This kind of study is helpful for physician for planing and decision making for treatment cells.

Figure 1.

$H(t)$ with fractal fractional derivative at

$\xi = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 2.

$M(t)$ with fractal fractional derivative at

$\xi = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 3.

$C(t)$ with fractal fractional derivative at

$\xi = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 4.

$D(t)$ with fractal fractional derivative at

$\xi = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 5.

$I_2(t)$ with fractal fractional derivative at

$\xi = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 6.

$I_{12}(t)$ with fractal fractional derivative at

$\xi = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 7.

$Z(t)$ with fractal fractional derivative at

$\xi = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 8.

$H(t)$ with fractal fractional derivative at

$\xi = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 9.

$M(t)$ with fractal fractional derivative at

$\xi = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 10.

$C(t)$ with fractal fractional derivative at

$\xi = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 11.

$D(t)$ with fractal fractional derivative at

$\xi = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 12.

$I_2(t)$ with fractal fractional derivative at

$\xi = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 13.

$I_{12}(t)$ with fractal fractional derivative at

$\xi = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 14.

$Z(t)$ with fractal fractional derivative at

$\xi = 0.9$ .

DownLoad: Full-Size Img PowerPoint

6. Conclusions

We developed a scheme for the analytical solution of the fractional order cancer model by using the Fractal Fractional operator. The model represents population dynamics during the disease as a set of non-linear fractional-order ordinary differential equations. The Fractional-order system is analyzed qualitatively to verify the state of the disease as well as determine the unique positive solutions. Uniqueness and boundedness for solution are proved using the fixed point theory. The solution is obtained for the fractional-order cancer model using the Fractal Fractional operator. Numerical simulations are carried out to check the actual behavior of the cancer outbreak using a developed scheme for fractional-order model, which will be helpful in future understanding of this disease and control strategics. The simulation easily observed that the fractional-order model examines that IL-12 cytokine and anti-PD-L1 inhibitor are used to check the relationship between the immune system and cancer cells, which cause to increase the immune system and decrease the cancer cell is due to IL-12 cytokine and anti-PD-L1 inhibitor with the impact of fractional parameters. Results show significant changes using different fractional values with different dimensions. Numerical simulations are important for developing an algorithm of treatment cells and the impact of vaccination to control the disease. In the future, the other cytokineses are also added to the mathematical model, which can further check with different types of derivatives and different types of techniques.

Conflict of interest

All authors declare no conflicts of interest in this paper.

References

[1]	Alkış Küçükaydın, M., Predictors of climate change conspiracy beliefs collective narcissism paranoid thoughts and distrust in science. Current Psychology, 2024, 43: 34717–34728. https://dx.doi.org/10.1007/s12144-024-06975-6 doi: 10.1007/s12144-024-06975-6
[2]	Amabile, T.M., A model of creativity and innovation in organizations. Research in Organizational Behavior, 1988, 10(1): 123-167.
[3]	Baer, J. and Kaufman, J.C., Gender differences in creativity. The Journal of Creative Behavior, 2008, 42(2): 75-105. https://doi.org/10.1002/j.2162-6057.2008.tb01289.x doi: 10.1002/j.2162-6057.2008.tb01289.x
[4]	Barbot, B., Besançon, M. and Lubart, T., Creative potential in educational settings: Its nature, measure, and nurture. Creativity and Creative Pedagogies in the Early and Primary Years, 2018, 12-22. https://doi.org/10.1080/03004279.2015.1020643
[5]	Beaton, D.E., Bombardier, C., Guillemin, F. and Ferraz, M.B., Guidelines for the process of cross-cultural adaptation of self-report measures. Spine, 2000, 25(24): 3186–3191. https://doi.org/10.1097/00007632-200012150-00014 doi: 10.1097/00007632-200012150-00014
[6]	Besançon, M., Lubart, T. and Barbot, B., Creative giftedness and educational opportunities. Educational and Child Psychology, 2013, 30(2): 79-88. https://doi.org/10.53841/bpsecp.2013.30.2.79 doi: 10.53841/bpsecp.2013.30.2.79
[7]	Brophy, D.R., Understanding, measuring, enhancing collective creative problem-solving efforts. Creativity Research Journal, 1998, 11(3): 199-229. https://doi.org/10.1207/s15326934crj1103_2 doi: 10.1207/s15326934crj1103_2
[8]	Castelijns, J., Vermeulen, M. and Kools, Q., Collective learning in primary schools and teacher education institutes. Journal of Educational Change, 2013, 14: 373-402. https://doi.org/10.1007/s10833-013-9209-6 doi: 10.1007/s10833-013-9209-6
[9]	Claxton, A.F., Pannells, T.C. and Rhoads, P.A., Developmental trends in the creativity of school-age children. Creativity Research Journal, 2005, 17(4): 327-335. https://doi.org/10.1207/s15326934crj1704_4 doi: 10.1207/s15326934crj1704_4
[10]	Davies, D., Jindal-Snape, D., Collier, C., Digby, R., Hay, P. and Howe, A., Creative learning environments in education—A systematic literature review. Thinking Skills and Creativity, 2013, 8: 80-91. https://doi.org/10.1016/j.tsc.2012.07.004 doi: 10.1016/j.tsc.2012.07.004
[11]	Dikici, A., The adaptation of creativity fostering primary teachers index scale into Turkish. Educational Sciences: Theory and Practice, 2013, 13(1): 318-323.
[12]	Ekvall, G., Organizational climate for creativity and innovation. European Journal of Work and Organizational Psychology, 1996, 5(1): 105–123. https://doi.org/10.1080/13594329608414845 doi: 10.1080/13594329608414845
[13]	Fan, M. and Cai, W., How does a creative learning environment foster student creativity? An examination on multiple explanatory mechanisms. Current Psychology, 2022, 41(7): 4667-4676. https://doi.org/10.1007/s12144-020-00974-z doi: 10.1007/s12144-020-00974-z
[14]	Ferguson, E. and Cox, T., Exploratory factor analysis: A users' guide. International Journal of Selection and Assessment, 1993, 1(2): 84–94. https://doi.org/10.1111/j.1468-2389.1993.tb00092.x doi: 10.1111/j.1468-2389.1993.tb00092.x
[15]	Fornell, C. and Larcker, D.F., Structural equation models with unobservable variables and measurement error: Algebra and statistics. Journal of Marketing Research, 1981, 18(3): 382–388. https://doi.org/10.1177/002224378101800313 doi: 10.1177/002224378101800313
[16]	Gefen, D., Straub, D. and Boudreau, M.C., Structural equation modeling and regression: Guidelines for research practice. Communications of the Association for Information Systems, 2000, 4(1): 1-70. https://doi.org/10.17705/1CAIS.00407 doi: 10.17705/1CAIS.00407
[17]	Hadzigeorgiou, Y., Fokialis, P. and Kabouropoulou, M., Thinking about creativity in science education. Creative Education, 2012, 3(5): 603–611. http://dx.doi.org/10.4236/ce.2012.350 doi: 10.4236/ce.2012.350
[18]	Hair, J.F., Anderson, R.E., Tatham, R.L. and Black W.C., Multivariate data analysis (5th edition), Prentice-Hall, 1998.
[19]	Hargadon, A.B. and Bechky, B.A., When collections of creatives become creative collectives: A field study of problem solving at work. Organization Science, 2006, 17(4): 484–500. https://doi.org/10.1287/orsc.1060.0200 doi: 10.1287/orsc.1060.0200
[20]	Harris, J.A., Measured intelligence, achievement, openness to experience, and creativity. Personality and Individual Differences, 2004, 36(4): 913–929. https://doi.org/10.1016/S0191-8869(03)00161-2 doi: 10.1016/S0191-8869(03)00161-2
[21]	Hayes, A.F. and Coutts, J.J., Use omega rather than Cronbach's alpha for estimating reliability. But….Communication Methods and Measures, 2020, 14(1): 1-24. https://doi.org/10.1080/19312458.2020.1718629
[22]	He, W.J. and Wong, W.C., Gender differences in creative thinking revisited: Findings from analysis of variability. Personality and Individual Differences, 2011, 51(7): 807-811. https://doi.org/10.1016/j.paid.2011.06.027 doi: 10.1016/j.paid.2011.06.027
[23]	Hong, O. and Song, J., A componential model of Science Classroom Creativity (SCC) for understanding collective creativity in the science classroom. Thinking Skills and Creativity, 2020, 37: 100698. https://doi.org/10.1016/j.tsc.2020.100698 doi: 10.1016/j.tsc.2020.100698
[24]	Hong, O., Park, M.H. and Song, J., The assessment of science classroom creativity: Scale development. International Journal of Science Education, 2022, 44(8): 1356-1377. https://doi.org/10.1080/09500693.2022.2077466 doi: 10.1080/09500693.2022.2077466
[25]	Hu, W. and Adey, P., A scientific creativity test for secondary school students. International Journal of Science Education, 2002, 24(4): 389-403. https://doi.org/10.1080/09500690110098912 doi: 10.1080/09500690110098912
[26]	Huang, X., Chin-Hsi, L., Mingyao, S. and Peng, X., What drives teaching for creativity? Dynamic componential modelling of the school environment, teacher enthusiasm, and metacognition. Teaching and Teacher Education, 2021,107: 103491. https://doi.org/10.1016/j.tate.2021.103491 doi: 10.1016/j.tate.2021.103491
[27]	Jo, S., A validation study of the creative school environment perceptions scale and a study of group differences. Journal of Gifted/Talented Education, 2012, 22(3): 663–677. http://dx.doi.org/10.9722/JGTE.2012.22.3.663
[28]	Kandemir, M.A. and Kaufman, J.C., The Kaufman domains of creativity scale: Turkish validation and relationship to academic major. The Journal of Creative Behavior, 2020, 54(4): 1002-1012. https://doi.org/10.1002/jocb.428 doi: 10.1002/jocb.428
[29]	Kang, E.J., A multilevel analysis of factors affecting kindergartners' creative dispositions in relations to child-level variables and teacher-level variables. International Journal of Child Care and Education Policy, 2020, 14(1): 1-17. https://doi.org/10.1186/s40723-020-00077-z doi: 10.1186/s40723-020-00077-z
[30]	Karaca, N.H. and Aral, N., Implementation of thinking creatively in action and movement test for Turkish children. Mehmet Akif Ersoy University Facult of Education Journal, 2017, 42: 240-253. https://doi.org/10.21764/efd.26968 doi: 10.21764/efd.26968
[31]	Karakelle, S. and Saraç, S., Creative mindsets scale: Turkish validity and reliability. Cyprus Turkish Journal of Psychiatry & Psychology, 2022, 4(2): 122-129. https://doi.org/10.35365/ctjpp.22.2.01 doi: 10.35365/ctjpp.22.2.01
[32]	Kaufman, J.C. and Beghetto, R.A., Beyond big and little: The four c model of creativity. Review of General Psychology, 2009, 13(1): 1-12. https://doi.org/10.1037/a0013688 doi: 10.1037/a0013688
[33]	Kelloway, E.K., Using Mplus for structural equation modeling (2nd ed.), Sage, 2015.
[34]	Kline, R.B., Principles and practice of structural equation modeling (3rd ed.), Guilford Press, 2011.
[35]	Lau, S. and Cheung, P.C., Developmental trends of creativity: What twists of turn do boys and girls take at different grades?. Creativity Research Journal, 2010, 22(3): 329-336. https://doi.org/10.1080/10400419.2010.503543 doi: 10.1080/10400419.2010.503543
[36]	Liu, S.C. and Lin, H.S., Primary teachers' beliefs about scientific creativity in the classroom context. International Journal of Science Education, 2014, 36(10): 1551-1567. https://doi.org/10.1080/09500693.2013.868619 doi: 10.1080/09500693.2013.868619
[37]	Lubart, T.I. and Lautrey, J., Relationships between creative development and cognitive development, Paper presented at the Seventh European Conference on Developmental Psychology, Krakow, Poland, 1995.
[38]	Matud, M.P., Rodríguez, C. and Grande, J., Gender differences in creative thinking. Personality and Individual Differences, 2007, 43(5): 1137-1147. https://doi.org/10.1016/j.paid.2007.03.006 doi: 10.1016/j.paid.2007.03.006
[39]	McDonald, R.P., Test theory: A unified treatment, Psychology Press, 2013.
[40]	McLean, L.D., Organizational culture's influence on creativity and innovation: A review of the literature and implications for human resource development. Advances in Developing Human Resources, 2005, 7(2): 226-246. https://doi.org/10.1177/1523422305274528 doi: 10.1177/1523422305274528
[41]	Miller, A.L., A self-report measure of cognitive processes associated with creativity. Creativity Research Journal, 2014, 26(2): 203-218. https://doi.org/10.1080/10400419.2014.901088 doi: 10.1080/10400419.2014.901088
[42]	Min, J. and Seo, E., The relationship among creative climate, motivation and creativity. The Korean Journal of Educational Psychology, 2009, 23(4): 787–800.
[43]	Organisation for Economic Co-operation and Development (OECD), The future of education and skills: Education 2030, OECD, 2018. Retrieved from: http://www.oecd.org/education/2030project/contact/E2030
[44]	Osborne, J., Collins, S., Ratcliffe, M., Millar, R. and Duschl, R., What? Ideas-about-science? Should be taught in school science? A Delphi study of the expert community. Journal of Research in Science Teaching, 2003, 40(7): 692–720. https://doi.org/10.1002/tea.10105 doi: 10.1002/tea.10105
[45]	Puccio, G.J., Murdock, M C. and Mance, M., Creative leadership: Skills that drive change, Sage Publications, 2007.
[46]	Rhodes, M., An analysis of creativity. The Phi Delta Kappan, 1961, 42(7): 305–310.
[47]	Richardson, C. and Mishra, P., Learning environments that support student creativity: Developing the SCALE. Thinking Skills and Creativity, 2018, 27: 45-54. https://doi.org/10.1016/j.tsc.2017.11.004 doi: 10.1016/j.tsc.2017.11.004
[48]	Saleem Masadeh, T., Teaching practices of EFL teachers and the enhancement of creative thinking skills among learners. International Journal of Asian Education, 2021, 2(2): 153–166.
[49]	Sawyer, R.K., Explaining creativity: The science of human innovation, Oxford University Press, 2011.
[50]	Sen, S. and Yö rük, S., Adaptation of the Runco ideational behavior scale into Turkish: A confirmatory factor analysis and Rasch study. Psychology of Aesthetics, Creativity, and the Arts, 2024, Advance online publication. https://doi.org/10.1037/aca0000676
[51]	Soh, K., Fostering student creativity through teacher behaviors. Thinking Skills and Creativity, 2017, 23: 58-66. https://doi.org/10.1016/j.tsc.2016.11.002 doi: 10.1016/j.tsc.2016.11.002
[52]	Sung, E., Kim, N. and Park, S., An exploratory study for the validation of a questionnaire of creative class environment- the case of high school students. The Journal of Creativity Education, 2014, 14(3): 1–15.
[53]	Swanzy-Impraim, E., Morris, J.E., Lummis, G.W. and Jones, A., Promoting creativity: Secondary visual art teachers' perceptions and understanding of creativity in Ghana. Thinking Skills and Creativity, 2022, 45: 101057. https://doi.org/10.1016/j.tsc.2022.101057 doi: 10.1016/j.tsc.2022.101057
[54]	Tabachnick, B.G. and Fidell, L.S., Using multivariate statistics (6th ed.), Pearson Education, 2013.
[55]	Tan, J.P.L., Caleon, I., Ng, H.L., Poon, C.L. and Koh, E. Collective creativity competencies and collaborative problem-solving outcomes: Insights from the dialogic interactions of Singapore student teams. In E. Care, P. Griffin, & M. Wilson (Eds.), Assessment and teaching of 21st century skills: Research and applications, 2018, 95–118. Springer.
[56]	Tang, T., Vezzani, V. and Eriksson, V., Developing critical thinking, collective creativity skills and problem solving through playful design jams. Thinking Skills and Creativity, 2020, 37: 100696. https://doi.org/10.1016/j.tsc.2020.100696 doi: 10.1016/j.tsc.2020.100696
[57]	Timotheou, S. and Ioannou, A., Collective creativity in STEAM making activities. The Journal of Educational Research, 2021,114(2): 130-138. https://doi.org/10.1080/00220671.2021.1873721 doi: 10.1080/00220671.2021.1873721
[58]	Urban, K.K. and Jellen, H.G., Test for Creative Thinking –Drawing Production (TCT-DP), Swets Test Services, 1996.
[59]	Vygotsky, L.S., Imagination and creativity in childhood. Journal of Russian and East European Psychology, 2004, 42(1): 7–97. https://doi.org/10.1080/10610405.2004.11059210 doi: 10.1080/10610405.2004.11059210
[60]	Woodman, R.W. and Schoenfeldt, L.F., An interactionist model of creative behavior. The Journal of Creative Behavior, 1990, 24(4): 279–290. https://doi.org/10.1002/j.2162-6057.1990.tb00525.x doi: 10.1002/j.2162-6057.1990.tb00525.x
[61]	Yıldız Ç içekler, C., Pirpir, D.A. and Aral, N., Turkish standardization of early childhood creativity scale. İlkö ğretim Online, 2020,817-830. https://doi.org/10.17051/ilkonline.2020.695260
[62]	Yıldız, N.O., Kaçay, Z. and Soyer, F., Adaptation of individual creativity scale to Turkish culture, validity and reliability study. Sakarya University Journal of Education, 2022, 12(3): 712-729. https://doi.org/10.19126/suje.1064436 doi: 10.19126/suje.1064436
[63]	Yin, Y., Wang, P. and Childs, P., Understanding creativity process through electroencephalography measurement on creativity-related cognitive factors. Frontiers in Neuroscience, 2022, 16: 951272. https://doi.org/10.3389/fnins.2022.951272 doi: 10.3389/fnins.2022.951272

This article has been cited by:

1.	Sugandha Arora, Trilok Mathur, Kamlesh Tiwari, A fractional-order model to study the dynamics of the spread of crime, 2023, 426, 03770427, 115102, 10.1016/j.cam.2023.115102
2.	Shahram Rezapour, Chernet Tuge Deressa, Robert G. Mukharlyamov, Sina Etemad, Watcharaporn Cholamjiak, On a Mathematical Model of Tumor-Immune Interaction with a Piecewise Differential and Integral Operator, 2022, 2022, 2314-4785, 1, 10.1155/2022/5075613
3.	Kottakkaran Sooppy Nisar, Muhammad Shoaib, Muhammad Asif Zahoor Raja, Rafia Tabassum, Ahmed Morsy, A novel design of evolutionally computing to study the quarantine effects on transmission model of Ebola Virus Disease, 2023, 22113797, 106408, 10.1016/j.rinp.2023.106408
4.	CHERNET TUGE DERESSA, ROBERT G. MUKHARLYAMOV, HOSSAM A. NABWEY, SINA ETEMAD, İBRAHIM AVCI, ON THE CHAOTIC NATURE OF A CAPUTO FRACTIONAL MATHEMATICAL MODEL OF CANCER AND ITS CROSSOVER BEHAVIORS, 2024, 32, 0218-348X, 10.1142/S0218348X2440053X
5.	Huda Alsaud, Muhammad Owais Kulachi, Aqeel Ahmad, Mustafa Inc, Muhammad Taimoor, Investigation of SEIR model with vaccinated effects using sustainable fractional approach for low immune individuals, 2024, 9, 2473-6988, 10208, 10.3934/math.2024499
6.	Kottakkaran Sooppy Nisar, Muhammad Farman, Mahmoud Abdel-Aty, Jinde Cao, A review on epidemic models in sight of fractional calculus, 2023, 75, 11100168, 81, 10.1016/j.aej.2023.05.071
7.	Anusmita Das, Kaushik Dehingia, Evren Hinçal, Fatma Özköse, Kamyar Hosseini, A study on the dynamics of a breast cancer model with discrete-time delay, 2024, 99, 0031-8949, 035235, 10.1088/1402-4896/ad2753
8.	ZULQURNAIN SABIR, DUMITRU BALEANU, MUHAMMAD ASIF ZAHOOR RAJA, ALI S. ALSHOMRANI, EVREN HINCAL, MEYER WAVELET NEURAL NETWORKS PROCEDURES TO INVESTIGATE THE NUMERICAL PERFORMANCES OF THE COMPUTER VIRUS SPREAD WITH KILL SIGNALS, 2023, 31, 0218-348X, 10.1142/S0218348X2340025X
9.	Muhammad Farman, Maryam Batool, Kottakkaran Sooppy Nisar, Abdul Sattar Ghaffari, Aqeel Ahmad, Controllability and analysis of sustainable approach for cancer treatment with chemotherapy by using the fractional operator, 2023, 51, 22113797, 106630, 10.1016/j.rinp.2023.106630
10.	Parvaiz Ahmad Naik, Muhammad Owais Kulachi, Aqeel Ahmad, Muhammad Farman, Faiza Iqbal, Muhammad Taimoor, Zhengxin Huang, Modeling different strategies towards control of lung cancer: leveraging early detection and anti-cancer cell measures, 2024, 1025-5842, 1, 10.1080/10255842.2024.2404540
11.	Muhammad Farman, Cicik Alfiniyah, A constant proportional caputo operator for modeling childhood disease epidemics, 2024, 10, 27726622, 100393, 10.1016/j.dajour.2023.100393
12.	Aqeel Ahmad, Muhammad Owais Kulachi, Muhammad Farman, Moin-ud-Din Junjua, Muhammad Bilal Riaz, Sidra Riaz, Muntazir Hussain, Mathematical modeling and control of lung cancer with IL2 cytokine and anti-PD-L1 inhibitor effects for low immune individuals, 2024, 19, 1932-6203, e0299560, 10.1371/journal.pone.0299560
13.	ZULQURNAIN SABIR, DUMITRU BALEANU, MUHAMMAD ASIF ZAHOOR RAJA, ALI S. ALSHOMRANI, EVREN HINCAL, COMPUTATIONAL PERFORMANCES OF MORLET WAVELET NEURAL NETWORK FOR SOLVING A NONLINEAR DYNAMIC BASED ON THE MATHEMATICAL MODEL OF THE AFFECTION OF LAYLA AND MAJNUN, 2023, 31, 0218-348X, 10.1142/S0218348X23400169
14.	Kaushik Dehingia, Salah Boulaaras, The Stability of a Tumor–Macrophages Model with Caputo Fractional Operator, 2024, 8, 2504-3110, 394, 10.3390/fractalfract8070394
15.	Rongrong Qiao, Yuhan Hu, Dynamic analysis of a $$\pmb {SI_{1}I_{2}AD}$$ information dissemination model considering the word of mouth, 2023, 111, 0924-090X, 22763, 10.1007/s11071-023-09021-5
16.	Kottakkaran Sooppy Nisar, Muhammad Owais Kulachi, Aqeel Ahmad, Muhammad Farman, Muhammad Saqib, Muhammad Umer Saleem, Fractional order cancer model infection in human with CD8+ T cells and anti-PD-L1 therapy: simulations and control strategy, 2024, 14, 2045-2322, 10.1038/s41598-024-66593-x
17.	Kaushik Dehingia, Sana Abdulkream Alharbi, Awatif Jahman Alqarni, Mounirah Areshi, Mona Alsulami, Reima Daher Alsemiry, Reem Allogmany, Homan Emadifar, Mati ur Rahman, Kavikumar Jacob, Exploring the combined effect of optimally controlled chemo-stem cell therapy on a fractional-order cancer model, 2025, 20, 1932-6203, e0311822, 10.1371/journal.pone.0311822
18.	Kaushik Dehingia, Bhagya Jyoti Nath, 2025, Chapter 5, 978-981-97-8714-2, 99, 10.1007/978-981-97-8715-9_5
19.	Chinwe Peace Igiri, Samuel Shikaa, Monkeypox Transmission Dynamics Using Fractional Disease Informed Neural Network: A Global and Continental Analysis, 2025, 13, 2169-3536, 77611, 10.1109/ACCESS.2025.3559005
20.	Jianping Li, Nan Liu, Danni Wang, Hongli Yang, Mathematical modeling and Hopf bifurcation analysis of tumor macrophage interaction with polarization delay, 2025, 19, 1751-3758, 10.1080/17513758.2025.2508240

Author's biography Doç. Dr. Menşure Alkış Küçükaydın is a researcher and faculty member at Necmettin Erbakan University. She holds a degree in science education. Her research focuses on empirical studies in science and science education. She includes critical learning theories in her studies. Dr.Alkış Küçükaydın is also interested in the scientific inquiry, computational thinking, scientist image, and science teacher education; Dr. Çiğdem Akkanat Avşar is a researcher and works in the Ministry of National Education in Turkey. She holds a degree in science education. Her research focuses on science education and creativity. Dr.Akkanat Avşar is interested in scientific thinking, special education, scientist history, and science teacher education

Reader Comments

Your name:*

Email:*
© 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

STEM Education

2.2

Metrics

Article views(761) PDF downloads(32) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Tables(4)

STEM Education

Collective creativity in science classrooms: Scale adaptation and an investigation in terms of demographic variables

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Cancer model with fractal-fractional operator

4. Proposed scheme with fractional operator

5. Numerical results and discussion

6. Conclusions

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

STEM Education

Collective creativity in science classrooms: Scale adaptation and an investigation in terms of demographic variables

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Cancer model with fractal-fractional operator

4. Proposed scheme with fractional operator

5. Numerical results and discussion

6. Conclusions

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog