In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element...This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element method based on piecewise linear approximating functions and we also use the barycentric quadrature rule to compute the stiffness matrix and the L2-norm.展开更多
To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltai...To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.展开更多
To solve the receding horizon control （RHC） problem in an online manner, a novel numerical method called the indirect Radau pseudospectral method （IRPM） is proposed in this paper. Based on calculus of variations a...To solve the receding horizon control （RHC） problem in an online manner, a novel numerical method called the indirect Radau pseudospectral method （IRPM） is proposed in this paper. Based on calculus of variations and the first-order necessary optimality condition, the RHC problem for linear time-varying （LTV） system is transformed into the two-point boundary value problem （TPBVP）. The Radau pseudospectral approximation is employed to discretize the TPBVP into well-posed linear algebraic equations. The resulting linear algebraic equations are solved via a matrix partitioning approach afterwards to obtain the optimal feedback control law. For the nonlinear system, the linearization method or the quasi linearization method is employed to approximate the RHC problem with successive linear approximations. Subsequently, each linear problem is solved via the similar method which is used to solve the RHC problem for LTV system. Simulation results of three examples show that the IRPM is of high accuracy and of high compu- tation efficiency to solve the RHC problem and the stability of closed-loop systems is guaranteed.展开更多
The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″（t） ＋ q（t）f（t,u（t）,u′（t）） = 0, 0 〈 t 〈 1, u（0） = u″（0） = u′（1） = 0 is studied....The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″（t） ＋ q（t）f（t,u（t）,u′（t）） = 0, 0 〈 t 〈 1, u（0） = u″（0） = u′（1） = 0 is studied. The iterative schemes for approximating the solutions are obtained by applying a monotone iterative method.展开更多
This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matr...This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.展开更多
To overcome the inaccuracy problem of the traditional wedge evaluation of steel plates and strips caused by the ran-domness of the thicknesses of two local points and improve the reliability of the wedge index,the dou...To overcome the inaccuracy problem of the traditional wedge evaluation of steel plates and strips caused by the ran-domness of the thicknesses of two local points and improve the reliability of the wedge index,the double-centroid method for the wedge evaluation was proposed,and a model based on the centroid theory was established.Meanwhile,an integral model for the discrete thickness values of the cross-section profiles was derived.The discussion focused on the distinct characteristics of the two-point method,asymmetric method,and double-centroid method in evaluating the asymmetric distribution of cross-sections.The three methods were employed to evaluate the wedge values of both the theoretical and measured cross-sections of steel plates and strips,and the accuracies of three wedge evaluation models were analyzed and discussed.The results showed that the double-centroid method objectively reflects the degree and variation characteristics of the wedge values of the cross-sections of steel plates and strips,and this method is feasible,reliable,and outstanding.展开更多
In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, w...In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.展开更多
This paper presents an enhanced version of the standard shooting method that enables problems with two unknown parameters to be solved.A novel approach is applied to the analysis of the natural vibrations of Euler-Ber...This paper presents an enhanced version of the standard shooting method that enables problems with two unknown parameters to be solved.A novel approach is applied to the analysis of the natural vibrations of Euler-Bernoulli beams.The proposed algorithm,named as two-parameter multiple shooting method,is a new powerful numerical tool for calculating the natural frequencies and modes of multi-segment prismatic and non-prismatic beams with different boundary conditions.The impact of the axial force and additional point masses is also taken into account.Due to the fact that the method is based directly on the fourth-order ordinary differential equation,the structures do not have to be divided into many small elements to obtain an accurate enough solution,even though the geometry is very complex.To verify the proposed method,three different examples are considered,i.e.,a three-segment non-prismatic beam,a prismatic column subject to non-uniformly distributed compressive loads,and a two-segment beam with an additional point mass.Numerical analyses are carried out with the software MATHEMATICA.The results are compared with the solutions computed by the commercial finite element program SOFiSTiK.Good agreement is achieved,which confirms the correctness and high effectiveness of the formulated algorithm.展开更多
In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis sh...In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.展开更多
For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has ...For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has the accuracy O（h^min{2k,k＋4}） The theoretical analysis coincides the reported numerical results.展开更多
We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the p...We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the problem into an equivalent system of weakly singular integral equations.Then,a Legendre-based spectral collocation method is developed for solving the transformed system.Therefore,we can make good use of the advantages of the Gauss quadrature rule.We present the construction and analysis of the collocation method.These results can be indirectly applied to solve fractional optimal control problems by considering the corresponding Euler–Lagrange equations.Two numerical examples are given to confirm the convergence analysis and robustness of the scheme.展开更多
Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly pertur...Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems （TPBVPs） with one boundary layer. First, the inhomogeneous ordinary differential equations （ODEs） are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.展开更多
This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singula...This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.展开更多
Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Br...Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Bratu’s equation, Troesch’s problems) occurs engineering and science, including the modeling of chemical reactions diffusion processes and heat transfer. An analytical expression pertaining to the concentration of substrate is obtained using Homotopy perturbation method for all values of parameters. These approximate analytical results were found to be in good agreement with the simulation results.展开更多
The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and...The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.展开更多
The objective of the present study is to explore the relation between the near-wall vortices and the shear stress on the wall in two-dimensional channel flows. A direct numerical simulation of an incompressible two-di...The objective of the present study is to explore the relation between the near-wall vortices and the shear stress on the wall in two-dimensional channel flows. A direct numerical simulation of an incompressible two-dimensional turbulent channel flow is performed with spectral method and the results are used to examine the relation between wall shear stress and near-wall vortices. The two-point correlation results indicate that the wall shear stress is associated with the vortices near the wall and the maximum correlation-value location of the near-wall vortices is obtained. The analysis of the instantaneous diagrams of fluctuation velocity vectors provides a further expression for the above conclusions. The results of this research provide a useful supplement for the control of turbulent boundary layers.展开更多
Matrix methods, now-a-days, are playing an important role in solving the real life problems governed by ODEs and/or by PDEs. Many differential models of sciences and engineers for which the existing methodologies do n...Matrix methods, now-a-days, are playing an important role in solving the real life problems governed by ODEs and/or by PDEs. Many differential models of sciences and engineers for which the existing methodologies do not give reliable results, these methods are solving them competitively. In this work, a matrix methods is presented for approximate solution of the second-order singularly-perturbed delay differential equations. The main characteristic of this technique is that it reduces these problems to those of solving a system of algebraic equations, thus greatly simplifying the problem. The error analysis and convergence for the proposed method is introduced. Finally some experiments and their numerical solutions are given.展开更多
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
文摘This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element method based on piecewise linear approximating functions and we also use the barycentric quadrature rule to compute the stiffness matrix and the L2-norm.
基金This research is supported by the Deputyship forResearch&Innovation,Ministry of Education in Saudi Arabia under Project Number(IFP-2022-35).
文摘To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.
基金supported by the National Natural Science Foundation of China(Nos.61174221 and 61402039)
文摘To solve the receding horizon control （RHC） problem in an online manner, a novel numerical method called the indirect Radau pseudospectral method （IRPM） is proposed in this paper. Based on calculus of variations and the first-order necessary optimality condition, the RHC problem for linear time-varying （LTV） system is transformed into the two-point boundary value problem （TPBVP）. The Radau pseudospectral approximation is employed to discretize the TPBVP into well-posed linear algebraic equations. The resulting linear algebraic equations are solved via a matrix partitioning approach afterwards to obtain the optimal feedback control law. For the nonlinear system, the linearization method or the quasi linearization method is employed to approximate the RHC problem with successive linear approximations. Subsequently, each linear problem is solved via the similar method which is used to solve the RHC problem for LTV system. Simulation results of three examples show that the IRPM is of high accuracy and of high compu- tation efficiency to solve the RHC problem and the stability of closed-loop systems is guaranteed.
基金Supported by the Natural Science Foundation of Zhejiang Province (Y605144)the XNF of Zhejiang University of Media and Communications (XN080012008034)
文摘The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″（t） ＋ q（t）f（t,u（t）,u′（t）） = 0, 0 〈 t 〈 1, u（0） = u″（0） = u′（1） = 0 is studied. The iterative schemes for approximating the solutions are obtained by applying a monotone iterative method.
基金Project supported by the National Natural Science Foundation of China(No.10672194)the China-Russia Cooperative Project(the National Natural Science Foundation of China and the Russian Foundation for Basic Research)(No.10811120012)
文摘This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.
基金supported by the National Natural Science Foundation of China (Grant No.52305405)the Natural Science Foundation Research Program of Shanxi Province (Grant No.202203021222121)+4 种基金the Major Project of Science and Technology of Shanxi Province (Grant No.20181102016)the Chinese Postdoctoral Science Foundation (Grant No.2021M702544)the Central Government Guides the Special Fund Projects of Local Scientific and Technological Development (YDZX20191400002149)the Open Project of Research Institute of Hai'an-Taiyuan University of Technology (Grant No.2023HA-TYUTKFYF008)the School Fund of Taiyuan University of Technology (Grant No.2022QN007).
文摘To overcome the inaccuracy problem of the traditional wedge evaluation of steel plates and strips caused by the ran-domness of the thicknesses of two local points and improve the reliability of the wedge index,the double-centroid method for the wedge evaluation was proposed,and a model based on the centroid theory was established.Meanwhile,an integral model for the discrete thickness values of the cross-section profiles was derived.The discussion focused on the distinct characteristics of the two-point method,asymmetric method,and double-centroid method in evaluating the asymmetric distribution of cross-sections.The three methods were employed to evaluate the wedge values of both the theoretical and measured cross-sections of steel plates and strips,and the accuracies of three wedge evaluation models were analyzed and discussed.The results showed that the double-centroid method objectively reflects the degree and variation characteristics of the wedge values of the cross-sections of steel plates and strips,and this method is feasible,reliable,and outstanding.
基金The Project was supported by National Natural Science Foundation of China
文摘In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.
文摘This paper presents an enhanced version of the standard shooting method that enables problems with two unknown parameters to be solved.A novel approach is applied to the analysis of the natural vibrations of Euler-Bernoulli beams.The proposed algorithm,named as two-parameter multiple shooting method,is a new powerful numerical tool for calculating the natural frequencies and modes of multi-segment prismatic and non-prismatic beams with different boundary conditions.The impact of the axial force and additional point masses is also taken into account.Due to the fact that the method is based directly on the fourth-order ordinary differential equation,the structures do not have to be divided into many small elements to obtain an accurate enough solution,even though the geometry is very complex.To verify the proposed method,three different examples are considered,i.e.,a three-segment non-prismatic beam,a prismatic column subject to non-uniformly distributed compressive loads,and a two-segment beam with an additional point mass.Numerical analyses are carried out with the software MATHEMATICA.The results are compared with the solutions computed by the commercial finite element program SOFiSTiK.Good agreement is achieved,which confirms the correctness and high effectiveness of the formulated algorithm.
基金supported partially by the innovation fund of Shanghai Normal Universitysupported partially by NSERC of Canada under Grant OGP0046726.
文摘In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.
基金Project supported by the National Natural Science Foundation of China (Nos. 10571046, 10371038)
文摘For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has the accuracy O（h^min{2k,k＋4}） The theoretical analysis coincides the reported numerical results.
基金The Russian Foundation for Basic Research(RFBR)Grant No.19-01-00019.
文摘We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the problem into an equivalent system of weakly singular integral equations.Then,a Legendre-based spectral collocation method is developed for solving the transformed system.Therefore,we can make good use of the advantages of the Gauss quadrature rule.We present the construction and analysis of the collocation method.These results can be indirectly applied to solve fractional optimal control problems by considering the corresponding Euler–Lagrange equations.Two numerical examples are given to confirm the convergence analysis and robustness of the scheme.
基金Project supported by the National Natural Science Foundation of China(Key Program)(Nos.11132004 and 51078145)
文摘Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems （TPBVPs） with one boundary layer. First, the inhomogeneous ordinary differential equations （ODEs） are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.
基金supported by the National Natural Science Foundation of China(Key Program)(Nos.11132004 and 51078145)
文摘This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.
文摘Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Bratu’s equation, Troesch’s problems) occurs engineering and science, including the modeling of chemical reactions diffusion processes and heat transfer. An analytical expression pertaining to the concentration of substrate is obtained using Homotopy perturbation method for all values of parameters. These approximate analytical results were found to be in good agreement with the simulation results.
文摘The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.
文摘The objective of the present study is to explore the relation between the near-wall vortices and the shear stress on the wall in two-dimensional channel flows. A direct numerical simulation of an incompressible two-dimensional turbulent channel flow is performed with spectral method and the results are used to examine the relation between wall shear stress and near-wall vortices. The two-point correlation results indicate that the wall shear stress is associated with the vortices near the wall and the maximum correlation-value location of the near-wall vortices is obtained. The analysis of the instantaneous diagrams of fluctuation velocity vectors provides a further expression for the above conclusions. The results of this research provide a useful supplement for the control of turbulent boundary layers.
文摘Matrix methods, now-a-days, are playing an important role in solving the real life problems governed by ODEs and/or by PDEs. Many differential models of sciences and engineers for which the existing methodologies do not give reliable results, these methods are solving them competitively. In this work, a matrix methods is presented for approximate solution of the second-order singularly-perturbed delay differential equations. The main characteristic of this technique is that it reduces these problems to those of solving a system of algebraic equations, thus greatly simplifying the problem. The error analysis and convergence for the proposed method is introduced. Finally some experiments and their numerical solutions are given.
