Next Article in Journal
Bayes Inference of Structural Safety under Extreme Wind Loads Based upon a Peak-Over-Threshold Process of Exceedances
Next Article in Special Issue
Continuum and Molecular Modeling of Chemical Vapor Deposition at Nano-Scale Fibrous Substrates
Previous Article in Journal
Asymptotic Consideration of Rayleigh Waves on a Coated Orthorhombic Elastic Half-Space Reinforced Using an Elastic Winkler Foundation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Explicit Integrating Factor Runge–Kutta Method for the Extended Fisher–Kolmogorov Equation

Fujian Province University Key Laboratory of Computation Science, School of Mathematical Sciences, Huaqiao University, Quanzhou 362021, China
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2023, 28(6), 110; https://doi.org/10.3390/mca28060110
Submission received: 8 October 2023 / Revised: 20 November 2023 / Accepted: 20 November 2023 / Published: 22 November 2023

Abstract

:
The extended Fisher–Kolmogorov (EFK) equation is an important model for phase transitions and bistable phenomena. This paper presents some fast explicit numerical schemes based on the integrating factor Runge–Kutta method and the Fourier spectral method to solve the EFK equation. The discrete global convergence of these new schemes is analyzed rigorously. Three numerical examples are presented to verify the theoretical analysis and the efficiency of the proposed schemes.

1. Introduction

The Fisher–Kolmogorov (FK) equation was first proposed by Fisher [1] and Kolmogorov [2] in 1937 to describe the interaction between the spread and adaptation of biological populations. By adding a stabilizing fourth-order derivative term to the FK equation, Collent [3] and Saarloos [4] proposed the extended FK (EFK) equation, which is a very important mathematical and physical model and has been widely used in many physics and engineering applications. In this paper, we consider the following EFK model with periodic boundary conditions:
t v = κ Δ 2 v + Δ v + g ( v ) , ( x , t ) Ω × ( 0 , T ] , v ( x , 0 ) = v 0 ( x ) , x Ω ,
where Ω [ a , b ] d , ( d = 1 , 2 ) is a bounded area and κ is a positive constant. The function g ( v ) = G ( v ) and G ( v ) = 1 4 ( v 2 1 ) 2 is a double-well potential. When κ = 0 in (1), the EFK model reduces to the classical FK model. We assume that the function g ( v ) exhibits Lipschitz continuity with respect to Ω , where the Lipschitz constant L is defined as follows
max v R | g ( v ) | L .
The EFK model (1) can be viewed as an L 2 gradient flow associated with the following energy
E ( v ) = Ω κ 2 | Δ v | 2 + 1 2 | v | 2 + G ( v ) d x ,
which is diminishing in time, i.e., d d t E ( v ) 0 .
All the above properties are determined by the inherent nature of the physical model. Thus, in order to avoid non-physical effects in the simulations over a long period of time, it is highly desirable to design a structure-preserving numerical scheme. There have been some excellent results in numerical research on the EFK model, such as the spline configuration method [5], nonlinear/linear finite difference schemes [6], and the local boundary integral method [7]. However, these works do not consider the physical properties of the EFK model. Recently, Sun et al. [8,9] proposed two convex splitting variable step BDF2/BDF3 schemes for the EFK model, and the proposed schemes preserve the modified discrete energy dissipation law.
The development of high-precision numerical schemes for PDEs of gradient-flow-type has attracted much attention. In this direction, it is worth mentioning the integrating factor Runge–Kutta (IFRK) method [10]. The IFRK method has demonstrated remarkable advantages for equations with stiff linear terms. The exponential function in this method provides the exact solution of the linear part, so the stiffness does not restrict the step size and the solution to the nonlinear part may be approximated explicitly. As a result, the IFRK method is often explicit [11,12]. In recent years, the IFRK method has attracted widespread attention for solving partial differential equations. Ju et al. [13] proposed the MBP-preserving IFRK method for semilinear parabolic equations. Li et al. [14] proposed a class of unconditional MBP-preserving IFRK schemes for the conservative Allen–Cahn equation. Zhang et al. [15,16] developed a class of high-order structure-preserving IFRK methods for the Allen–Cahn equation. Then, they further [17] proposed and analyzed a series of temporal up to fourth-order unconditionally structure-preserving single-step methods to solve the Allen–Cahn equation, and introduced parametric IFRK (pIFRK) schemes which can be used to construct higher-order parametric single-step methods. To the authors’ best knowledge, there are very few works in the literature on high-precision numerical methods for the EFK model. The objective of this paper is to develop a class of efficient, high-order-accurate schemes for the EFK model based on the explicit IFRK method coupled with non-decreasing abscissas (eIFRK+) [13,18].
The rest of the article is arranged as follows: In Section 2, a class of eIFRK+ schemes is proposed for solving the EFK model, and the corresponding theoretical analysis is given in Section 3. Numerical experiments are presented to test the performance of the proposed numerical schemes in Section 4, and some concluding remarks are given in Section 5.

2. eIFRK+ Fourier-Spectral Schemes for EFK Model

Define a periodic spatial grid h = ( x i , y j ) = a + i h , a + j h , 0 i , j N 1 ) with h = b a N (N is even). All of the 2D periodic grid functions defined on h are denoted by M h .
The discrete Fourier transform ϕ ˜ = Q ϕ and the corresponding inverse transform ϕ = Q 1 ϕ ˜ are defined by [19]
ϕ ˜ p q = h 2 ( b a ) 2 i = 0 N 1 j = 0 N 1 ϕ i j exp i 2 p π ( x i a ) b a i 2 q π ( y j a ) b a , p , q = 0 , ± 1 , ± 2 , , N 2 ,
ϕ i j = p = N / 2 N / 2 q = N / 2 N / 2 ϕ ˜ p q ( t ) exp i 2 p π ( x i a ) b a + i 2 q π ( y j a ) b a , ϕ ˜ N 2 , q = ϕ ˜ N 2 , q , ϕ ˜ p , N 2 = ϕ ˜ p , N 2 ,
respectively. We further define the operators D ^ x and D ^ y on M h ^ = { Q ϕ | ϕ M h } as
( D ^ x ϕ ˜ ) p q = 2 p π i b a ϕ ˜ p q , ( D ^ y ϕ ˜ ) p q = 2 q π i b a ϕ ˜ p q ,
and then, the Laplace operator can be approximated by
Δ N = Q 1 ( D ^ x 2 + D ^ y 2 ) Q .
Next, we give some fully discrete eIFRK+ schemes for solving the EFK equation. Let τ = T M and define the time node as t k = k τ ( k = 0 , 1 , M ). According to the definitions of discrete spatial operators presented in Section 2, the EFK model (1) is reduced to the following nonlinear ordinary differential equations (ODEs) in time
d d t v ( t ) = L N v ( t ) + g ( v ( t ) ) ,
where L N = κ Δ N 2 + Δ N and v ( t ) = { v i j ( t ) } M h . Let v ( t ) = e t L N u ( t ) , then Equation (4) can be transformed into
d d t u ( t ) = e t L N g e t L N u ( t ) ,
which can be solved by the following s-stage classical RK scheme [18]
u ( 0 ) = u k , u ( i ) = j = 0 i 1 [ η i j u ( j ) + τ γ i j e t k ( j ) L N g ( e t k ( j ) L N u ( j ) ) ] , 1 i n , u k + 1 = u ( n ) ,
where t k ( j ) = t k + d j τ , d j is the abscissa at each j-th layer. η i j 0 and j = 0 i 1 η i j = 1 , γ i j is a real constant.
Since v ( j ) = e t k ( j ) L N u ( j ) , then (6) becomes
v ( 0 ) = v k , v ( i ) = j = 0 i 1 e ( d i d j ) τ L N [ η i j v ( j ) + τ γ i j g ( v j ) ] , 1 i n , v k + 1 = v ( n ) ,
where d i d j 0 , d 0 = 0 and d n = 1 . By a simple derivation, the above scheme can be rewritten in the following equivalent form [13]:
v ( 0 ) = v k , v ( i ) = e d i τ L N v k + τ j = 0 i 1 ε i j e ( d i d j ) τ L N g ( v ( j ) ) , 1 i n , v k + 1 = v ( n ) ,
where j = 0 i 1 ε i j 1 , γ i j = ε i j s = j + 1 i 1 η i s ε s j .
Some specific cases of system (7) are presented, which are denoted by eIFRK+ ( n , q ) , where q denotes the time accuracy:
eIFRK + ( 1 , 1 ) : v k + 1 = e τ L N [ v k + τ g ( v k ) ] .
eIFRK + ( 2 , 2 ) : v ( 1 ) = e τ L N [ v k + τ g ( v k ) ] , v k + 1 = 1 2 e τ L N v k + 1 2 [ v ( 1 ) + τ g ( v ( 1 ) ) ] .
eIFRK + ( 3 , 3 ) : v ( 1 ) = 1 2 e 2 τ 3 L N v k + 1 2 e 2 τ 3 L N [ v k + 4 τ 3 g ( v k ) ] , v ( 2 ) = 2 3 e 2 τ 3 L N v k + 1 3 [ v ( 1 ) + 4 τ 3 g ( v ( 1 ) ) ] , v k + 1 = 59 128 e τ L N v k + 15 128 e τ L N [ v k + 4 τ 3 g ( v k ) ] + 27 64 e τ 3 L N [ v ( 2 ) + 4 τ 3 g ( v ( 2 ) ) ] .
eIFRK + ( 4 , 4 ) : v ( 1 ) = e τ 2 L N [ v k + τ 2 g ( v k ) ] , v ( 2 ) = 1 2 e τ 2 L N [ v k τ 2 g ( v k ) ] + 1 2 [ v ( 1 ) + τ g ( v ( 1 ) ) ] , v ( 3 ) = 1 9 e τ L N [ v k τ g ( v k ) ] + 2 9 e τ 2 L N [ v ( 1 ) 3 τ 2 g ( v ( 1 ) ) ] + 2 3 e τ 2 L N [ v ( 2 ) + 3 τ 2 g ( v ( 2 ) ) ] , v k + 1 = 1 3 e τ 2 L N [ v ( 1 ) + τ 2 g ( v ( 1 ) ) ] + 1 3 e τ 2 L N v ( 2 ) + 1 3 [ v ( 3 ) + τ 2 g ( v ( 3 ) ) ] .

3. Discrete Error Estimate

In this section, the discrete error estimation of the scheme (8) is discussed. We first give one lemma that helps to prove the discrete error estimation for the scheme.
Lemma 1 
([20]). If a matrix A is negative semi-definite, then we have e t A 2 1 for t > 0 , where · 2 denotes the standard vector or matrix L 2 -norm.
Theorem 1. 
(Error estimate.) Assume that v 0 H p e r p ( Ω ) = { v H p e r p ( Ω ) a n d v s . i s Ω p e r i o d i c } with p > 4 and the solution to the EFK model belongs to C q + 1 ( 0 , T ; H p e r p ( Ω ) ) . Under the conditions of Lemma 1, then the numerical solutions { V k } M h generated by the eIFRK+ schemes (8) with V 0 = I h v 0 satisfy the following error estimate
V k I h v ( t k ) 2 c ˜ h p + τ q .
where c ˜ is a positive constant depending on τ and h. I h : H p e r p ( Ω ) M h is a sample operator as I h ( v ( x , y ) ) = { v ( x i , y j ) } .
Proof. 
Let Z ( t ) = V ( t ) I h v ( x , y , t ) with Z ( 0 ) = 0 ; the difference between Equations (1) and (4) leads to
d d t Z ( t ) = L N Z ( t ) + g ( V ( t ) ) g ( I h ( v ) ) + R N ( v ) ,
where the truncated error R N ( v ) = L ( I h v ) L N ( I h v ) satisfies R N 2 c h p .
Multiplying both sides of (14) by e t L N , then we have
d d t e t L N Z ( t ) = e t L N g ( V ( t ) ) g ( I h ( v ) ) + R N ( v ) .
Using the eIFRK+ scheme (8), we have the following fully discrete scheme for the above equation
Z ( i ) = e d i τ L N Z k + τ j = 0 i 1 ε i j e ( d i d j ) τ L N g ( V j ) g ( I h v ( j ) ) + R N ( v ( j ) ) , 1 i n 1 , Z k + 1 = e τ L N Z k + τ i = 0 n 1 ρ i e ( 1 d i ) τ L N g ( V j ) g ( I h v ( j ) ) + R N ( v ( j ) ) + ξ k ,
where j = 0 i 1 ε i j 1 , i = 0 n 1 ε n i = i = 0 n 1 ρ i = 1 , ξ k is the truncation error satisfying ξ k 2 c τ q + 1 .
By the Lipschitz condition for g ( v ) , we have
g ( V j ) g ( I h v ( j ) ) 2 L V j I h v ( j ) 2 = L Z ( j ) 2 .
Combination of the above estimates leads to
Z ( i ) 2 e d i τ L N Z k 2 + τ j = 0 i 1 ε i j e ( d i d j ) τ L N ( g ( V j ) g ( I h v ( j ) ) ) 2 + τ j = 0 i 1 ε i j e ( d i d j ) τ L N R N ( v j ) 2 Z k 2 + L τ j = 0 i 1 Z j 2 + c h p τ .
By using the Gronwall’s inequality, we further have
Z ( i ) 2 e L ( i + 1 ) τ ( Z k 2 + c h p τ ) .
Similar to the estimation of Z k + 1 2 , we have
Z k + 1 2 Z k 2 + L τ i = 0 n 1 Z ( i ) 2 + c h p τ + ξ k 2 .
Combining (17) with (18), we have
Z k + 1 2 Z k 2 + L τ i = 0 n 1 e L ( i + 1 ) τ ( Z k 2 + c h p τ ) + c h p τ + ξ k 2 Z k 2 + L τ n e L n τ Z k 2 + c h p τ + c τ q + 1 .
Summing the above inequality from k = 0 to m and using Gronwall’s inequality, we obtain the desired result
Z m + 1 | 2 e L n T e L T ( c h p T + c τ q T ) c ˜ ( h p + τ q ) .

4. Numerical Experiments

We examine the accuracy and efficiency of the proposed method through three numerical experiments, including convergence and diminishing energy.
Example 1. 
The purpose of the first numerical experiment is to verify the role of parameter κ. Consider the 1D EFK model in Ω = [ 4 , 4 ] with the initial condition
v 0 x = sin π x .
We choose N = 256 and τ = 0.01 . Figure 1 shows the numerical results using scheme eIFRK+(4,4) at different times t = 0 , 0.05 , 0.1 , 0.15 , 0.2 . It is observed that the dynamic evolution of κ = 0 is almost identical to that of κ = 10 4 . However, when κ = 0.1 , the solution rapidly evolves from the initial state to zeros, which illustrates that the coefficient κ in the EFK model is a stable parameter.
Example 2. 
To show the temporal accuracy and convergence of the proposed schemes, we consider the 2D EFK model in Ω = [ 16 , 16 ] 2 with κ = 0.01 and the initial condition
v 0 ( x , y ) = 0.1 0.2 cos π ( x 12 ) 16 sin π ( y 12 ) 16 + 0.1 cos 2 π ( x + 10 ) 32 sin 2 π ( y + 3 ) 32 0.2 sin 2 π x 8 cos π ( y 6 ) 8 .
We adopt the uniform 1024 × 1024 spatial mesh, which is sufficiently fine so that the errors caused by the spatial approximation can be ignored. The error between two different time steps τ and τ 2 is calculated. The computational results are presented in Table 1. One may see that the numerical orders of time accuracy are close to the optimal order.
Example 3. 
Consider the EFK model in Ω = [ 1 , 1 ] 2 with the initial condition
v 0 ( x , y ) = 1 + tanh 0.4 x 2 + y 2 ε 2 tanh 0.3 x 2 + y 2 ε 2 .
We choose ε = 5 h 2 t a n h 1 ( 0.9 ) , N = 256 , τ = 0.01 and different κ, Figure 2 shows the evolutions of the snapshots and the energies of the numerical solutions obtained by the eIFRK+(4,4) scheme. It can be seen that the evolution speed of the numerical solution is largely clipped with the increase in κ, and the energy obtained from the energy image gradually reaches a steady state with decreasing time.

5. Conclusions

Based on the explicit IFRK method coupled with nondecreasing abscissas, we have obtained a class of fast and effective numerical schemes for the EFK model. The optimal error estimates of the fully discrete schemes have been analyzed. Three numerical experiments were carried out to test the accuracy and applicability of the proposed schemes.

Author Contributions

Y.W.: Conceptualization, Methodology, Software, Writing—original draft. S.Z.: Methodology, Writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data will be made available on reasonable request.

Acknowledgments

The authors would like to thank the editor and referees for their valuable comments and suggestions which helped us to improve the results of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fisher, R.A. The wave of advance of advantageous genes. Ann. Eugen. 1937, 7, 355–369. [Google Scholar] [CrossRef]
  2. Kolmogorov, A.N.; Petrovskii, I.G.; Piskunov, N.S. A study of the diffusion equation with increase in the quantity of matter, and its application to a biological problem. Mosc. Univ. Math. Bull. 1937, 1, 1–25. [Google Scholar]
  3. Coullet, P.; Elphick, C.; Repaux, D. Nature of spatial chaos. Phys. Rev. Lett. 1987, 58, 431. [Google Scholar] [CrossRef]
  4. Saarloos, W.V. Dynamical velocity selection: Marginal stability. Phys. Rev. Lett. 1987, 58, 2571. [Google Scholar] [CrossRef]
  5. Danumjaya, P.; Pani, A.K. Orthogonal cubic spline collocation method for the extended Fisher–Kolmogorov equation. J. Comput. Appl. Math. 2005, 174, 101–117. [Google Scholar] [CrossRef]
  6. Li, S.G.; Xu, D.; Zhang, J.; Sun, C.J. A new three-level fourth-order compact finite difference scheme for the extended Fisher–Kolmogorov equation. Appl. Numer. Math. 2022, 178, 41–51. [Google Scholar] [CrossRef]
  7. Ilati, M.; Dehghan, M. Direct local boundary integral equation method for numerical solution of extended Fisher–Kolmogorov equation. Eng. Comput. 2018, 34, 203–213. [Google Scholar] [CrossRef]
  8. Sun, Q.H.; Ji, B.Q.; Zhang, L. A convex splitting BDF2 method with variable time-steps for the extended Fisher–Kolmogorov equation. Comput. Math. Appl. 2022, 114, 73–82. [Google Scholar] [CrossRef]
  9. Sun, Q.H.; Wang, J.D.; Zhang, L.M. A new third-order energy stable technique and error estimate for the extended Fisher–Kolmogorov equation. Comput. Math. Appl. 2023, 142, 198–207. [Google Scholar] [CrossRef]
  10. Isherwood, L.; Grant, Z.J.; Gottlieb, S. Strong stability preserving integrating factor Runge–Kutta methods. SIAM J. Numer. Anal. 2018, 58, 3276–3307. [Google Scholar] [CrossRef]
  11. Zhai, S.Y.; Weng, Z.F.; Yang, Y.F. A high order operator splitting method based on spectral deferred correction for the nonlocal viscous Cahn–Hilliard equation. J. Comput. Phys. 2021, 446, 110632. [Google Scholar] [CrossRef]
  12. Zhai, S.Y.; Weng, Z.F.; Feng, X.L.; He, Y.N. Stability and error estimate of the operator splitting method for the phase field crystal equation. J. Sci. Comput. 2021, 86, 8. [Google Scholar] [CrossRef]
  13. Ju, L.L.; Li, X.; Qiao, Z.H.; Yang, J. Maximum bound principle preserving integrating factor Runge–Kutta methods for semilinear parabolic equations. J. Comput. Phys. 2021, 439, 110405. [Google Scholar] [CrossRef]
  14. Li, J.W.; Li, X.; Ju, L.L.; Feng, X.L. Stabilized integrating factor Runge–Kutta method and unconditional preservation of maximum bound principle. SIAM J. Sci. Comput. 2021, 43, A1780–A1802. [Google Scholar] [CrossRef]
  15. Zhang, H.; Yan, J.Y.; Qian, X.; Chen, X.W.; Song, S.H. Explicit third-order unconditionally structure-preserving schemes for conservative Allen–Cahn equations. J. Sci. Comput. 2022, 90, 1–29. [Google Scholar] [CrossRef]
  16. Zhang, H.; Yan, J.Y.; Qian, X.; Song, S.H. Numerical analysis and applications of explicit high order maximum principle preserving integrating factor Runge–Kutta schemes for Allen–Cahn equation. Appl. Numer. Math. 2021, 161, 372–390. [Google Scholar] [CrossRef]
  17. Zhang, H.; Yan, J.Y.; Qian, X.; Song, S.H. Up to fourth-order unconditionally structure-preserving parametric single-step methods for semilinear parabolic equations. Comput. Methods Appl. Mech. Eng. 2022, 393, 114817. [Google Scholar] [CrossRef]
  18. Shu, C.W.; Osher, S. Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys. 1988, 77, 439–471. [Google Scholar] [CrossRef]
  19. Shen, J.; Tang, T.; Wang, L.L. Spectral Methods: Algorithms, Analysis and Applications; Springer Science: London, UK; New York, NY, USA, 2011; Volume 41, pp. 1–45. [Google Scholar]
  20. Nan, C.X.; Song, H.L. The high-order maximum-principle-preserving integrating factor Runge–Kutta methods for nonlocal Allen–Cahn equation. J. Comput. Phys. 2022, 456, 111028. [Google Scholar] [CrossRef]
Figure 1. Numerical solutions with different κ values for Example 1. (a) κ = 0 ; (b) κ = 10 4 ; (c) κ = 0.1 .
Figure 1. Numerical solutions with different κ values for Example 1. (a) κ = 0 ; (b) κ = 10 4 ; (c) κ = 0.1 .
Mca 28 00110 g001
Figure 2. The numerical solution u at different times for Example 3. Top: κ = 10 4 . Bottom: κ = 0.004 . (a) T = 0 ; (b) T = 0.004 ; (c) T = 0.012 ; (d) T = 0.02 ; (e) T = 0.028 ; (f) E ( v ) ; (g) T = 0 ; (h) T = 0.004 ; (i) T = 0.012 ; (j) T = 0.016 ; (k) T = 0.02 ; (l) E ( v ) .
Figure 2. The numerical solution u at different times for Example 3. Top: κ = 10 4 . Bottom: κ = 0.004 . (a) T = 0 ; (b) T = 0.004 ; (c) T = 0.012 ; (d) T = 0.02 ; (e) T = 0.028 ; (f) E ( v ) ; (g) T = 0 ; (h) T = 0.004 ; (i) T = 0.012 ; (j) T = 0.016 ; (k) T = 0.02 ; (l) E ( v ) .
Mca 28 00110 g002
Table 1. Temporal errors and convergence orders at T = 1 with κ = 0.01 and N = 1024 for Example 2.
Table 1. Temporal errors and convergence orders at T = 1 with κ = 0.01 and N = 1024 for Example 2.
eIFRK+MErrmaxRateErr L 2 Rate
( 1 , 1 ) 16 3.9212 × 10 3 -1.9491 × 10 2 -
32 1.9980 × 10 3 0.97289.85420.9840
64 1.0086 × 10 3 0.98624.9540 × 10 3 0.9921
128 5.0671 × 10 4 0.9931 2.4837 × 10 3 0.9961
( 2 , 2 ) 16 1.2908 × 10 4 - 6.5886 × 10 4 -
32 3.3107 × 10 5 1.9630 1.6909 × 10 4 1.9622
64 8.3832 × 10 6 1.9816 4.2826 × 10 5 1.9812
128 2.1092 × 10 6 1.9908 1.0776 × 10 5 1.9906
( 3 , 3 ) 16 1.4809 × 10 6 - 5.0559 × 10 6 -
32 1.8992 × 10 7 2.9631 6.5053 × 10 7 2.9583
64 2.4046 × 10 8 2.9815 8.2502 × 10 8 2.9791
128 3.0251 × 10 9 2.9907 1.0388 × 10 8 2.9896
( 4 , 4 ) 16 3.1132 × 10 8 - 1.6382 × 10 7 -
32 1.9696 × 10 9 3.9824 1.0359 × 10 8 3.9830
64 1.2395 × 10 10 3.9901 6.5125 × 10 10 3.9915
128 7.7799 × 10 12 3.9939 4.0823 × 10 11 3.9958
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, Y.; Zhai, S. Explicit Integrating Factor Runge–Kutta Method for the Extended Fisher–Kolmogorov Equation. Math. Comput. Appl. 2023, 28, 110. https://doi.org/10.3390/mca28060110

AMA Style

Wang Y, Zhai S. Explicit Integrating Factor Runge–Kutta Method for the Extended Fisher–Kolmogorov Equation. Mathematical and Computational Applications. 2023; 28(6):110. https://doi.org/10.3390/mca28060110

Chicago/Turabian Style

Wang, Yanan, and Shuying Zhai. 2023. "Explicit Integrating Factor Runge–Kutta Method for the Extended Fisher–Kolmogorov Equation" Mathematical and Computational Applications 28, no. 6: 110. https://doi.org/10.3390/mca28060110

APA Style

Wang, Y., & Zhai, S. (2023). Explicit Integrating Factor Runge–Kutta Method for the Extended Fisher–Kolmogorov Equation. Mathematical and Computational Applications, 28(6), 110. https://doi.org/10.3390/mca28060110

Article Metrics

Back to TopTop