
Fractal ordinary differential equations are successfully established by He's fractal derivative in a fractal space, and their variational principles are obtained by semi-inverse transform method.Taylor series method is used to solve the given fractal equations with initial boundary value conditions, and sometimes Ying Buzu algorithm play an important role in this process. Examples show the Taylor series method and Ying Buzu algorithm are powerful and simple tools.
Citation: Yi Tian. Approximate solution of initial boundary value problems for ordinary differential equations with fractal derivative[J]. Mathematical Modelling and Control, 2022, 2(2): 75-80. doi: 10.3934/mmc.2022009
[1] | Ji-Huan He, Shuai-Jia Kou, Hamid M. Sedighi . An ancient Chinese algorithm for two-point boundary problems and its application to the Michaelis-Menten kinetics. Mathematical Modelling and Control, 2021, 1(4): 172-176. doi: 10.3934/mmc.2021016 |
[2] | Abduljawad Anwar, Shayma Adil Murad . On the Ulam stability and existence of Lp-solutions for fractional differential and integro-differential equations with Caputo-Hadamard derivative. Mathematical Modelling and Control, 2024, 4(4): 439-458. doi: 10.3934/mmc.2024035 |
[3] | Ihtisham Ul Haq, Nigar Ali, Hijaz Ahmad . Analysis of a chaotic system using fractal-fractional derivatives with exponential decay type kernels. Mathematical Modelling and Control, 2022, 2(4): 185-199. doi: 10.3934/mmc.2022019 |
[4] | Min Zeng, Yongxin Yuan . On the solutions of the dual matrix equation A⊤XA=B. Mathematical Modelling and Control, 2023, 3(3): 210-217. doi: 10.3934/mmc.2023018 |
[5] | Mrutyunjaya Sahoo, Dhabaleswar Mohapatra, S. Chakraverty . Wave solution for time fractional geophysical KdV equation in uncertain environment. Mathematical Modelling and Control, 2025, 5(1): 61-72. doi: 10.3934/mmc.2025005 |
[6] | Wenxv Ding, Ying Li, Dong Wang, AnLi Wei . Constrainted least squares solution of Sylvester equation. Mathematical Modelling and Control, 2021, 1(2): 112-120. doi: 10.3934/mmc.2021009 |
[7] | Ihteram Ali, Imtiaz Ahmad . Applications of the nonlinear Klein/Sinh-Gordon equations in modern physics: a numerical study. Mathematical Modelling and Control, 2024, 4(3): 361-373. doi: 10.3934/mmc.2024029 |
[8] | C. Kavitha, A. Gowrisankar . Fractional integral approach on nonlinear fractal function and its application. Mathematical Modelling and Control, 2024, 4(3): 230-245. doi: 10.3934/mmc.2024019 |
[9] | Rashid Jan, Normy Norfiza Abdul Razak, Sania Qureshi, Imtiaz Ahmad, Salma Bahramand . Modeling Rift Valley fever transmission: insights from fractal-fractional dynamics with the Caputo derivative. Mathematical Modelling and Control, 2024, 4(2): 163-177. doi: 10.3934/mmc.2024015 |
[10] | Ruxin Zhang, Zhe Yin, Ailing Zhu . Numerical simulations of a mixed finite element method for damped plate vibration problems. Mathematical Modelling and Control, 2023, 3(1): 7-22. doi: 10.3934/mmc.2023002 |
Fractal ordinary differential equations are successfully established by He's fractal derivative in a fractal space, and their variational principles are obtained by semi-inverse transform method.Taylor series method is used to solve the given fractal equations with initial boundary value conditions, and sometimes Ying Buzu algorithm play an important role in this process. Examples show the Taylor series method and Ying Buzu algorithm are powerful and simple tools.
The ordinary differential equations (ODEs) arise in many fields like the physics, mechanics, economics and management, etc. There are many methods for solving nonlinear ODEs, for example, the He's frequency formulation [1,2], homotopy perturbation method [3,4], variational iteration method [5,6], and variational-based method [7,8]. In addition, there have been some research results on differential equations recently [9,10,11,12,13,14,15]. Each method has its advantages and disadvantages. For example, the Exp-function method can lead to the analytical solutions, but its complex calculation makes those inaccessible who are not familiar some mathematics software. The variational-based methods can obtain a globally valid solution, however, it is extremely difficult to establish a needed variational principle for a complex nonlinear problem, etc.
In my memory, most researchers in the wold pay little attention to Chinese mathematics, especially on ancient Chinese mathematics. So the present author feels strongly necessary to do some work let the world know that China has 5000 years of civilization, not only in social sciences, but also in natural sciences. This paper concerns briefly a famous ancient Chinese algorithm, named Ying Buzu. Every student knows Newton's iteration method from a textbook, which is widely used in numerical simulation, what few may know is that its ancient Chinese partner, Ying Buzu algorithm, in about second century BC has much advantages over Newton's method [16].
In this paper, we mainly study using Taylor series method [17,18,19,20,21,22] and Ying Buzu algorithm [23] to solve ordinary differential equations with initial boundary value conditions in a fractal space. To the best of our knowledge, the hybrid method has not been studied.
The remainder of this paper is organized as follows. In section 2, we briefly study the initial value problems for ODEs. In section 3, boundary value problems for ODEs with Neumann boundary conditions are studied. In section 4, boundary value problems for ODEs with Dirichlet boundary conditions are considered. In section 5, we conclude this paper and some discussions are also given there.
Consider the initial value problem [24]
y″+y=2x−1,y(1)=1,y′(1)=3. | (2.1) |
In a fractal space, Eq.(2.1) can be described as follows
ddxα(dydxα)+y=2x−1,y(1)=1,dydxα(1)=3, | (2.2) |
where dydxα is the fractal derivative [25,26]
dydxα(x0)=Γ(1+α)limx−x0→ΔxΔx≠0y(x)−y(x0)x−x0. | (2.3) |
The variational principle of Eq (2.2) can be given by semi-inverse method as follows
J=∫{−12(dydxα)2+12y2−2xy+y}dxα. | (2.4) |
Using the two-scale transform [27,28]
s=xα, | (2.5) |
we can convert Eq (2.2) approximately into the following one
d2yds2+y=2s−1,y(1)=1,dyds(1)=3. | (2.6) |
From Eq.(2.6), we have
y″(1)=0,y(3)(1)=−1,y(4)(1)=0. | (2.7) |
The Taylor series solution of Eq (2.6) is
y(s)=y(1)+(s−1)y′(1)+12(s−1)2y″(1) |
+16(s−1)3y(3)(1)+124(s−1)4y(4)(1) |
=1+3(s−1)−16(s−1)3, | (2.8) |
which converges to the exact solution, which is
y(s)=2s−1+sin(s−1). | (2.9) |
Figure 1 shows the Taylor series solutions, which is very close to the exact solutions.
Consider the boundary value problem [24]
d2θdx2=ϵθ4,θ′(0)=0,θ(1)=1. | (3.1) |
In a fractal space, Eq (3.1) can be described as follows
ddxα(dθdxα)=ϵθ4,dθdxα(0)=0,θ(1)=1, | (3.2) |
where dθdxα is the fractal derivative [25,26]
dθdxα(x0)=Γ(1+α)limx−x0→ΔxΔx≠0θ(x)−θ(x0)x−x0. | (3.3) |
The variational principle of Eq (3.2) can be given by semi-inverse method as follows
J=∫{−12(dθdxα)2−15ϵθ5}dxα. | (3.4) |
Using the two-scale transform
s=xα, | (3.5) |
Equation (3.2) can be converted approximately into the following one
d2θds2=ϵθ4,dθds(0)=0,θ(1)=1. | (3.6) |
For simplicity, we assume
ϵ=1,θ(0)=c, | (3.7) |
from (3.6), we have
θ″(0)=c4, | (3.8) |
θ(3)(0)=4c3, | (3.9) |
θ(4)(0)=12c2. | (3.10) |
The 4th-order Taylor series solution is
θ(s)=θ(0)+sθ′(0)+12s2θ″(0)+16s3θ(3)(0)+124s4θ(4)(0) |
=c+12c4s2+23c3s3+12c2s4, | (3.11) |
incorporating the boundary condition θ(1)=1, we have
c+12c4+23c3+12c2=1. | (3.12) |
We use the Ying Buzu algorithm [23] to solve c and write (3.12) in the form
R(c)=c+12c4+23c3+12c2−1. | (3.13) |
Assume the two initial solutions are
θ1(0)=0.6,θ2(0)=0.8, | (3.14) |
we obtain the following residuals
R1(0.6)=−0.0112,R2(0.8)=0.666133, | (3.15) |
the inital guess θ(0) can be updated as
θ3(0)=R2θ1−R1θ2R2−R1=0.603307, | (3.16) |
the shooting process using (3.16) results in
θ3(1)=0.997931, | (3.17) |
which deviates the exact value of θ(1)=1 with a relative error of 0.2%.
The 4th-order Taylor series solution is
θ(s)=0.603307+0.066241s2+0.146394s3+0.18199s4. | (3.18) |
Figure 2 shows the Taylor series solutions, which approximately meet the requirement of the boundary condition.
Consider the boundary value problem [29]
u″(t)+tu(t)=t3+2,u(0)=0,u(1)=1. | (4.1) |
In a fractal space, Equation (4.1) can be described as follows
ddtα(dudtα)+tu(t)=t3+2,u(0)=0,u(1)=1, | (4.2) |
where dudtα is the fractal derivative [25,26]
dudtα(t0)=Γ(1+α)limt−t0→ΔtΔt≠0u(t)−u(t0)t−t0. | (4.3) |
The variational principle of Eq (4.2) can be given by semi-inverse method as follows
J=∫{−12(dudtα)2+12tu(t)2−(t3+2)u(t)}dtα. | (4.4) |
Using the two-scale transform
s=tα, | (4.5) |
Equation (4.2) can be converted approximately into the following one
d2uds2+su(s)=s3+2,u(0)=0,u(1)=1. | (4.6) |
We assume
u′(0)=α, | (4.7) |
from (4.6), we have
u″(0)=2, | (4.8) |
u(3)(0)=0, | (4.9) |
u(4)(0)=−2α. | (4.10) |
The 4th-order Taylor series solution is
u(s)=u(0)+su′(0)+12s2u″(0)+16s3u(3)(0)+124s4u(4)(0) |
=αs+s2−112αs4. | (4.11) |
by using the boundary condition u(1)=1, we have
α+1−112α=1, | (4.12) |
and α=0. The 4th-order Taylor series solution is
u(s)=s2. | (4.13) |
which is also the exact solution of Eq (4.1).
Fractal ordinary differential equations are successfully established by He's fractal derivative in a fractal space, and their variational principles are obtained by semi-inverse transform method.The two-scale transform method and Taylor series method are adopted to solve the fractal ODEs with initial boundary value conditions.The examples show the Taylor series method is simple and effective and the ancient Chinese Ying Buzu algorithm is a simple and straightforward tool to two-point boundary value problems. In the future, we will study how to extend this hybrid method to PDEs.
The author has no conflicts of interest to declare that are relevant to the content of this article.
[1] |
J. He, W. Hou, N. Qie, K. Gepreel, A. Shirazi, H. Mohammad-Sedighi, Hamiltonian-based frequency-amplitude formulation for nonlinear oscillators, Facta Univ.-Ser. Mech., 19 (2021), 199–208. https://doi.org/10.22190/FUME201205002H doi: 10.22190/FUME201205002H
![]() |
[2] |
J. He, Q. Yang, C. He, Y. Khan, A simple frequency formulation for the tangent oscillator, Axioms, 10 (2021), 320. https://doi.org/10.3390/axioms10040320 doi: 10.3390/axioms10040320
![]() |
[3] |
N. Anjum, J. He, Q. Ain, D. Tian, Li-He's modified homotopy perturbation method for doubly-clamped electrically actuated microbeams-based microelectromechanical system, Facta Univ.-Ser. Mech., 4 (2021), 601–612. https://doi.org/10.22190/FUME210112025A doi: 10.22190/FUME210112025A
![]() |
[4] |
J. He, Y. El-Dib, The enhanced homotopy perturbation method for axial vibration of strings, Facta Univ.-Ser. Mech., 4 (2021), 735–750. https://doi.org/10.22190/FUME210125033H doi: 10.22190/FUME210125033H
![]() |
[5] |
M. Nadeem, J. He, He-Laplace variational iteration method for solving the nonlinear equations arising in chemical kinetics and population dynamics, J. Math. Chem., 59 (2021), 1234–1245. https://doi.org/10.1007/s10910-021-01236-4 doi: 10.1007/s10910-021-01236-4
![]() |
[6] |
M. Nadeem, J. He, Analysis of nonlinear vibration of nano/microelectromechanical system switch induced by electromagnetic force under zero initial conditions, Alex. Eng. J., 59 (2020), 4343–4352. https://doi.org/10.1016/j.aej.2020.07.039 doi: 10.1016/j.aej.2020.07.039
![]() |
[7] |
J. He, Lagrange crisis and generalized variational principle for 3D unsteady flow, International Journal of Numerical Methods for Heat and Fluid Flow, 30 (2020), 1189–1196. https://doi.org/10.1108/HFF-07-2019-0577 doi: 10.1108/HFF-07-2019-0577
![]() |
[8] |
J. He, C. Sun, A variational principle for a thin film equation, J. Math. Chem., 57 (2019), 2075–2081. https://doi.org/10.1007/s10910-019-01063-8 doi: 10.1007/s10910-019-01063-8
![]() |
[9] |
E. Bonyah, M. Yavuz, D. Baleanu, S. Kumar, A robust study on the listeriosis disease by adopting fractal-fractional operators, Alex. Eng. J., 61 (2022), 2016–2028. https://doi.org/10.1016/j.aej.2021.07.010 doi: 10.1016/j.aej.2021.07.010
![]() |
[10] |
P. Naik, Z. Eskandari, H. Shahkari, Flip and generalized flip bifurcations of a two-dimensional discrete-time chemical model, Mathematical Modelling and Numerical Simulation with Applications, 1 (2021), 95–101. https://doi.org/10.53391/mmnsa.2021.01.009 doi: 10.53391/mmnsa.2021.01.009
![]() |
[11] |
Z. Hammouch, M. Yavuz, N. Ozdemir, Numerical solutions and synchronization of a variable-order fractional chaotic system, Mathematical Modelling and Numerical Simulation with Applications, 1 (2021), 11–23. https://doi.org/10.53391/mmnsa.2021.01.002 doi: 10.53391/mmnsa.2021.01.002
![]() |
[12] |
R.Jena, S. Chakraverty, M.Yavuz, Two-Hybrid techniques coupled with an integral transformation for caputo time-fractional navier-stokes equations, Progress in Fractional Differentiation and Applications, 6 (2020), 201–213. https://doi.org/10.18576/pfda/060304 doi: 10.18576/pfda/060304
![]() |
[13] |
P. Kumar, V. Erturk, Dynamics of cholera disease by using two recentfractional numerical methods, Mathematical Modelling and Numerical Simulation with Applications, 1 (2021), 102–111. https://doi.org/10.53391/mmnsa.2021.01.010 doi: 10.53391/mmnsa.2021.01.010
![]() |
[14] |
B. Daşbaşl, Stability analysis of an incommensurate fractional-order SIR model, Mathematical Modelling and Numerical Simulation with Applications, 1 (2021), 44–55. https://doi.org/10.53391/mmnsa.2021.01.005 doi: 10.53391/mmnsa.2021.01.005
![]() |
[15] |
F. Ozkose, M. Senel, R. Habbireeh, Fractional-order mathematical modelling of cancercells-cancer stem cells-immune system interaction with chemotherapy, Mathematical Modelling and Numerical Simulation with Applications, 1 (2021), 67–83. https://doi.org/10.53391/mmnsa.2021.01.007 doi: 10.53391/mmnsa.2021.01.007
![]() |
[16] |
C. He, An introduction to an ancient Chinese algorithm and its modification, International Journal of Numerical Methods for Heat and Fluid Flow, 26 (2016), 2486–2491. https://doi.org/10.1108/HFF-09-2015-0377 doi: 10.1108/HFF-09-2015-0377
![]() |
[17] |
J. He, A simple approach to one-dimensional convection-diffusion equation and its fractional modification for E reaction arising in rotating disk electrodes, Electroanal. Chem., 854 (2019), 113565. https://doi.org/10.1016/j.jelechem.2019.113565 doi: 10.1016/j.jelechem.2019.113565
![]() |
[18] |
C. He, Y. Shen, F. Ji, J. He, Taylor series solution for fractal bratu-type equation arising in electrospinning process, Fractals, 28 (2020), 2050011. https://doi.org/10.1142/S0218348X20500115 doi: 10.1142/S0218348X20500115
![]() |
[19] |
J. He, A short review on analytical methods for a fully fourth-order nonlinear integral boundary value problem with fractal derivatives, International Journal of Numerical Methods for Heat and Fluid Flow, 30 (2020), 4933–4943. https://doi.org/10.1108/HFF-01-2020-0060 doi: 10.1108/HFF-01-2020-0060
![]() |
[20] |
J. He, F. Ji, Taylor series solution for Lane-Emden equation, J. Math. Chem., 57 (2019), 1932–1934. https://doi.org/10.1007/s10910-019-01048-7 doi: 10.1007/s10910-019-01048-7
![]() |
[21] | J. He, X. Jin, A short review on analytical methods for the capillary oscillator in a nanoscale deformable tube, Math. Method. Appl. Sci., (2020). https://doi.org/10.1002/mma.6321 |
[22] |
J. He, Taylor series solution for a third order boundary value problem arising in Architectural Engineering, Ain Shams Eng. J., 11 (2020), 1411–1414. https://doi.org/10.1016/j.asej.2020.01.016 doi: 10.1016/j.asej.2020.01.016
![]() |
[23] |
J. He, S. Kou, H. Sedighi, An ancient Chinese algorithm for two-point boundary problems and its application to the Michaelis-Menten kinetics, Mathematical Modelling and Control, 1 (2021), 172–176. https://doi.org/10.3934/mmc.2021016 doi: 10.3934/mmc.2021016
![]() |
[24] | R. Agarwal, D. O'Regan, Ordinary and partial differential equations: with special functions, fourier series, and boundary value problems, 1 Eds., New York: Springer-Verlag, 2009. https://doi.org/10.1007/978-0-387-79146-3 |
[25] |
J. He, Fractal calculus and its geometrical explanation, Results Phys., 10 (2018), 272–276. https://doi.org/10.1016/j.rinp.2018.06.011 doi: 10.1016/j.rinp.2018.06.011
![]() |
[26] |
J. He, A tutorial review on fractal spacetime and fractional calculus, Int. J. Theor. Phys., 53 (2014), 3698–3718. https://doi.org/10.1007/s10773-014-2123-8 doi: 10.1007/s10773-014-2123-8
![]() |
[27] |
T. Ain, J. He, On two-scale dimension and its applications, Therm. Sci., 23 (2019), 1707–1712. https://doi.org/10.2298/TSCI190408138A doi: 10.2298/TSCI190408138A
![]() |
[28] |
J. He, F. Ji, Two-scale mathematics and fractional calculus for thermodynamics, Therm. Sci., 23 (2019), 2131–2133. https://doi.org/10.2298/TSCI1904131H doi: 10.2298/TSCI1904131H
![]() |
[29] |
E. Hacloglu, F. Gursoy, S. Maldar, Y. Atalan, G. Milovanovic, Iterative approximation of fixed points and applications to two-point second-order boundary value problems and to machine learning, Appl. Numer. Math., 167 (2021), 143–172. https://doi.org/10.1016/j.apnum.2021.04.020 doi: 10.1016/j.apnum.2021.04.020
![]() |