Mathematical model derived from multi-dimensional differential equations play a crucial role in modeling a variety of scientific and engineering application problems. The study of vibration of membrane of different geometries has been the subject of interest to scientists and engineers since long back owing to their wide applications in every field of modern technology. Different shapes or geometries of membranes are the fundamental constituents in many scientific states. Many researchers concentrated on membrane structures and studied about the vibration theory of membrane employing numerous analytical or numerical techniques. However, to the best of our knowledge, none had attempted to provide approximate eigenvalues of boundary value problems employing Bernstein-Galerkin method of weighted residual (MWR).
Some physical problems arise in science and engineering are modeled by the elliptic, parabolic and hyperbolic partial differential equations. The process of solving certain linear evolution equations such as the heat or wave equations have received increased interest from several fields of science and engineering, either on bounded domains. Certain class of linear evolution equations for instance the heat or wave equations have a second order linear differential eigenvalue problem with two boundary conditions but no unique solution exists. A code for second order problems has been implemented and special attention has been paid to the approximation of the boundary conditions in the two-dimensional case.
A new numerical method based on Bernstein polynomials expansion is proposed for solving one dimensional (1-D) elliptic interface problems . At present, formulations exploiting Galerkin method of weighted residual are constructed to determine the coefficient matrix derived from linear system of equations. Helmholtz equation appears in vibrating-string problems and in finding the temperature distribution in a rod. Calculations of eigenvalues of Helmholtz equation using boundary method are presented in . Numerical technique utilizing Bernstein polynomials basis to give the approximate solution of a parabolic partial differential equation illustrated in . The large deflection theory of membrane is applied  to derive the governing vibration equations of orthotropic membrane. The Bernstein polynomials are used to find numerical solution of elliptic boundary value problems with a point collocation method . The authors solved the problem and obtained the power series formula of nonlinear vibration frequency of rectangular membrane with four edges fixed. Bernstein Ritz-Galerkin method for solving an initial boundary value problem for two-dimensional (2-D) wave equation has been studied by . A spectral element model for the transverse vibration of a membrane has been developed by  using the boundary splitting method and the waveguide FEM-based spectral super element method (SSEM). Exact solutions are well known for the free vibration of a membrane with simple geometry such as elliptical , circular , general convex polygonal membranes , rectangular membranes  - . Nevertheless, closed-form solutions rarely exist which leads many researchers to have numerical solution applying several techniques in view to minimize the error with less computational difficulties. Free vibrations of membranes have been studied in the literature by applying the method of superposition using wave-type base functions , the optimized Galerkin-Kantorovich method , classical Rayleigh-Ritz method , finite difference method , the Hybrid method , the discrete singular advance convolution method , the Kantorovich method . Furthermore, the external and internal excitation methods , are also investigated. Two-dimensional (2-D) Helmholtz equation is solved for non-homogeneous membrane exploiting Sinc collocation method .
Double ultraspherical spectral Galerkin method with Chebychev polynomials is applied to solve elliptic partial differential equations (PDEs) . Chebychev spectral quasi-inverse matrix diagonalization Galerkin method  is utilized to solve multidimensional PDE’s. Finite difference methods have been offered for the numerical solution of the one-dimensional parabolic equation in the study . Eigenvalue analysis for second order boundary value problems arises in several engineering application areas has been investigated using analytical or numerical techniques . Bernstein basis has been exploited to solve ordinary and partial differential to solve eigenvalues as well as boundary value problems employing numerous techniques are in  - .
In this study, we have presented Bernstein polynomials based Galerkin method of weighted residual (MWR) technique that offers accurate solutions, are put up with in terms of truncated series of smooth polynomial functions. The key advantages of applying the Galerkin MWR are its superior accuracy which can be achieved even with a few numbers of polynomials compared to the other techniques reported in the proposed work.
The significant benefit of the present method is that the trial Bernstein basis function can be modified effortlessly to any desired form so as to satisfy the Dirichlet and derivative boundary conditions for both one and two-dimensional problems.
Our goal is to construct an appropriate method with Dirichlet and derivative type boundary conditions and to reduce the errors as well as the computational difficulties. The main reason why the Galerkin MWR is chosen is its flexibility and simple implementation. The notable properties of Bernstein polynomials vanishing at its endpoints over the finite interval inspire us to execute it in the Galerkin technique. Furthermore, the smallest eigenvalue which characterizes potentially the most visual structures of the dynamical systems can be computed very accurately applying the said technique which converges to large significant digits and fairly close to the analytical results.
Some analytical methods involve a reduced amount of memory and arithmetic, but usually they are unable to fix the maximal dimension for the spectrum. Perturbation methods limit the extent of its application due to disadvantages in its theory in its existing form and all types of techniques required small parameter and consequently require some skill to implement. Finite difference methods (FDM) and finite elements methods (FEM) are well-known discretization techniques that comprise extensive calculation and huge memory storage in machine for storing matrices having larger dimensions. Additionally, the precision rapidly worsens for the upper eigenvalues. Furthermore, FDM gives value at particular points ineffective to evaluate the values at the desired points between two grid points and takes more computational cost for getting higher accuracy. The co-efficient matrix in Galerkin MWR is sparse and has symmetric banded matrix, which minimizes the computational effort and attains relatively smaller condition numbers. Although discretization methods attain high accuracy, in the end produce ill-conditioned system. On the other hand, the obtained matrices are full and non-symmetric which intensify the computational cost.
In section 2, we provide an outline of Bernstein polynomial basis in single and multi-dimensions. Approximation of functions to Bernstein basis functions and their derivatives required for numerical computations are also illustrated. Section 3 is offered to the construction of 2-D Helmholtz eigenvalue problems and formulation of one dimensional transient heat equation. In section 4, we give short description of hyperbolic equation, analytical solution of 2-D wave equation and Galerkin formulation of inhomogeneous rectangular membrane. Convergence of Galerkin method and error estimation is discussed in section 5. In section 6, we represent some numerical examples eigenvalue problems in one and two dimensions. In section 7, we illustrate some results with detailed discussions. Moreover, we validate our code by the offered method and eigenvalues are compared to the analytical ones. We have drawn conclusions of this study in section 8.
2. Bernstein Polynomials
2.1. Bernstein Polynomials in One Dimension  
, , (1)
For each positive integer n, there is a sequence of Bernstein polynomials over the finite interval . Well convergence of Bernstein polynomials for some class of functions is found in . Bernstein polynomials which approximates a functions of one real variable is defined as
2.2. Two Dimensional Bernstein Polynomials
2-D Bernstein polynomials approximates a functions of two real variables is defined as
where , , ; with .
There are and -th degree Bernstein polynomials. For suitability, we usually set , if , or , .
Some interesting and useful properties of 2-D Bernstein polynomials are demonstrated in .
2.3. The Bernstein Approximation Theorem 1 
Every continuous function f defined on [0, 1] can be uniformly approximated as closely as desired by a polynomial function. For any , there exists a positive integer N such that for all , an integer we have
Hence given any power-form polynomial of degree N, it can be uniquely converted into a Bernstein polynomial of degree n for .
Let is a continuous function space and is a linear norm space on with the norm and all partial derivatives of exists, while is a sequence on the space , .
which implies the uniform convergence to the function . Let
and , .
Let ; . Since V is a finite dimensional vector space, f has a unique best approximation from V say such that .
Lemma: Suppose that is a continuous differentiable function and all partial derivatives of exist and are continuous. If is the best approximation for f out of V then the error bound is represented as
as , and the partial derivatives of f are bounded in magnitude by . By applying the Taylor’s expansion in two variables we define
2.4. Derivatives of Multivariable Bernstein Polynomials
We can define p-th order partial derivative in x-axis direction in a following way
Also q-th order partial derivative in y-axis direction in a following way
Finally derivative of order is defined as follows
3. Description of 2-D Helmholtz Problem
We consider the governing transient heat transfer on a two dimensional region  .
with the boundary conditions (b. c’s)
on . (11a)
and the initial condition
in . (11b)
where t specifies for the time and c, , , , , f, , are regarded as functions of positions and/or time.
The corresponding homogeneous form of the above Equation (11) is
on . (12)
Approximate solution of Equation (12) is given as
Thus Equation (12) becomes
and boundary conditions take the form , . (13a)
Here the eigenvalues ’s and eigenfunctions U’s satisfy the time independent form of Equation (13). Then Equation (13) is known as Helmholtz equation which can be solved using separation of variables technique. Now the approximate solution of Equation (13) is
where denotes spatial variable and parameter defines the function of time only.
Galerkin Formulation for One Dimensional Heat Equation
Here we applied a modified weighted residual technique which involves only first derivative instead of second derivative terms in the residual equation. This can be performed applying integration by parts in which the operational matrices for integration and the product is utilized to convert the given differential eigenvalue problems to a system of algebraic equations.
One dimensional parabolic partial differential equation
with the boundary and initial conditions
which governs the transient heat transfer in one-dimensional system (i.e., a plane wall), here u denotes the temperature, k the thermal conductivity, the density, A the cross-sectional area, c the specific heat, q the heat generation per unit length. Here represents the temperature of a rod for both space and time and is trial or basis function. If non-homogeneous heat Equation (15) with a lateral heat loss reduces to homogeneous equation (q = 0).
Let us assume the approximate solution of Equation (16)
The solution of Equation (15) can be written as
The Galerkin weighted residual equation can be written as
On simplification the by using the conditions the above residual Equation (17) takes the matrix form
where refers to the time derivative.
Hence the matrix form of Equation (18)
where refers to the time derivative.
The solution of Equation (18) can be given as
Hence the matrix form of (18)
Solving (19) we can determine and nonzero such that (19) holds and the boundary conditions of the problem met. Equation (19) describe an eigenvalue problem; being the eigenvalues and be the eigenfunction.
4. Brief Description of Hyperbolic Equations
The transverse motion of a membrane, is governed by a partial differential equation of the form as given in (14), where, denotes the transverse deflection, c is the material density of the membrane, and are the tensions in the x and y directions of the membrane, is the modulus of elastic foundation on which the membrane is stretched (often, ) i.e., there is no foundation), and is the transversely distributed force. The above equation is known as wave equation, and is classified mathematically as a hyperbolic equation. The function must be determined such that it satisfies Equation (19) in a region together with specified boundary and initial conditions. The problem of finding the solution as given by Equation (14) holds for homogeneous boundary and initial conditions and is called an eigenvalue problem.
4.1. Solution of Two-Dimensional Wave Equation in Rectangular Membrane
We consider the equation for the vibration of a tightly stretched membrane as the membrane of a drum. We assume that the membrane is of uniform tension and the tension per unit length is same in all directions at every point and m be the mass of per unit of area of the membrane.
The fact that we are keeping the edges of the membrane fixed is expressed by the boundary conditions
We must also specify how the membrane is initially deformed and set into motion. This is done via the initial conditions
The solution of Equation (21) applying the separation of variables gives the particular values of l and n.
These are the solutions of the wave Equation (22) which are zero on the boundary of the rectangular membrane. These functions are known as eigenfunctions and the numbers are known eigenvalues.
The general solution is
4.2. Vibration of Distributed Inhomogeneous System in a Rectangular Domain
We consider free vibrations of a uniformly stretched inhomogeneous membrane of rectangular shape with clamped edges in an elastic medium. We consider the following eigenvalue problem
where the boundary conditions are
and . (24a)
Let denotes an approximation to the exact solution.
and , is specified by the Dirichlet boundary conditions, of Equation (24a).
where is the residual function.
The Green’s theorem states that if A is a region in the xy plane bounded by a closed curve then we suppose:
is a continuous vector field defined on a region A in Ω. Moreover, suppose F and G have continuous partial derivatives and that the boundary is a simple closed curve C with positive orientation. Then
We can write
Using Equations (26a), (27) in Equation (26b)
We finally obtain
5. Error and Convergence of Galerkin MWR
Let be the set of normalized eigenfunctions of the Sturm-Liouville system. If f be a function in Hilbert space H, then
We consider the linear operator L with a field of definition , that is, Lu is defined for . The inner product is
The operator is symmetric if for elements
The operator is positive definite if for any function in , not identically zero:
and is positive bounded below if for any ;
, for .
Uniform Convergence: For any , we can find n such that
Convergence in energy requires where energy is defined as
The corresponding error bounds are point wise error in Equation (31). A sequence , converges weakly to an element u of a space if holds for all in the space. Sequences produced by the Galerkin method sometimes converge weakly to a generalized solution.
A set of trial functions form a complete set of functions and they are complete in a space if any function in the space can be expanded in terms of the set of functions, for sufficiently large n .
Since the powers of Bernstein polynomials ( ) have sufficient continuity and they form an infinite sequence of functions which is complete. Also the convergence of the eigenvalues by Galerkin WRM method is measured by the relative error.
depends upon the problems.
We consider a two dimensional boundary value problem.
is non-self-adjoint unless . For convergence of the Galerkin method, we assume there exists a constant such that for any point in V and real numbers ,
The coefficients , and their first derivatives are continuous in and the coefficients C and F are continuous in . We assume the problem has a unique solution. The proof of the theorem is illustrated in and the references therein . Here is a finite bounded domain and .
where the operator is positive bounded below for the set of functions which vanishes on S. The convergence proofs are done in terms of the operator and the remaining terms in the differential equation are bounded in terms of this operator and its inverse.
Theorem2: Assume the problem (11) is unique for the boundary conditions on , on , on .
Exploiting the Galerkin method using trial functions which are complete in energy of , are linearly independent, and satisfy on . Then the approximate solution , converges in energy of to the solution.
Let us consider the two-dimensional the following eigenvalue problem
on S where , satisfies Equation (34).
Theorem 3: If the Galerkin method is applied to Equation (36) using a system of functions which is complete in energy and linearly independent, then the approximate eigenvalues converge to the true eigenvalue. This proves convergence of the Galerkin method in section 5 when applied to Equation (36).
6. Numerical Experiments
In this section we apply Galerkin MWR exploiting Bernstein polynomial only to the eigenvalue problems frequently encountered in science and engineering in one and two dimensions. For validation of our proposed technique formulated in section 4, we demonstrate three linear problems in one dimension and two linear problems in two dimensions.
In this section we apply the current techhnique exploiting Bernstein polynomial only to the eigenvalue problems frequently encountered in science and engineering in one and two dimensions. For validation of our proposed technique formulated in section 4, we demonstrate three linear problems in one dimension and two linear problems in two dimensions.
Example 1: To determine the boundaries of parametric instability regions for the following modified Mathieu equation using the method of accelerated convergence which is described in details in .
, , .
For , Equation (37) differs from the classical Mathieu equation.
To compute the eigenvalues depending on parameter e which determines the boundaries of the parametric instability , we need to construct periodic solutions of that:
, . (37a)
Since the function is even with respect to t, in (the conditions of perioiodicity in (37a) are equivalent to the conditions of the first kind given as
and the second kind given as
here the vibration of elastic Crankshaft with concentrated load yield the following equation with periodic co-efficient and modulation depth e illustrated in  is considered.
, . (38)
Example 2: Let us consider one dimensional Helmholtz equation .
We consider a plane wall, initially at a uniform temperature .
Set 1: , , (40a)
Set 2: , . (40b)
Table 2 displays the first ten eigenvalues for n = 15, exploiting Bernstein polynomials.
Example 4: The well-known transient heat conduction   in a square region with boundary and initial conditions is considered.
, . (41)
Subject to the boundary conditions
The boundary conditions are:
, . (41a)
and . (41b)
and the initial condition is
Homogeneous form of (41) for is
The derivative boundary conditions will be applied integrand of the residual equation. Hence for executing only the Dirichlet boundary conditions in Equation (41a), we modify the polynomial basis as
The exact solution of (42) is , where refers to the required eigenvalue which corresponds to the eigenfunctions . The analytical eigenvalue ,
Example 5: The homogeneous membrane of sides 4 ft by 2 ft, fixed on all its four edges i.e., on is considered . We assume that the tension of the membrane is (i.e. ). The exact natural frequencies of the rectangular membrane are obtained by considering the full domain
The fact that we are keeping the edges of the membrane fixed is expressed by the boundary conditions
and the initial conditions
For implementing boundary conditions in Equations (43a), (43b) we modify the polynomial basis as
7. Result and Discussions
In Table 1, we computed even and odd eigen solutions for different modulation depth e ranges from 0.1 to 0.9. Analytical eigen-solutions using perturbation method are listed in . Our computed relative errors for the smallest modes are less than 5% which shows that the present technique produced significantly small errors. From Table 2, the smallest eigenvalues attains the accuracy up to 10−16 and error increases rapidly for higher eigenvalues than the lower values which is better than boundary method. As we increase the grid points or nodes from n = 15 to n = 30, the error decays fast for all the eigenvalues and consistent accuracy is obtained up to 10−13.
We observed that increasing of nodes reveal the stable behaviour of all the eigenvalues for n = 30. It is also observed that our present approach accomplishes accurate results and is compatible to the existence new boundary approach for one dimensional Helmholtz equation. Comparison of relative errors obtained for Gal. Relative errors between our proposed Galerkin MWR and FEM  are compared which are displayed in Figure 1. In Table 3(a) and Table 3(b) we computed seven eigen frequencies for the two sets of boundary conditions exploiting current formulation and comparison to the exact and finite element results are shown with relative errors. From these tables, we observed that the smallest error decreases up to which is much smaller than that of FEM . Eigen frequencies using different degree of polynomial are displayed in Table 4 of example 4. We noticed that the smallest order of the accuracies for the largest eigenvalues are 10−2, 10−3, 10−4 respectively with 4, 8, 10, 12 Bernstein polynomials and the first three eigenvalues converge to the exact results for n = 13 illustrated in Table 4. This confirms that the accuracy increases with the increased degree of polynomials. Our current method produces small errors in percentages (%) than those attained by FEM .
Plot of the relative errors between exact and approximate eigenvalues are depicted in Figure 2 for different values of n. Also, the relative errors obtained by Galerkin MWR for different degree of Bernstein polynomials are depicted in Figure 3. Table 5 shows that the maximum relative error for the smallest and the largest eigenvalues using present method are 10−6 and 10−7 respectively, whereas applying FEM , the errors obtained as 10−3, respectively. From the results in tabular form we conclude that the Galerkin MWR is much accurate, efficient and compatible with other techniques.
Table 1. Comparison of eigenvalues for vibration of crankshaft of Example 1.
Table 2. Comparison of relative errors obtained by Galerkin MWR with for Example 2.
Table 3. (a) Table of calculated eigenvalues and their relative errors computed by the current method for example 3; (b). List of eigenvalues with their relative errors attained by present method and FEM for Example 3.
Table 4. List of eigenvalues with their relative errors attained by present method and FEM for Example 4.
Table 5. Comparison of natural frequencies and relative errors using Gal. WRM with FEM for Example 5.
Figure 1. Comparison of relative errors obtained for Gal. MWR and FEM.
Figure 2. comparison of eigenvalues using GAL MWR and exact result.
Figure 3. Relative errors attained by GAL MWR for increasing polynomial.
For the vibration of crankshaft computed results utilizing Galerkin MWR are very close to the analytical results  and error increased as the modulation depth closer to one. In the case of 1-D heat equation for the plane wall  and 2-D heat equation , the smallest errors by present method attains superior accuracy than finite element methods using linear and quadrilateral elements worked out in those studies. Eigenvalues achieved from 1D Helmholtz equation for smaller eigenvalues converge faster as the degree of Bernstein polynomials increases and attains many accurate results. For homogeneous rectangular membrane, relative errors in the present method are smaller than that of finite element method . This proves that our present method is much more efficient than various numerical methods available in the literature and well suited which has much applicability in physical and engineering models.
The first author is very much thankful to Muhammad Sajjad Hossain, Assistant Professor of Mathematics, Department of Arts & Sciences, Ahsanullah University of Science and Technology, Dhaka, for his knowledgeable input, inexorable assistance and for always being keen to help in editing this article.
 Liu J., Zheng, Z. and Xu, Q. (2012) Bernstein-Polynomials-Based Highly Accurate Methods for One-Dimensional Interface Problems. Journal of Applied Mathematics, 2012, 1-11.
 Reutskiy, S.Y. (2006) The Method of Fundamental Solutions for Helmholtz Eigenvalues Problems in Simply and Multiply Connected Domains. Engineering Analysis with Boundary Elements, 30, 150-159.
 Yousefi, S.A., Behroozifar, M. and Dehghan, M. (2011) The Operational Matrices of Bernstein Polynomials for Solving the Parabolic Equation Subject to Specification of Mass. Journal of Computational and Applied Mathematics, 235, 5272-5283.
 Zhou-Lian, Z., Chang-Jiang L.H., Xiao, T. and Chen, S-L. (2009) Free Vibration Analysis of Rectangular Orthotropic Membranes in Large Deflection. Mathematical Problems in Engineering, 2009, 1-9.
 Mirkov, N. and Rasuo, B. (2013) Bernstein Polynomial Collocation Method for Elliptic Boundary Value Problems. Proceedings in Applied Mathematics and Mechanics, 13, 421-422.
 Yousefi, S.A. and Barikbiz, M.Z. (2010) Bernstein Ritz-Galerkin Method for Solving a Initial Boundary Value Problem that Combines Neumann and Integral Condition for the Wave Equation. Numerical Methods for Partial Differential Equations, 26, 1236-1246.
 Park, J., Park, I. and Lee, U. (2014) Transverse Vibration and Waves in a Membrane: Frequency Do Main Spectral Element Modeling and Analysis. Mathematical Problems in Engineering, 2014, 1-14.
 Sato, K. (1974) Free Vibration Analysis of a Composite Elliptical Membrane Consisting of Confocal Elliptical Parts. Journal of Sound and Vibration, 34, 161-171.
 Chen, S.-L. and Zheng, Z.-L. (2003) Large Deformation of Circular Membrane under the Concentrated Force. Applied Mathematics and Mechanics (English Edition), 24, 28-31.
 Kim, Y.Y. and Kang, J.H. (1996) Free Vibration Analysis of Membranes Using Wave-Type Base Functions. The Journal of the Acoustical Society of America, 99, 2938-2946.
 Laura, P.A.A., Rossi, R.E. and Gutierrez, R.H. (1997) The Fundamental Frequency of Non-Homogeneous Rectangular Membrane. Journal of Sound and Vibration, 204, 229-306.
 Vega, D.A., Vera, S.A. and Laura, P.A.A. (1999) Fundamental Frequency of Vibration of Rectangular Membranes with an Internal Oblique Support. Journal of Sound and Vibration, 224, 780-783.
 Ho, S.H. and Chen, C.K. (2000) Free Vibration Analysis of Non-Homogeneous Rectangular Membranes Using a Hybrid Methods. Journal of Sound and Vibration, 233, 547-555.
 Ersoy, H., Civalekb, O. and Ozpolatb, L. (2010) Free Vibration Analysis of Rectangular Membranes with Variable Density Using the Discrete Singular Convolution Approach. Asian Journal of Civil Engineering, 11, 83-94.
 Reutskiy, S.Y. (2007) The Methods of External and Internal Excitation for Problems of Free Vibrations of Non-Homogeneous Membranes. Engineering Analysis with Boundary Elements, 31, 906-918.
 Doha, E.H. and Abd-Elhameed, W.M. (2005) Accurate Spectral Solutions for the Parabolic and Elliptic Partial Differential Equations by the Ultra-Spherical Tau Method. Journal of Computational and Applied Mathematics , 181, 24-45.
 Sahuck, O.H. (2019) An Efficient Spectral Method to Solve Multi-Dimensional Linear Partial Differential Equations Using Chebyshev Polynomials. Mathematics, 7, 1-21.
 Dehghan, M. (2006) A Computational Study of the One-Dimensional Parabolic Equation Subject to Non-Classical Boundary Specifications. Numerical Methods for Partial Differential Equations, 22, 220-257.
 Alshbool, M.H.T., Bataineh, A.S., Hashim, I. and Isik, O.R. (2015) Approximate Solutions of Singular Differential Equations with Estimation Error by Using Bernstein Polynomials. International Journal of Pure and Applied Mathematics, 100, 109-125.
 Zhou, J. and Li, H. (2015) A Ritz-Galerkin Approximation to the Solution of Parabolic Equation with Moving Boundaries. Boundary Value Problems, 236, 1-17.
 Doha, E.H., Bhrawy, A.H. and Saker, M.A. (2011) On the Derivatives of Bernstein Polynomials, an Application for the Solution of High Even-Order Differential Equations. Boundary Value Problems, 2011, Article ID: 829543.
 Farzana, H., Islam, M.S. and Bhowmik, S.K. (2015) Computation of Eigenvalues of the Fourth Order Sturm-Liouville BVP by Galerkin Weighted Residual Method. British Journal of Mathematics and Computer Science, 9, 73-85.
 Islam, M.S., Farzana, H. and Bhowmik, S.K. (2017) Numerical Solutions of Sixth Order Eigenvalue Problems Using Galerkin Weighted Residual Method. Differential Equations and Dynamical Systems, 25, 187-205.
 Farzana, H. and Bhowmik, S.K. (2019) Comparative Study on Sixth Order Boundary Value Problems with Application to Linear Hydrodynamic Stability Problem and Benard Layer Eigenvalue Problem. Differential Equations and Dynamical Systems.