International Journal of Computational Mathematics
http://www.hindawi.com
The latest articles from Hindawi Publishing Corporation
© 2014 , Hindawi Publishing Corporation . All rights reserved.

Computation of a Canonical Form for Linear 2D Systems
Mon, 15 Dec 2014 09:04:54 +0000
http://www.hindawi.com/journals/ijcm/2014/487465/
Symbolic computation techniques are used to obtain a canonical
form for polynomial matrices arising from discrete 2D linear statespace
systems. The canonical form can be regarded as an extension of the
companion form often encountered in the theory of 1D linear systems. Using
previous results obtained by Boudellioua and Quadrat (2010) on the reduction by
equivalence to Smith form, the exact connection between the original polynomial
matrix and the reduced canonical form is set out. An example is given to
illustrate the computational aspects involved.
Mohamed Salah Boudellioua
Copyright © 2014 Mohamed Salah Boudellioua. All rights reserved.

Computational Modelling of Couette Flow of Nanofluids with Viscous Heating and Convective Cooling
Sun, 14 Dec 2014 07:24:01 +0000
http://www.hindawi.com/journals/ijcm/2014/631749/
The combined effect of viscous heating and convective cooling on Couette flow and heat transfer characteristics of water base nanofluids containing Copper Oxide (CuO) and Alumina (Al2O3) as nanoparticles is investigated. It is assumed that the nanofluid flows in a channel between two parallel plates with the channel’s upper plate accelerating and exchange heat with the ambient surrounding following the Newton’s law of cooling, while the lower plate is stationary and maintained at a constant temperature. Using appropriate similarity transformation, the governing NavierStokes and the energy equations are reduced to a set of nonlinear ordinary differential equations. These equations are solved analytically by regular perturbation method with series improvement technique and numerically by an efficient RungeKuttaFehlberg integration technique coupled with shooting method. The effects of the governing parameters on the dimensionless velocity, temperature, skin friction, pressure drop and Nusselt number are presented graphically, and discussed quantitatively.
Oluwole Daniel Makinde, Ahmada Omar, and M. Samuel Tshehla
Copyright © 2014 Oluwole Daniel Makinde et al. All rights reserved.

Multiresolution Analysis Based on Coalescence HiddenVariable Fractal Interpolation Functions
Thu, 11 Dec 2014 06:39:48 +0000
http://www.hindawi.com/journals/ijcm/2014/531562/
Multiresolution analysis arising from Coalescence Hiddenvariable Fractal Interpolation Functions (CHFIFs) is developed. The availability of a larger set of free variables and constrained variables with CHFIF in multiresolution analysis based on CHFIFs provides more control in reconstruction of functions in than that provided by multiresolution analysis based only on Affine Fractal Interpolation Functions (AFIFs). Our approach consists of introduction of the vector space of CHFIFs, determination of its dimension and construction of Riesz bases of vector subspaces , consisting of certain CHFIFs in .
G. P. Kapoor and Srijanani Anurag Prasad
Copyright © 2014 G. P. Kapoor and Srijanani Anurag Prasad. All rights reserved.

Extended Vector Equilibrium Problem
Wed, 03 Dec 2014 12:04:33 +0000
http://www.hindawi.com/journals/ijcm/2014/643230/
We introduce and study extended vector equilibrium problem. By using KKMFan Theorem as basic tool, we prove existence theorem in the setting of Hausdorff topological vector space and reflexive Banach space. Some examples are also given.
Khushbu and Zubair Khan
Copyright © 2014 Khushbu and Zubair Khan. All rights reserved.

Combined Effect of Surface Roughness and Slip Velocity on Jenkins Model Based Magnetic Squeeze Film in Curved Rough Circular Plates
Wed, 03 Dec 2014 09:08:58 +0000
http://www.hindawi.com/journals/ijcm/2014/367618/
This paper aims to discuss the effect of slip velocity and surface roughness on the performance of Jenkins model based magnetic squeeze film in curved rough circular plates. The upper plate’s curvature parameter is governed by an exponential expression while a hyperbolic form describes the curvature of lower plates. The stochastic model of Christensen and Tonder has been adopted to study the effect of transverse surface roughness of the bearing surfaces. Beavers and Joseph’s slip model has been employed here. The associated Reynolds type equation is solved to obtain the pressure distribution culminating in the calculation of load carrying capacity. The computed results show that the Jenkins model modifies the performance of the bearing system as compared to NeuringerRosensweig model, but this model provides little support to the negatively skewed roughness for overcoming the adverse effect of standard deviation and slip velocity even if curvature parameters are suitably chosen. This study establishes that for any type of improvement in the performance characteristics the slip parameter is required to be reduced even if variance (−ve) occurs and suitable magnetic strength is in force.
Jimit R. Patel and Gunamani Deheri
Copyright © 2014 Jimit R. Patel and Gunamani Deheri. All rights reserved.

Matrices in Fuzzy Linear Systems
Sun, 30 Nov 2014 07:32:00 +0000
http://www.hindawi.com/journals/ijcm/2014/394135/
We consider a class of fuzzy linear system of equations and demonstrate some of the existing challenges. Furthermore, we explain the efficiency of this model when the coefficient matrix is an matrix. Numerical experiments are illustrated to show the applicability of the theoretical analysis.
H. Saberi Najafi and S. A. Edalatpanah
Copyright © 2014 H. Saberi Najafi and S. A. Edalatpanah. All rights reserved.

On the Dynamics of Laguerre’s Iteration Method for Finding the th Roots of Unity
Wed, 26 Nov 2014 11:44:16 +0000
http://www.hindawi.com/journals/ijcm/2014/321585/
Previous analyses of Laguerre’s iteration method have provided results on the behavior of this popular method when applied to the polynomials , . In this paper, we summarize known analytical results and provide new results. In particular, we study symmetry properties of the Laguerre iteration function and clarify the dynamics of the method. We show analytically and demonstrate computationally that for each the basin of attraction to the roots is a subset of an annulus that contains the unit circle and whose Lebesgue measure shrinks to zero as . We obtain a good estimate of the size of the bounding annulus. We show that the boundary of the basin of convergence exhibits fractal nature and quasi selfsimilarity. We also discuss the connectedness of the basin for large values of . We also numerically find some short finite cycles on the boundary of the basin of convergence for . Finally, we demonstrate that when using the floating point arithmetic and the general formulation of the method, convergence occurs even from starting values outside of the basin of convergence due to the loss of significance during the evaluation of the iteration function.
Pavel Bělík, HeeChan Kang, Andrew Walsh, and Emma Winegar
Copyright © 2014 Pavel Bělík et al. All rights reserved.

Plane Elastostatic Solution in an Infinite Functionally Graded Layer Weakened by a Crack Lying in the Middle of the Layer
Tue, 25 Nov 2014 16:57:14 +0000
http://www.hindawi.com/journals/ijcm/2014/358617/
This paper is concerned with an internal crack problem in an infinite functionally graded elastic layer. The crack is opened by an internal uniform pressure along its surface. The layer surfaces are supposed to be acted on by symmetrically applied concentrated forces of magnitude with respect to the centre of the crack. The applied concentrated force may be compressive or tensile in nature. Elastic parameters λ and μ are assumed to vary along the normal to the plane of crack. The problem is solved by using integral transform technique. The solution of the problem has been reduced to the solution of a Cauchytype singular integral equation, which requires numerical treatment. The stressintensity factors and the crack opening displacements are determined and the effects of graded parameters on them are shown graphically.
R. Patra, S. P. Barik, M. Kundu, and P. K. Chaudhuri
Copyright © 2014 R. Patra et al. All rights reserved.

A Numerical Method for 1D Parabolic Equation with Nonlocal Boundary Conditions
Thu, 20 Nov 2014 09:44:34 +0000
http://www.hindawi.com/journals/ijcm/2014/923693/
This paper is concerned with a local method for the solution of onedimensional parabolic equation with nonlocal boundary conditions. The method uses a coordinate transformation. After the coordinate transformation, it is then possible to obtain exact solutions for the resulting equations in terms of the local variables. These exact solutions are in terms of constants of integration that are unknown. By imposing the given boundary conditions and smoothness requirements for the solution, it is possible to furnish a set of linearly independent conditions that can be used to solve for the constants of integration. A number of examples are used to study the applicability of the method. In particular, three nonlinear problems are used to show the novelty of the method.
M. Tadi and Miloje Radenkovic
Copyright © 2014 M. Tadi and Miloje Radenkovic. All rights reserved.

A Numerical Test of Padé Approximation for Some Functions with Singularity
Thu, 20 Nov 2014 00:00:00 +0000
http://www.hindawi.com/journals/ijcm/2014/587430/
The aim of this study is to examine some numerical tests of Padé approximation for some typical functions with singularities such as simple pole, essential singularity, brunch cut, and natural boundary. As pointed out by Baker, it was shown that the simple pole and the essential singularity can be characterized by the poles of the Padé approximation. However, it was not fully clear how the Padé approximation works for the functions with the branch cut or the natural boundary. In the present paper, it is shown that the poles and zeros of the Padé approximated functions are alternately lined along the branch cut if the test function has branch cut, and poles are also distributed around the natural boundary for some lacunary power series and random power series which rigorously have a natural boundary on the unit circle. On the other hand, Froissart doublets due to numerical errors and/or external noise also appear around the unit circle in the Padé approximation. It is also shown that the residue calculus for the Padé approximated functions can be used to confirm the numerical accuracy of the Padé approximation and quasianalyticity of the random power series.
Hiroaki S. Yamada and Kensuke S. Ikeda
Copyright © 2014 Hiroaki S. Yamada and Kensuke S. Ikeda. All rights reserved.

Block Hybrid Step Backward Differentiation Formulas for Large Stiff Systems
Mon, 20 Oct 2014 07:21:02 +0000
http://www.hindawi.com/journals/ijcm/2014/162103/
This paper presents a generalized high order block hybrid step backward differentiation formula (HBDF) for solving stiff systems, including large systems resulting from the semidiscretization parabolic partial differential equations (PDEs). A block scheme in which two offgrid points are specified by the zeros of the second degree Chebyshev polynomial of the first kind is examined for convergence, and stabilities. Numerical simulations that illustrate the accuracy of a Chebyshev based method are given for selected stiff systems and partial differential equations.
S. N. Jator and E. Agyingi
Copyright © 2014 S. N. Jator and E. Agyingi. All rights reserved.

Modelling Hepatotoxicity of Antiretroviral Therapy in the Liver during HIV Monoinfection
Sun, 19 Oct 2014 00:00:00 +0000
http://www.hindawi.com/journals/ijcm/2014/659675/
Liver related complications are currently the leading cause of morbidity and mortality among human immunodeficiency virus (HIV) infected individuals. In HIV monoinfected individuals on therapy, liver injury has been associated with the use of antiretroviral agents as most of them exhibit some degree of toxicity. In this study we proposed a mathematical model with the aim of investigating hepatotoxicity of combinational therapy of antiretroviral drugs. Therapy efficacy and toxicity were incorporated in the model as doseresponse functions. With the parameter values used in the study, protease inhibitorsbased regimens were found to be more toxic than nonnucleoside reverse transcriptase inhibitorsbased regimens. In both regimens, the combination of stavudine and zidovudine was the most toxic baseline nucleoside reverse transcriptase inhibitors followed by didanosine with stavudine. However, the least toxic combinations were zidovudine and lamivudine followed by didanosine and lamivudine. The study proposed that, under the same second line regimens, the most toxic first line combination gives the highest viral load and vice versa.
Hasifa Nampala, Livingstone S. Luboobi, Joseph Y. T. Mugisha, and Celestino Obua
Copyright © 2014 Hasifa Nampala et al. All rights reserved.

Mathematical Modeling of Multienzyme Biosensor System
Wed, 01 Oct 2014 12:16:02 +0000
http://www.hindawi.com/journals/ijcm/2014/694037/
A mathematical model of hybrid inhibitor biosensor system is discussed. This model consists of five nonlinear partial differential equations for bisubstrate sensitive amperometric system. Simple and closed form of analytical expressions for concentration of glucose6phosphate (substrate), potassium dihydrogen phosphate (inhibitor), oxygen (cosubstrate), glucose (product 1), and hydrogen peroxide (product 3) is obtained in terms of rate constant using modified Adomian decomposition method (MADM). In this study, behavior of biokinetic parameters is analyzed using this theoretical result. The obtained analytical results (concentrations) are compared with the numerical results and are found to be in satisfactory agreement.
SP. Ganesan, K. Saravanakumar, and L. Rajendran
Copyright © 2014 SP. Ganesan et al. All rights reserved.

Approximate Periodic Solution for the Nonlinear HelmholtzDuffing Oscillator via Analytical Approaches
Mon, 29 Sep 2014 11:52:45 +0000
http://www.hindawi.com/journals/ijcm/2014/939623/
The conservative HelmholtzDuffing oscillator is analyzed by means of three analytical techniques. The maxmin, secondorder of the Hamiltonian, and the global error minimization approaches are applied to achieve natural frequencies. The obtained results are compared with the homotopy perturbation method and numerical solutions. The results show that secondorder of the global error minimization method is very accurate, so it can be widely applicable in engineering problems.
A. Mirzabeigy, M. K. Yazdi, and M. H. Nasehi
Copyright © 2014 A. Mirzabeigy et al. All rights reserved.

A Collocation Method for Numerical Solution of Hyperbolic Telegraph Equation with Neumann Boundary Conditions
Sun, 28 Sep 2014 11:46:26 +0000
http://www.hindawi.com/journals/ijcm/2014/526814/
We present a technique based on collocation of cubic Bspline basis functions to solve second order onedimensional hyperbolic telegraph equation with Neumann boundary conditions. The use of cubic Bspline basis functions for spatial variable and its derivatives reduces the problem into system of first order ordinary differential equations. The resulting system subsequently has been solved by SSPRK54 scheme. The accuracy of the proposed approach has been confirmed with numerical experiments, which shows that the results obtained are acceptable and in good agreement with the exact solution.
R. C. Mittal and Rachna Bhatia
Copyright © 2014 R. C. Mittal and Rachna Bhatia. All rights reserved.

Hermitian Positive Definite Solution of the Matrix Equation
Wed, 24 Sep 2014 08:25:11 +0000
http://www.hindawi.com/journals/ijcm/2014/727093/
We consider the Hermitian positive definite solution of the nonlinear matrix equation . Some new sufficient conditions and necessary conditions for the existence of Hermitian positive definite solutions are derived. An iterative method is proposed to compute the Hermitian positive definite solution. In the end, an example is used to illustrate the correctness and application of our results.
ChunMei Li and JingJing Peng
Copyright © 2014 ChunMei Li and JingJing Peng. All rights reserved.

Solving Operator Equation Based on Expansion Approach
Sun, 07 Sep 2014 07:08:56 +0000
http://www.hindawi.com/journals/ijcm/2014/671965/
To date, researchers usually use spectral and pseudospectral methods for only numerical approximation of ordinary and partial differential equations and also based on polynomial basis. But the principal importance of this paper is to develop the expansion approach based on general basis functions (in particular case polynomial basis) for solving general operator equations, wherein the particular cases of our development are integral equations, ordinary differential equations, difference equations, partial differential equations, and fractional differential equations. In other words, this paper presents the expansion approach for solving general operator equations in the form , with respect to boundary condition , where , and are linear, nonlinear, and boundary operators, respectively, related to a suitable Hilbert space, is the domain of approximation, is an arbitrary constant, and is an arbitrary function. Also the other importance of this paper is to introduce the general version of pseudospectral method based on general interpolation problem. Finally some experiments show the accuracy of our development and the error analysis is presented in norm.
A. Aminataei, S. AhmadiAsl, and M. Pakbaz
Copyright © 2014 A. Aminataei et al. All rights reserved.

On the Spectrum and Spectral Norms of Circulant Matrices with Generalized Horadam Numbers Entries
Sun, 31 Aug 2014 10:56:35 +0000
http://www.hindawi.com/journals/ijcm/2014/795175/
This work is concerned with the spectrum and spectral norms of circulant matrices with generalized Horadam numbers entries. By using Abel transformation and some identities we obtain an explicit formula for the eigenvalues of them. In addition, a sufficient condition for an circulant matrix to be normal is presented. Based on the results we obtain the precise value for spectral norms of normal circulant matrix with generalized Horadam numbers, which generalize and improve the known results.
Lele Liu
Copyright © 2014 Lele Liu. All rights reserved.

Solving the Generalized Regularized Long Wave Equation Using a Distributed Approximating Functional Method
Tue, 12 Aug 2014 12:53:08 +0000
http://www.hindawi.com/journals/ijcm/2014/178024/
The generalized regularized long wave (GRLW) equation is solved numerically by using a distributed approximating functional (DAF) method realized by the regularized Hermite local spectral kernel. Test problems including propagation of single solitons, interaction of two and three solitons, and conservation properties of mass, energy, and momentum of the GRLW equation are discussed to test the efficiency and accuracy of the method. Furthermore, using the Maxwellian initial condition, we show that the number of solitons which are generated can be approximately determined. Comparisons are made between the results of the proposed method, analytical solutions, and numerical methods. It is found that the method under consideration is a viable alternative to existing numerical methods.
Edson Pindza and Eben Maré
Copyright © 2014 Edson Pindza and Eben Maré. All rights reserved.

The Convergence of Geometric Mesh Cubic Spline Finite Difference Scheme for Nonlinear Higher Order TwoPoint Boundary Value Problems
Wed, 23 Jul 2014 09:30:37 +0000
http://www.hindawi.com/journals/ijcm/2014/527924/
An efficient algorithm for the numerical solution of higher (even) orders twopoint nonlinear boundary value problems has been developed. The method is third order accurate and applicable to both singular and nonsingular cases. We have used cubic spline polynomial basis and geometric mesh finite difference technique for the generation of this new scheme. The irreducibility and monotone property of the iteration matrix have been established and the convergence analysis of the proposed method has been discussed. Some numerical experiments have been carried out to demonstrate the computational efficiency in terms of convergence order, maximum absolute errors, and root mean square errors. The numerical results justify the reliability and efficiency of the method in terms of both order and accuracy.
Navnit Jha, R. K. Mohanty, and Vinod Chauhan
Copyright © 2014 Navnit Jha et al. All rights reserved.