In this study,we aimto investigate certain triple integral transformand its application to a class of partial differentialequations.We discuss various properties of the new transformincluding inversion, linearity, exi...In this study,we aimto investigate certain triple integral transformand its application to a class of partial differentialequations.We discuss various properties of the new transformincluding inversion, linearity, existence, scaling andshifting, etc. Then,we derive several results enfolding partial derivatives and establish amulti-convolution theorem.Further, we apply the aforementioned transform to some classical functions and many types of partial differentialequations involving heat equations,wave equations, Laplace equations, and Poisson equations aswell.Moreover,wedraw some figures to illustrate 3-D contour plots for exact solutions of some selected examples involving differentvalues in their variables.展开更多
Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational...Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.展开更多
This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)metho...This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.展开更多
Machine learning of partial differential equations(PDEs)from data is a potential breakthrough for addressing the lack of physical equations in complex dynamic systems.Recently,sparse regression has emerged as an attra...Machine learning of partial differential equations(PDEs)from data is a potential breakthrough for addressing the lack of physical equations in complex dynamic systems.Recently,sparse regression has emerged as an attractive approach.However,noise presents the biggest challenge in sparse regression for identifying equations,as it relies on local derivative evaluations of noisy data.This study proposes a simple and general approach that significantly improves noise robustness by projecting the evaluated time derivative and partial differential term into a subspace with less noise.This method enables accurate reconstruction of PDEs involving high-order derivatives,even from data with considerable noise.Additionally,we discuss and compare the effects of the proposed method based on Fourier subspace and POD(proper orthogonal decomposition)subspace.Generally,the latter yields better results since it preserves the maximum amount of information.展开更多
We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural netw...We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural network frameworks.Concretely,we first perform Lagrange interpolation in front of the deep feedforward neural network.The Lagrange basis function has a neat structure and a strong expression ability,which is suitable to be a preprocessing tool for pre-fitting and feature extraction.Second,we introduce small sample learning into training,which is beneficial to guide themodel to be corrected quickly.Taking advantages of the theoretical support of traditional numerical method and the efficient allocation of modern machine learning,LaNets achieve higher predictive accuracy compared to the state-of-the-artwork.The stability and accuracy of the proposed algorithmare demonstrated through a series of classical numerical examples,including one-dimensional Burgers equation,onedimensional carburizing diffusion equations,two-dimensional Helmholtz equation and two-dimensional Burgers equation.Experimental results validate the robustness,effectiveness and flexibility of the proposed algorithm.展开更多
The nonlinearity inmany problems occurs because of the complexity of the given physical phenomena.The present paper investigates the non-linear fractional partial differential equations’solutions using the Caputo ope...The nonlinearity inmany problems occurs because of the complexity of the given physical phenomena.The present paper investigates the non-linear fractional partial differential equations’solutions using the Caputo operator with Laplace residual power seriesmethod.It is found that the present technique has a direct and simple implementation to solve the targeted problems.The comparison of the obtained solutions has been done with actual solutions to the problems.The fractional-order solutions are presented and considered to be the focal point of this research article.The results of the proposed technique are highly accurate and provide useful information about the actual dynamics of each problem.Because of the simple implementation,the present technique can be extended to solve other important fractional order problems.展开更多
The Laplace transformation is a very important integral transform,and it is extensively used in solving ordinary differential equations,partial differential equations,and several types of integro-differential equation...The Laplace transformation is a very important integral transform,and it is extensively used in solving ordinary differential equations,partial differential equations,and several types of integro-differential equations.Our purpose in this study is to introduce the notion of fuzzy double Laplace transform,fuzzy conformable double Laplace transform(FCDLT).We discuss some basic properties of FCDLT.We obtain the solutions of fuzzy partial differential equations(both one-dimensional and two-dimensional cases)through the double Laplace approach.We demonstrate through numerical examples that our proposed method is very successful and convenient for resolving partial differential equations.展开更多
Due to the fact that the fourth-order partial differential equation (PDE) for noise removal can provide a good trade-off between noise removal and edge preservation and avoid blocky effects often caused by the secon...Due to the fact that the fourth-order partial differential equation (PDE) for noise removal can provide a good trade-off between noise removal and edge preservation and avoid blocky effects often caused by the second-order PDE, a domain-based fourth-order PDE method for noise removal is proposed. First, the proposed method segments the image domain into two domains, a speckle domain and a non-speckle domain, based on the statistical properties of isolated speckles in the Laplacian domain. Then, depending on the domain type, different conductance coefficients in the proposed fourth-order PDE are adopted. Moreover, the frequency approach is used to determine the optimum iteration stopping time. Compared with the existing fourth-order PDEs, the proposed fourth-order PDE can remove isolated speckles and keeps the edges from being blurred. Experimental results show the effectiveness of the proposed method.展开更多
In this paper, we extend the mapping deformation method proposed by Lou. It is used to find new exacttravelling wave solutions of nonlinear partial differential equation or coupled nonlinear partial differential equat...In this paper, we extend the mapping deformation method proposed by Lou. It is used to find new exacttravelling wave solutions of nonlinear partial differential equation or coupled nonlinear partial differential equations(PDEs). Based on the idea of the homogeneous balance method, we construct the general mapping relation betweenthe solutions of the PDEs and those of the cubic nonlinear Klein-Gordon (NKG) equation. By using this relation andthe abundant solutions of the cubic NKG equation, many explicit and exact travelling wave solutions of three systemsof coupled PDEs, which contain solitary wave solutions, trigonometric function solutions, Jacobian elliptic functionsolutions, and rational solutions, are obtained.展开更多
Nonlinear partial differetial equation(NLPDE) is converted into ordinary differential equation(ODE) via a new ansatz.Using undetermined function method,the ODE obtained above is replaced by a set of algebraic equation...Nonlinear partial differetial equation(NLPDE) is converted into ordinary differential equation(ODE) via a new ansatz.Using undetermined function method,the ODE obtained above is replaced by a set of algebraic equations which are solved out with the aid of Mathematica.The exact solutions and solitary solutions of NLPDE are obtained.展开更多
In this paper, sane sufficient conditions are obtained for the oscillation for solutions of systems of high order partial differential equations of neutral type.
In this paper,by making use of the calculous technique and some results of the impulsive differential inequality,oscillatory properties of the solutions of certain nonlinear impulsive delay hyperbolic partial differen...In this paper,by making use of the calculous technique and some results of the impulsive differential inequality,oscillatory properties of the solutions of certain nonlinear impulsive delay hyperbolic partial differential equations with nonlinear diffusion coefficient are investigated.Sufficient conditions for oscillations of such equations are obtained.展开更多
The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.
This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut +αuux +βu^nuz +γuxx +δuzxx +ζxxxx = 0. This model...This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut +αuux +βu^nuz +γuxx +δuzxx +ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.展开更多
Novel distributed parameter neural networks are proposed for solving partial differential equations, and their dynamic performances are studied in Hilbert space. The locally connected neural networks are obtained by s...Novel distributed parameter neural networks are proposed for solving partial differential equations, and their dynamic performances are studied in Hilbert space. The locally connected neural networks are obtained by separating distributed parameter neural networks. Two simulations are also given. Both theoretical and computed results illustrate that the distributed parameter neural networks are effective and efficient for solving partial differential equation problems.展开更多
In this paper, with the aid of the symbolic computation, a further extended tanh function method was presented. Based on the new general ansatz, many nonlinear partial differential equation(s)(NPDE(s)) can he so...In this paper, with the aid of the symbolic computation, a further extended tanh function method was presented. Based on the new general ansatz, many nonlinear partial differential equation(s)(NPDE(s)) can he solved. Especially, as applications, a compound KdV-mKdV equation and the Broer-Kaup equations are considered successfully, and many solutions including periodic solutions, triangle solutions, and rational solutions are obtained. The method can also be applied to other NPDEs.展开更多
In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear...In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear equations by expanding the approximate solution in terms of shifted Legendre polynomials with unknown coefficients. The performance of presented method has been compared with other methods, namely Sinc-Galerkin, quadratic spline collocation and LiuLin method. Numerical examples show better accuracy of the proposed method. Moreover, the computation cost decreases at least by a factor of 6 in this method.展开更多
The first part of this article is an overview on some recent major developments in the field of analysis and partial different equations.It is a brief presentation given by the author at a round table discussion.The s...The first part of this article is an overview on some recent major developments in the field of analysis and partial different equations.It is a brief presentation given by the author at a round table discussion.The second part is a supplement of various details provided by several outstanding researchers on subjects.展开更多
In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomi...In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].展开更多
Partial differential equations arise in formulations of problems involving functions of several variables such as the propagation of sound or heat, electrostatics, electrodynamics, fluid flow, and elasticity, etc. The...Partial differential equations arise in formulations of problems involving functions of several variables such as the propagation of sound or heat, electrostatics, electrodynamics, fluid flow, and elasticity, etc. The present paper deals with a general introduction and classification of partial differential equations and the numerical methods available in the literature for the solution of partial differential equations.展开更多
文摘In this study,we aimto investigate certain triple integral transformand its application to a class of partial differentialequations.We discuss various properties of the new transformincluding inversion, linearity, existence, scaling andshifting, etc. Then,we derive several results enfolding partial derivatives and establish amulti-convolution theorem.Further, we apply the aforementioned transform to some classical functions and many types of partial differentialequations involving heat equations,wave equations, Laplace equations, and Poisson equations aswell.Moreover,wedraw some figures to illustrate 3-D contour plots for exact solutions of some selected examples involving differentvalues in their variables.
基金supported by the National Key R&D Program of China under Grant No.2021ZD0110400.
文摘Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.
基金the National Natural Science Foundation of China(Grant Nos.71961022,11902163,12265020,and 12262024)the Natural Science Foundation of Inner Mongolia Autonomous Region of China(Grant Nos.2019BS01011 and 2022MS01003)+5 种基金2022 Inner Mongolia Autonomous Region Grassland Talents Project-Young Innovative and Entrepreneurial Talents(Mingjing Du)2022 Talent Development Foundation of Inner Mongolia Autonomous Region of China(Ming-Jing Du)the Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region Program(Grant No.NJYT-20-B18)the Key Project of High-quality Economic Development Research Base of Yellow River Basin in 2022(Grant No.21HZD03)2022 Inner Mongolia Autonomous Region International Science and Technology Cooperation High-end Foreign Experts Introduction Project(Ge Kai)MOE(Ministry of Education in China)Humanities and Social Sciences Foundation(Grants No.20YJC860005).
文摘This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.
基金the support of the National Natural Science Foundation of China(Grant No.92152301)。
文摘Machine learning of partial differential equations(PDEs)from data is a potential breakthrough for addressing the lack of physical equations in complex dynamic systems.Recently,sparse regression has emerged as an attractive approach.However,noise presents the biggest challenge in sparse regression for identifying equations,as it relies on local derivative evaluations of noisy data.This study proposes a simple and general approach that significantly improves noise robustness by projecting the evaluated time derivative and partial differential term into a subspace with less noise.This method enables accurate reconstruction of PDEs involving high-order derivatives,even from data with considerable noise.Additionally,we discuss and compare the effects of the proposed method based on Fourier subspace and POD(proper orthogonal decomposition)subspace.Generally,the latter yields better results since it preserves the maximum amount of information.
基金supported by NSFC(No.11971296)National Key Research and Development Program of China(No.2021YFA1003004).
文摘We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural network frameworks.Concretely,we first perform Lagrange interpolation in front of the deep feedforward neural network.The Lagrange basis function has a neat structure and a strong expression ability,which is suitable to be a preprocessing tool for pre-fitting and feature extraction.Second,we introduce small sample learning into training,which is beneficial to guide themodel to be corrected quickly.Taking advantages of the theoretical support of traditional numerical method and the efficient allocation of modern machine learning,LaNets achieve higher predictive accuracy compared to the state-of-the-artwork.The stability and accuracy of the proposed algorithmare demonstrated through a series of classical numerical examples,including one-dimensional Burgers equation,onedimensional carburizing diffusion equations,two-dimensional Helmholtz equation and two-dimensional Burgers equation.Experimental results validate the robustness,effectiveness and flexibility of the proposed algorithm.
基金Supporting Project No.(RSP-2021/401),King Saud University,Riyadh,Saudi Arabia.
文摘The nonlinearity inmany problems occurs because of the complexity of the given physical phenomena.The present paper investigates the non-linear fractional partial differential equations’solutions using the Caputo operator with Laplace residual power seriesmethod.It is found that the present technique has a direct and simple implementation to solve the targeted problems.The comparison of the obtained solutions has been done with actual solutions to the problems.The fractional-order solutions are presented and considered to be the focal point of this research article.The results of the proposed technique are highly accurate and provide useful information about the actual dynamics of each problem.Because of the simple implementation,the present technique can be extended to solve other important fractional order problems.
基金Manar A.Alqudah would like to thank Princess Nourah bint Abdulrahman University Researchers Supporting Project No.(PNURSP2022R14),Princess Nourah bint Abdulrahman University,Riyadh,Saudi Arabia。
文摘The Laplace transformation is a very important integral transform,and it is extensively used in solving ordinary differential equations,partial differential equations,and several types of integro-differential equations.Our purpose in this study is to introduce the notion of fuzzy double Laplace transform,fuzzy conformable double Laplace transform(FCDLT).We discuss some basic properties of FCDLT.We obtain the solutions of fuzzy partial differential equations(both one-dimensional and two-dimensional cases)through the double Laplace approach.We demonstrate through numerical examples that our proposed method is very successful and convenient for resolving partial differential equations.
基金The National Natural Science Foundation of China(No.60972001)the National Key Technology R&D Program of China during the 11th Five-Year Period(No.2009BAG13A06)
文摘Due to the fact that the fourth-order partial differential equation (PDE) for noise removal can provide a good trade-off between noise removal and edge preservation and avoid blocky effects often caused by the second-order PDE, a domain-based fourth-order PDE method for noise removal is proposed. First, the proposed method segments the image domain into two domains, a speckle domain and a non-speckle domain, based on the statistical properties of isolated speckles in the Laplacian domain. Then, depending on the domain type, different conductance coefficients in the proposed fourth-order PDE are adopted. Moreover, the frequency approach is used to determine the optimum iteration stopping time. Compared with the existing fourth-order PDEs, the proposed fourth-order PDE can remove isolated speckles and keeps the edges from being blurred. Experimental results show the effectiveness of the proposed method.
文摘In this paper, we extend the mapping deformation method proposed by Lou. It is used to find new exacttravelling wave solutions of nonlinear partial differential equation or coupled nonlinear partial differential equations(PDEs). Based on the idea of the homogeneous balance method, we construct the general mapping relation betweenthe solutions of the PDEs and those of the cubic nonlinear Klein-Gordon (NKG) equation. By using this relation andthe abundant solutions of the cubic NKG equation, many explicit and exact travelling wave solutions of three systemsof coupled PDEs, which contain solitary wave solutions, trigonometric function solutions, Jacobian elliptic functionsolutions, and rational solutions, are obtained.
基金Supported by the Natural Science Foundation of Zhejiang Province(1 0 2 0 3 7)
文摘Nonlinear partial differetial equation(NLPDE) is converted into ordinary differential equation(ODE) via a new ansatz.Using undetermined function method,the ODE obtained above is replaced by a set of algebraic equations which are solved out with the aid of Mathematica.The exact solutions and solitary solutions of NLPDE are obtained.
文摘In this paper, sane sufficient conditions are obtained for the oscillation for solutions of systems of high order partial differential equations of neutral type.
基金Supported by the Natural Science Foundation of China(10471086)Supported by the Science Research Foundation of Department of Education of Hunan Province(07C164)
文摘In this paper,by making use of the calculous technique and some results of the impulsive differential inequality,oscillatory properties of the solutions of certain nonlinear impulsive delay hyperbolic partial differential equations with nonlinear diffusion coefficient are investigated.Sufficient conditions for oscillations of such equations are obtained.
文摘The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10661005)Fujian Province Science and Technology Plan Item (Grant No. 2008F5019)
文摘This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut +αuux +βu^nuz +γuxx +δuzxx +ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.
文摘Novel distributed parameter neural networks are proposed for solving partial differential equations, and their dynamic performances are studied in Hilbert space. The locally connected neural networks are obtained by separating distributed parameter neural networks. Two simulations are also given. Both theoretical and computed results illustrate that the distributed parameter neural networks are effective and efficient for solving partial differential equation problems.
基金supported by the Science Foundation of Shanghai Municipal Commission of Education (Grant No.06AZ081)the Science Foundation of Key Laboratory of Mathematics Mechanization (Grant No.KLMM0806)the shanghai Leading Academic Discipline Project (Grant No.J50101)
文摘In this paper, with the aid of the symbolic computation, a further extended tanh function method was presented. Based on the new general ansatz, many nonlinear partial differential equation(s)(NPDE(s)) can he solved. Especially, as applications, a compound KdV-mKdV equation and the Broer-Kaup equations are considered successfully, and many solutions including periodic solutions, triangle solutions, and rational solutions are obtained. The method can also be applied to other NPDEs.
文摘In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear equations by expanding the approximate solution in terms of shifted Legendre polynomials with unknown coefficients. The performance of presented method has been compared with other methods, namely Sinc-Galerkin, quadratic spline collocation and LiuLin method. Numerical examples show better accuracy of the proposed method. Moreover, the computation cost decreases at least by a factor of 6 in this method.
基金This article is initiated during a round table discussion at International Center for Mathematics in Southern University of Science and Technology(SUSTech)held in the January 2020.I am grateful to the wonderful organization and hospitality provided by Professors Tao Tang,Xiaoming Wang and Linlin Su.In the course of the preparation of this work,several colleagues have contributed to the scientific part of this review article.Here I would like to thank particularly:Chenjie Fan from University of Chicago,Jun Geng from Lanzhou University,Yanlin Liu from The Chinese Academy of Sciences,Xinan Ma from University of Sciences and Technology of China,Shuang Miao and Kelei Wang from Wuhan UniversityPing Zhang from The Chinese Academy of Sciences and Zhifei Zhang from Beijing University.Their enthusiasm and passion have greatly improved the presentation and the quality of this work.
文摘The first part of this article is an overview on some recent major developments in the field of analysis and partial different equations.It is a brief presentation given by the author at a round table discussion.The second part is a supplement of various details provided by several outstanding researchers on subjects.
文摘In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].
文摘Partial differential equations arise in formulations of problems involving functions of several variables such as the propagation of sound or heat, electrostatics, electrodynamics, fluid flow, and elasticity, etc. The present paper deals with a general introduction and classification of partial differential equations and the numerical methods available in the literature for the solution of partial differential equations.