magnetic field and rotation. However when e * 0 the inhibiting effect of both magnetic field and rotation increases as the parameter c increases. An energy analysis of the Benard problem is considered when the magnetic field is in the direction of the vertical. The stabilizing effect of the magnetic field is investigated via a generalized energy functional for a non-linear magnetohydrodynamic fluid. It is shown that the magnetic field has a stabilizing effect for this model. The Compound Matrix Technique and the Chebyshev polynomials method have been used to solve the related eigenvalue problems. These methods shall be described in detail in chapter three and four respectively-Chapter Two Minimisation / Maximisation Procedures The Compound Matrix Method Let us examine the problem of determining an eigenvalue X of the nth order system-A (X , x) Y (3.1) subject to the boundary conditions B Y-0 on x-0 , (3.2) CY-0 on x-1 (3.3) where Y is an n vector, A(X,x) is an nxn matrix and B and C are respectively (n-m)xn and mxn matrices of full rank. i.e. (n-m) conditions are given at x-0 and m conditions at x = 1. We may assume without loss of generality that m < n/2. Notionally the general solution of the nth order system (3.1) has n degrees of freedom. Clearly (n-m) of these degrees of freedom can be removed by requiring that conditions (3.2) be satisfied and so we may conceptually construct m functionally independent solutions of (3.1) which also satisfy conditions (3.2). Let these solutions be w-i (X ,x) , w 2(X,x)........ wm (X,x) where n / ' (n" 0 n T Wi(X,x) = (wif w 1(W i ....... , w I) L Let us define the nxm matrix M to be the matrix whose rth column is w r (X,x) i.e. M = [w, , w 2, * 1 3-(2,3,6), d> 1 4 = (2,4,5), *1 5 = (2,4,6), $ 6 = (2,5,6),