The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundar...The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundary value problem of rectangular plates is proposed. The key concept behind this method is to transform the nonlinear or non-homogeneous part on the boundary into a lateral force within the governing function by the Dirac operator, which linearizes and homogenizes the original boundary, allowing one to employ the modal superposition method for obtaining solutions to reconstructive governing equations. Once projected into the modal space, the harmonic balance method(HBM) is utilized to solve coupled ordinary differential equations(ODEs)of truncated systems with nonlinearity. To validate the convergence and accuracy of the proposed Dirac method, the results of typical examples, involving nonlinearly restricted boundaries, moment excitation, and displacement excitation, are compared with those of the differential quadrature element method(DQEM). The results demonstrate that when dealing with nonlinear boundaries, the Dirac method exhibits more excellent accuracy and convergence compared with the DQEM. However, when facing displacement excitation, there exist some discrepancies between the proposed approach and simulations;nevertheless, the proposed method still accurately predicts resonant frequencies while being uniquely capable of handling nonuniform displacement excitations. Overall, this methodology offers a convenient way for addressing nonlinear and non-homogenous plate boundaries.展开更多
Clifford analysis is an important branch of modern analysis;it has a very important theoretical significance and application value,and its conclusions can be applied to the Maxwell equation,Yang-Mill field theory,quan...Clifford analysis is an important branch of modern analysis;it has a very important theoretical significance and application value,and its conclusions can be applied to the Maxwell equation,Yang-Mill field theory,quantum mechanics and value problems.In this paper,we first give the definition of a quasi-Cauchy type integral in complex Clifford analysis,and get the Plemelj formula for it.Second,we discuss the H?lder continuity for the Cauchy-type integral operators with values in a complex Clifford algebra.Finally,we prove the existence of solutions for a class of linear boundary value problems and give the integral representation for the solution.展开更多
In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be r...In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.展开更多
To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection techniq...To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).展开更多
We consider the singular Dirichlet problem for the Monge-Ampère type equation■=0,whereΩis a strictly convex and bounded smooth domain in■is positive and strictly decreasing in(0,∞)with■is positive inΩ.We ob...We consider the singular Dirichlet problem for the Monge-Ampère type equation■=0,whereΩis a strictly convex and bounded smooth domain in■is positive and strictly decreasing in(0,∞)with■is positive inΩ.We obtain the existence,nonexistence and global asymptotic behavior of the convex solution to such a problem for more general b and g.Our approach is based on the Karamata regular variation theory and the construction of suitable sub-and super-solutions.展开更多
This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of suff...This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of sufficient regularity of boundary conditions and coefficients, as long as the angle is sufficiently small, the regularity of the solution to the mixed boundary value problem of the second-order elliptic equation can reach any order.展开更多
Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with...Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.展开更多
We are concerned with global solutions of multidimensional(M-D)Riemann problems for nonlinear hyperbolic systems of conservation laws,focusing on their global configurations and structures.We present some recent devel...We are concerned with global solutions of multidimensional(M-D)Riemann problems for nonlinear hyperbolic systems of conservation laws,focusing on their global configurations and structures.We present some recent developments in the rigorous analysis of two-dimensional(2-D)Riemann problems involving transonic shock waves through several prototypes of hyperbolic systems of conservation laws and discuss some further M-D Riemann problems and related problems for nonlinear partial differential equations.In particular,we present four different 2-D Riemann problems through these prototypes of hyperbolic systems and show how these Riemann problems can be reformulated/solved as free boundary problems with transonic shock waves as free boundaries for the corresponding nonlinear conservation laws of mixed elliptic-hyperbolic type and related nonlinear partial differential equations.展开更多
k holomorphic functions are a type of generation of holomorphic functions.In this paper,a nonlinear boundary value problem for k holomorphic functions is primarily discussed on generalized polycylinders in C^(2).The e...k holomorphic functions are a type of generation of holomorphic functions.In this paper,a nonlinear boundary value problem for k holomorphic functions is primarily discussed on generalized polycylinders in C^(2).The existence of the solution for the problem is studied in detail with the help of the boundary properties of Cauchy type singular integral operators with a k holomorphic kernel.Furthermore,the integral representation for the solution is obtained.展开更多
In this paper,the weak pre-orthogonal adaptive Fourier decomposition(W-POAFD)method is applied to solve fractional boundary value problems(FBVPs)in the reproducing kernel Hilbert spaces(RKHSs)W_(0)^(4)[0,1] and W^(1)[...In this paper,the weak pre-orthogonal adaptive Fourier decomposition(W-POAFD)method is applied to solve fractional boundary value problems(FBVPs)in the reproducing kernel Hilbert spaces(RKHSs)W_(0)^(4)[0,1] and W^(1)[0,1].The process of the W-POAFD is as follows:(i)choose a dictionary and implement the pre-orthogonalization to all the dictionary elements;(ii)select points in[0,1]by the weak maximal selection principle to determine the corresponding orthonormalized dictionary elements iteratively;(iii)express the analytical solution as a linear combination of these determined dictionary elements.Convergence properties of numerical solutions are also discussed.The numerical experiments are carried out to illustrate the accuracy and efficiency of W-POAFD for solving FBVPs.展开更多
We study a special class of lower trigonometric matrix value boundary value problems on hyperbolas. Firstly, the pseudo-orthogonal polynomial on hyperbola is given in bilinear form and it is shown that it is the only ...We study a special class of lower trigonometric matrix value boundary value problems on hyperbolas. Firstly, the pseudo-orthogonal polynomial on hyperbola is given in bilinear form and it is shown that it is the only one. Secondly, a special boundary value problem of lower triangular matrix is presented and transformed into four related boundary value problems. Finally, Liouville theorem and Painlevé theorem and pseudo-orthogonal polynomials are used to give solutions.展开更多
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.展开更多
Aim The general arbitrary cracked problem in an elastic plane was discussed. Methods For the purpose of acquiring the solution of the problem, a new formulation on the problem was proposed. Compared with the classic...Aim The general arbitrary cracked problem in an elastic plane was discussed. Methods For the purpose of acquiring the solution of the problem, a new formulation on the problem was proposed. Compared with the classical plane elastic crack model, only the known conditions were revised in the new formulation, which are greatly convenient to solve the problem, and no other new condition was given. Results and Conclusion The general exact analytic solution is given here based on the formulation though the problem is very complicated. Furthermore, the stress intensity factors K Ⅰ, K Ⅱ of the problem are also given.展开更多
The singularly perturbed elliptic equation boundary value problem with turning point is considered. Using the method of multiple scales and the comparison theorem, the asymptotic behavior of solution for the boundary ...The singularly perturbed elliptic equation boundary value problem with turning point is considered. Using the method of multiple scales and the comparison theorem, the asymptotic behavior of solution for the boundary value problem is studied.展开更多
The existence of solutions for singular nonlinear two point boundary value problems subject to Sturm Liouville boundary conditions with p Laplacian operators is studied by the method of upper and lower solution...The existence of solutions for singular nonlinear two point boundary value problems subject to Sturm Liouville boundary conditions with p Laplacian operators is studied by the method of upper and lower solutions. The proof is based on an application of Schauder’s fixed point theorem to a modified problem whose solutions are that of the original one. At the same time, Arzela Ascoli theorem is used to prove that the defined operator N is a compact map.展开更多
By establishing equivalent fixed point theorem, the boundary value problems of p Laplace equations with finite time delay are studied. It’s the first time that the functional differential equation is discussed w...By establishing equivalent fixed point theorem, the boundary value problems of p Laplace equations with finite time delay are studied. It’s the first time that the functional differential equation is discussed with p Laplacian. The topological degree and fixed point theorem on cone are used to prove the existence of solution and positive solution. The conditions are all easy to check.展开更多
Aim To study a class of boundary value problem of hyperbolic partial functional differential equations with continuous deviating arguments. Methods An averaging technique was used. The multi dimensional problem was...Aim To study a class of boundary value problem of hyperbolic partial functional differential equations with continuous deviating arguments. Methods An averaging technique was used. The multi dimensional problem was reduced to a one dimensional oscillation problem for ordinary differential equations or inequalities. Results and Conclusion The known results of oscillation of solutions for a class of boundary value problem of hyperbolic partial functional differential equations with discrete deviating arguments are generalized, and the oscillatory criteria of solutions for such equation with two kinds of boundary value conditions are obtained.展开更多
The existence of multiple positive solutions for a class of higher order p Laplacian boundary value problem is studied. By means of the Leggett Williams fixed point theorem in cones, existence criteria which e...The existence of multiple positive solutions for a class of higher order p Laplacian boundary value problem is studied. By means of the Leggett Williams fixed point theorem in cones, existence criteria which ensure the existence of at least three positive solutions of the boundary value problem are established.展开更多
A class of nonlocal boundary value probl em s for elliptic systems in the unbounded domains are considered. Under suitable c onditions, the existence of solution and the comparison theorem for the boundary value prob...A class of nonlocal boundary value probl em s for elliptic systems in the unbounded domains are considered. Under suitable c onditions, the existence of solution and the comparison theorem for the boundary value problems are studied.展开更多
The present paper is concerned with the existence of positive solutions of the (k,n-k) conjugate boundary value problems(-1) n-k u (h) (t)=λa(t)f(u(t)),t∈(0,1), u (i) (0)=0,0≤i≤k-1, u (j) (0)=0,0...The present paper is concerned with the existence of positive solutions of the (k,n-k) conjugate boundary value problems(-1) n-k u (h) (t)=λa(t)f(u(t)),t∈(0,1), u (i) (0)=0,0≤i≤k-1, u (j) (0)=0,0≤j≤n-k-1,where λ is a positive parmeter. Krasnoselsii’s fixed point theorem is employed to obtain the existence criteria for positive solution.展开更多
基金Project supported by the National Natural Science Foundation of China (No. 12002195)the National Science Fund for Distinguished Young Scholars (No. 12025204)the Program of Shanghai Municipal Education Commission (No. 2019-01-07-00-09-E00018)。
文摘The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundary value problem of rectangular plates is proposed. The key concept behind this method is to transform the nonlinear or non-homogeneous part on the boundary into a lateral force within the governing function by the Dirac operator, which linearizes and homogenizes the original boundary, allowing one to employ the modal superposition method for obtaining solutions to reconstructive governing equations. Once projected into the modal space, the harmonic balance method(HBM) is utilized to solve coupled ordinary differential equations(ODEs)of truncated systems with nonlinearity. To validate the convergence and accuracy of the proposed Dirac method, the results of typical examples, involving nonlinearly restricted boundaries, moment excitation, and displacement excitation, are compared with those of the differential quadrature element method(DQEM). The results demonstrate that when dealing with nonlinear boundaries, the Dirac method exhibits more excellent accuracy and convergence compared with the DQEM. However, when facing displacement excitation, there exist some discrepancies between the proposed approach and simulations;nevertheless, the proposed method still accurately predicts resonant frequencies while being uniquely capable of handling nonuniform displacement excitations. Overall, this methodology offers a convenient way for addressing nonlinear and non-homogenous plate boundaries.
基金supported by the NSF of Hebei Province(A2022208007)the NSF of China(11571089,11871191)the NSF of Henan Province(222300420397)。
文摘Clifford analysis is an important branch of modern analysis;it has a very important theoretical significance and application value,and its conclusions can be applied to the Maxwell equation,Yang-Mill field theory,quantum mechanics and value problems.In this paper,we first give the definition of a quasi-Cauchy type integral in complex Clifford analysis,and get the Plemelj formula for it.Second,we discuss the H?lder continuity for the Cauchy-type integral operators with values in a complex Clifford algebra.Finally,we prove the existence of solutions for a class of linear boundary value problems and give the integral representation for the solution.
基金supported by the National Natural Science Foundation of China (No.12172154)the 111 Project (No.B14044)+1 种基金the Natural Science Foundation of Gansu Province (No.23JRRA1035)the Natural Science Foundation of Anhui University of Finance and Economics (No.ACKYC20043).
文摘In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.
基金supported by the the National Science and Technology Council(Grant Number:NSTC 112-2221-E239-022).
文摘To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).
基金supported by Shandong Provincial NSF(ZR2022MA020).
文摘We consider the singular Dirichlet problem for the Monge-Ampère type equation■=0,whereΩis a strictly convex and bounded smooth domain in■is positive and strictly decreasing in(0,∞)with■is positive inΩ.We obtain the existence,nonexistence and global asymptotic behavior of the convex solution to such a problem for more general b and g.Our approach is based on the Karamata regular variation theory and the construction of suitable sub-and super-solutions.
文摘This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of sufficient regularity of boundary conditions and coefficients, as long as the angle is sufficiently small, the regularity of the solution to the mixed boundary value problem of the second-order elliptic equation can reach any order.
文摘Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.
基金The research of Gui-Qiang G.Chen was supported in part by the UK Engineering and Physical Sciences Research Council Awards EP/L015811/1,EP/V008854/1,EP/V051121/1the Royal Society-Wolfson Research Merit Award WM090014.
文摘We are concerned with global solutions of multidimensional(M-D)Riemann problems for nonlinear hyperbolic systems of conservation laws,focusing on their global configurations and structures.We present some recent developments in the rigorous analysis of two-dimensional(2-D)Riemann problems involving transonic shock waves through several prototypes of hyperbolic systems of conservation laws and discuss some further M-D Riemann problems and related problems for nonlinear partial differential equations.In particular,we present four different 2-D Riemann problems through these prototypes of hyperbolic systems and show how these Riemann problems can be reformulated/solved as free boundary problems with transonic shock waves as free boundaries for the corresponding nonlinear conservation laws of mixed elliptic-hyperbolic type and related nonlinear partial differential equations.
基金the NSF of China(11571089,11871191)the NSF of Henan Province(222300420397)+1 种基金the NSF of Hebei Province(A2022208007)the Key Foundation of Hebei Normal University(L2018Z01)。
文摘k holomorphic functions are a type of generation of holomorphic functions.In this paper,a nonlinear boundary value problem for k holomorphic functions is primarily discussed on generalized polycylinders in C^(2).The existence of the solution for the problem is studied in detail with the help of the boundary properties of Cauchy type singular integral operators with a k holomorphic kernel.Furthermore,the integral representation for the solution is obtained.
基金University of Macao Multi-Year Research Grant Ref.No MYRG2016-00053-FST and MYRG2018-00168-FSTthe Science and Technology Development Fund,Macao SAR FDCT/0123/2018/A3.
文摘In this paper,the weak pre-orthogonal adaptive Fourier decomposition(W-POAFD)method is applied to solve fractional boundary value problems(FBVPs)in the reproducing kernel Hilbert spaces(RKHSs)W_(0)^(4)[0,1] and W^(1)[0,1].The process of the W-POAFD is as follows:(i)choose a dictionary and implement the pre-orthogonalization to all the dictionary elements;(ii)select points in[0,1]by the weak maximal selection principle to determine the corresponding orthonormalized dictionary elements iteratively;(iii)express the analytical solution as a linear combination of these determined dictionary elements.Convergence properties of numerical solutions are also discussed.The numerical experiments are carried out to illustrate the accuracy and efficiency of W-POAFD for solving FBVPs.
文摘We study a special class of lower trigonometric matrix value boundary value problems on hyperbolas. Firstly, the pseudo-orthogonal polynomial on hyperbola is given in bilinear form and it is shown that it is the only one. Secondly, a special boundary value problem of lower triangular matrix is presented and transformed into four related boundary value problems. Finally, Liouville theorem and Painlevé theorem and pseudo-orthogonal polynomials are used to give solutions.
文摘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.
文摘Aim The general arbitrary cracked problem in an elastic plane was discussed. Methods For the purpose of acquiring the solution of the problem, a new formulation on the problem was proposed. Compared with the classical plane elastic crack model, only the known conditions were revised in the new formulation, which are greatly convenient to solve the problem, and no other new condition was given. Results and Conclusion The general exact analytic solution is given here based on the formulation though the problem is very complicated. Furthermore, the stress intensity factors K Ⅰ, K Ⅱ of the problem are also given.
文摘The singularly perturbed elliptic equation boundary value problem with turning point is considered. Using the method of multiple scales and the comparison theorem, the asymptotic behavior of solution for the boundary value problem is studied.
文摘The existence of solutions for singular nonlinear two point boundary value problems subject to Sturm Liouville boundary conditions with p Laplacian operators is studied by the method of upper and lower solutions. The proof is based on an application of Schauder’s fixed point theorem to a modified problem whose solutions are that of the original one. At the same time, Arzela Ascoli theorem is used to prove that the defined operator N is a compact map.
文摘By establishing equivalent fixed point theorem, the boundary value problems of p Laplace equations with finite time delay are studied. It’s the first time that the functional differential equation is discussed with p Laplacian. The topological degree and fixed point theorem on cone are used to prove the existence of solution and positive solution. The conditions are all easy to check.
文摘Aim To study a class of boundary value problem of hyperbolic partial functional differential equations with continuous deviating arguments. Methods An averaging technique was used. The multi dimensional problem was reduced to a one dimensional oscillation problem for ordinary differential equations or inequalities. Results and Conclusion The known results of oscillation of solutions for a class of boundary value problem of hyperbolic partial functional differential equations with discrete deviating arguments are generalized, and the oscillatory criteria of solutions for such equation with two kinds of boundary value conditions are obtained.
文摘The existence of multiple positive solutions for a class of higher order p Laplacian boundary value problem is studied. By means of the Leggett Williams fixed point theorem in cones, existence criteria which ensure the existence of at least three positive solutions of the boundary value problem are established.
文摘A class of nonlocal boundary value probl em s for elliptic systems in the unbounded domains are considered. Under suitable c onditions, the existence of solution and the comparison theorem for the boundary value problems are studied.
文摘The present paper is concerned with the existence of positive solutions of the (k,n-k) conjugate boundary value problems(-1) n-k u (h) (t)=λa(t)f(u(t)),t∈(0,1), u (i) (0)=0,0≤i≤k-1, u (j) (0)=0,0≤j≤n-k-1,where λ is a positive parmeter. Krasnoselsii’s fixed point theorem is employed to obtain the existence criteria for positive solution.