Symmetric methods such as the implicit midpoint rule (IMR), implicit trapezoidal
rule (ITR) and 2-stage Gauss method are beneficial in solving Hamiltonian problems
since they are also symplectic. Symplectic methods have advantages over non-symplectic
methods in the long term integration of Hamiltonian problems. The study is to show
the efficiency of IMR, ITR and the 2-stage Gauss method in solving simple harmonic
oscillators (SHO). This study is done theoretically and numerically on the simple harmonic
oscillator problem. The theoretical analysis and numerical results on SHO problem
showed that the magnitude of the global error for a symmetric or symplectic method
with stepsize h is linearly dependent on time t. This gives the linear error growth when
a symmetric or symplectic method is applied to the simple harmonic oscillator problem.
Passive and active extrapolations have been implemented to improve the accuracy of the
numerical solutions. Passive extrapolation is observed to show quadratic error growth
after a very short period of time. On the other hand, active extrapolation is observed to
show linear error growth for a much longer period of time.
Analyzed the effects of thermal radiation, chemical reaction, heat gener-
ation/absorption, magnetic and electric fields on unsteady flow and heat transfer of
nanofluid. The transport equations used passively controlled. A similarity solution is
employed to transformed the governing equations from partial differential equations to
a set of ordinary differential equations, and then solve using Keller box method. It was
found that the temperature is a decreasing function with the thermal stratification due to
the fact the density of the fluid in the lower vicinity is much higher compared to the upper
region, whereas the thermal radiation, viscous dissipation and heat generation enhanced
the nanofluid temperature and thermal layer thickness.
Rainfall is an interesting phenomenon to investigate since it is directly related
to all aspects of life on earth. One of the important studies is to investigate and under-
stand the rainfall patterns that occur throughout the year. To identify the pattern, it
requires a rainfall curve to represent daily observation of rainfall received during the year.
Functional data analysis methods are capable to convert discrete data intoa function that
can represent the rainfall curve and as a result, try to describe the hidden patterns of the
rainfall. This study focused on the distribution of daily rainfall amount using functional
data analysis. Fourier basis functions are used for periodic rainfall data. Generalized
cross-validation showed 123 basis functions were sufficient to describe the pattern of daily
rainfall amount. North and west areas of the peninsula show a significant bimodal pattern
with the curve decline between two peaks at the mid-year. Meanwhile,the east shows uni-
modal patterns that reached a peak in the last three months. Southern areas show more
uniform trends throughout the year. Finally, the functional spatial method is introduced
to overcome the problem of estimating the rainfall curve in the locations with no data
recorded. We use a leave one out cross-validation as a verification method to compare
between the real curve and the predicted curve. We used coefficient of basis functions
to get the predicted curve. It was foundthatthe methods ofspatial prediction can match
up with theexistingspatialpredictionmethodsin terms of accuracy,but it isbetterasthe new
approach provides a simpler calculation.
Price stability is one of the main policy objectives that is targeted by policymakers in many countries. Price uncertainty occurs due to the changes in market structure and consumer preference and expectation, which may affect price stability. In this study, the researchers aimed to examine the effects of price uncertainty of consumer price disaggregation in Malaysian sectors. To be specific, the researchers were seeking to discover on how domestic and global commodity prices can affect sectoral Consumer Price Index (CPI) on price inflation in Malaysia and most importantly, whether the effect is different for economic sectors in Malaysia. In addition, the effects of other factors (i.e., internal and external factors) on Malaysian sectoral CPI inflation were also studied. The threshold generalized autoregressive conditional heteroscedasticity (TGARCH) model was used to generate the price uncertainties. For the purpose of analysis, the threshold regression approach was applied based on time series of each single sector, followed by a combination of panel data of all sectors. The results differed across sectors, revealing that the impact of price uncertainties was determined by the sensitivity of each sector towards the price uncertainties. The effect of price increase is larger than the effect of price decrease. Price fluctuations were obvious in sectors that were dependent on consumer price or commodity price. Exchange rate and oil price inflation had also greatly influenced the CPI inflation.
Logging activity is one of the most important activities for tropical countries
including Malaysia, as it produces quality trees for papers. One of the important tree
species is the Acacia Mangium which it produces a soft tree for papermaking enterprises.
The papers are exported to Europe and countries which have high demand for paper
due to the rapid development of the printing industry. Thus we analyzed the height for
individual trees. We investigate the maximum height of the trees from 1990 to 2006
and we fit the data using extreme value model. Some of the data are missing and three
imputation methods we used to solve this problem.
In this paper, extended Runge-Kutta fourth order method for directly solving the fuzzy logistic problem is presented. The extended Runge-Kutta method has lower number of function evaluations, compared with the classical Runge-Kutta method. The numerical robustness of the method in parameter estimation is enhanced via error minimization in predicting growth rate and carrying capacity. The results of fuzzy logistic model with the estimated parameters have been compared with population growth data in Malaysia, which indicate that this method is more accurate that the data population. Numerical example is given to illustrate the efficiency of the proposed model. It is concluded that robust parameter estimation technique is efficient in modelling population growth.
In this paper, we consider the system of Volterra-Fredholm integral equations
of the second kind (SVFI-2). We proposed fixed point method (FPM) to solve
SVFI-2 and improved fixed point method (IFPM) for solving the problem. In addition,
a few theorems and two new algorithms are introduced. They are supported by
numerical examples and simulations using Matlab. The results are reasonably good
when compared with the exact solutions.
In this paper, the problem of forced convection flow of micropolar fluid of
lighter density impinging orthogonally on another heavier density of micropolar fluid
on a stretching surface is investigated. The boundary layer governing equations are
transformed from partial differential equations into a system of nonlinear ordinary
differential equations using similarity transformation and solved numerically using dsolve
function in maple software version 2016. The velocity, microrotation and temperature of
micropolar fluid are analyzed. It is found that both upper fluid and lower fluid display
opposite behaviour when micropolar parameter k various with strong concentration
n = 0, pr = 7 and stretching parameter = 0.5. The results also show that stretching
surface exert the force that increasing the velocity of micropolar fluid.
An accurate forecasting of tropospheric ozone (O3) concentration is benefi-
cial for strategic planning of air quality. In this study, various forecasting techniques are
used to forecast the daily maximum O3 concentration levels at a monitoring station
in the Klang Valley, Malaysia. The Box-Jenkins autoregressive integrated movingaverage
(ARIMA) approach and three types of neural network models, namely, backpropagation
neural network, Elman recurrent neural network and radial basis function
neural network are considered. The daily maximum data, spanning from 1 January
2011 to 7 August 2011, was obtained from the Department of Environment, Malaysia.
The performance of the four methods in forecasting future values of ozone concentrations
is evaluated based on three criteria, which are root mean square error (RMSE),
mean absolute error (MAE) and mean absolute percentage error (MAPE). The findings
show that the Box-Jenkins approach outperformed the artificial neural network
methods.
A domination polynomial is a type of graph polynomial in which its coefficients represent the number of dominating sets in the graph. There are many researches being done on the domination polynomial of some common types of graphs but not yet for graphs associated to finite groups. Two types of graphs associated to finite groups are the conjugate graph and the conjugacy class graph. A graph of a group G is called a conjugate graph if the vertices are non-central elements of G and two distinct vertices are adjacent if they are conjugate to each other. Meanwhile, a conjugacy class graph of a group G is a graph in which its vertices are the non-central conjugacy classes of G and two distinct vertices are connected if and only if their class cardinalities are not coprime. The conjugate and conjugacy class graph of dihedral groups can be expressed generally as a union of complete graphs on some vertices. In this paper, the domination polynomials are computed for the conjugate and conjugacy class graphs of the dihedral groups.
In DNA splicing system, the potential effect of sets of restriction enzymes and
a ligase that allow DNA molecules to be cleaved and re-associated to produce further
molecules is modelled mathematically. This modelling is done in the framework of formal
language theory, in which the nitrogen bases, nucleotides and restriction sites are modelled
as alphabets, strings and rules respectively. The molecules resulting from a splicing system
is depicted as the splicing language. In this research, the splicing language resulting from
DNA splicing systems with one palindromic restriction enzyme for one and two (nonoverlapping)
cutting sites are generalised as regular expressions.
This study presents a two-strain deterministic model which incorporates Dengvaxia vaccine and insecticide (adulticide) control strategies to forecast the dynamics of transmission and control of dengue in Madeira Island if there is a new outbreak with a different virus serotypes after the first outbreak in 2012. We construct suitable Lyapunov functions to investigate the global stability of the disease-free and boundary equilibrium points. Qualitative analysis of the model which incorporates time-varying controls with the specific goal of minimizing dengue disease transmission and the costs related to the control implementation by employing the optimal control theory is carried out. Three strategies, namely the use of Dengvaxia vaccine only, application of adulticide only, and the combination of Dengvaxia vaccine and adulticide are considered for the controls implementation. The necessary conditions are derived for the optimal control of dengue. We examine the impacts of the control strategies on the dynamics of infected humans and mosquito population by simulating the optimality system. The disease-free equilibrium is found to be globally asymptotically stable whenever the basic reproduction numbers associated with virus serotypes 1 and j (j ∈ {2,3,4}), respectively, satisfy R01,R0j ≤ 1, and the boundary equilibrium is globally asymptotically stable when the related R0i (i = 1,j) is above one. It is shown that the strategy based on the combination of Dengvaxia vaccine and adulticide helps in an effective control of dengue spread in the Island.
Blood flow through a bifurcated artery with the presence of an overlapping stenosis located at parent’s arterial lumen under the action of a uniform external magnetic field is studied in this paper. Blood is treated as an electrically conducting fluid which exhibits the Magnetohydrodynamics principle and it is characterized by a Newtonian fluid model. The governing equations are discretized using a stabilization technique of finite element known as Galerkin least-squares. The maximum velocity and pressure drop evaluated in this present study are compared with the results found in previous literature and COMSOL Multiphysics. The solutions found in a satisfactory agreement, thus verify the source code is working properly. The effects of dimensionless parameters of Hartmann and Reynolds numbers in the fluid’s velocity and pressure are examined in details with further scientific discussions.
Heat and mass transfer of MHD boundary-layer flow of a viscous incompress-
ible fluid over an exponentially stretching sheet in the presence of radiation is investi-
gated. The two-dimensional boundary-layer governing partial differential equations are
transformed into a system of nonlinear ordinary differential equations by using similarity
variables. The transformed equations of momentum, energy and concentration are solved
by Homotopy Analysis Method (HAM). The validity of HAM solution is ensured by com-
paring the HAM solution with existing solutions. The influence of physical parameters
such as magnetic parameter, Prandtl number, radiation parameter, and Schmidt num-
ber on velocity, temperature and concentration profiles are discussed. It is found that
the increasing values of magnetic parameter reduces the dimensionless velocity field but
enhances the dimensionless temperature and concentration field. The temperature dis-
tribution decreases with increasing values of Prandtl number. However, the temperature
distribution increases when radiation parameter increases. The concentration boundary
layer thickness decreases as a result of increase in Schmidt number.
Cardiovascular disease (CVD) includes coronary heart disease, cerebrovascular disease (stroke), peripheral artery disease, and atherosclerosis of the aorta. All females face the threat of CVD. But becoming aware of symptoms and signs is a great challenge since most adults at increased risk of cardiovascular disease (CVD) have no symptoms or obvious signs especially in females. The symptoms may be identified by the assessment of their risk factors. The Bayesian approach is a specific way in dealing with this kind of problem by formalizing a priori beliefs and of combining them with the available observations. This study aimed to identify associated risk factors in CVD among female patients presenting with ST Elevation Myocardial Infarction (STEMI) using Bayesian logistic regression and obtain a feasible model to describe the data. A total of 874 STEMI female patients in the National Cardiovascular Disease Database-Acute Coronary Syndrome (NCVD-ACS) registry year 2006-2013 were analysed. Bayesian Markov Chain Monte Carlo (MCMC) simulation approach was applied in the univariate and multivariate analysis. Model performance was assessed through the model calibration and discrimination. The final multivariate model of STEMI female patients consisted of six significant variables namely smoking, dyslipidaemia, myocardial infarction (MI), renal disease, Killip class and age group. Females aged 65 years and above have higher incidence of CVD and mortality is high among female patients with Killip class IV. Also, renal disease was a strong predictor of CVD mortality. Besides, performance measures for the model was considered good. Bayesian logistic regression model provided a better understanding on the associated risk factors of CVD for female patients which may help tailor prevention or treatment plans more effectively.
Recent studies have shown that independent identical distributed Gaussian
random variables is not suitable for modelling extreme values observed during extremal
events. However, many real life data on extreme values are dependent and stationary
rather than the conventional independent identically distributed data. We propose a stationary
autoregressive (AR) process with Gumbel distributed innovation and characterise
the short-term dependence among maxima of an (AR) process over a range of sample
sizes with varying degrees of dependence. We estimate the maximum likelihood of the
parameters of the Gumbel AR process and its residuals, and evaluate the performance
of the parameter estimates. The AR process is fitted to the Gumbel-generalised Pareto
(GPD) distribution and we evaluate the performance of the parameter estimates fitted
to the cluster maxima and the original series. Ignoring the effect of dependence leads to
overestimation of the location parameter of the Gumbel-AR (1) process. The estimate
of the location parameter of the AR process using the residuals gives a better estimate.
Estimate of the scale parameter perform marginally better for the original series than the
residual estimate. The degree of clustering increases as dependence is enhance for the AR
process. The Gumbel-AR(1) fitted to the threshold exceedances shows that the estimates
of the scale and shape parameters fitted to the cluster maxima perform better as sample
size increases, however, ignoring the effect of dependence lead to an underestimation of
the parameter estimates of the scale parameter. The shape parameter of the original
series gives a superior estimate compare to the threshold excesses fitted to the Gumbel
distributed Generalised Pareto ditribution.
Hantaviruses are etiological agents of zoonotic diseases and certain other dis-
eases, which pose a serious threat to human health. When rodent and predator popula-
tions share in an ecology, the competitive force of the populations can lead to a reduction
or elimination of a hantavirus outbreak. The effect of the predator eliminating rodents
and predator populations that tends to reduce or eliminate hantavirus infection is investi-
gated. The existence of several equilibrium points of the model is identified and local and
global stabilities of the model at these equilibrium points are analysed in detail. Numerical
simulations are carried out to illustrate our model results.
Inventory Routing Problem (IRP) has been continuously developed and improved due to pressure from global warming issue particularly related to greenhouse gases (GHGs) emission. The burning of fossil fuel for transportations such as cars, trucks, ships, trains, and planes primarily emits GHGs. Carbon dioxide (CO2) from burning of fossil fuel to power transportation and industrial process is the largest contributor to global GHGs emission. Therefore, the focus of this study is on solving a multi-period inventory routing problem (MIRP) involving carbon emission consideration based on carbon cap and offset policy. Hybrid genetic algorithm (HGA) based on allocation first and routing second is used to compute a solution for the MIRP in this study. The objective of this study is to solve the proposed MIRP model with HGA then validate the effectiveness of the proposed HGA on data of different sizes. Upon validation, the proposed MIRP model and HGA is applied on real-world data. The HGA is found to be able to solve small size and large size instances effectively by providing near optimal solution in relatively short CPU execution time.
Subsea cable laying process is a difficult task for an engineer due to many
uncertain situations which occur during the operation. It is very often that the cable being
laid out is not perfectly fit on the route being planned, which results in the formation of
slack. In order to control wastages during installation, the slack needs to be minimized
and the movement of a ship/vessel needs to be synchronized with the cable being laid out.
The current problem was addressed using a mathematical model by considering a number
of defining parameters such as the external forces, the cable properties and geometry. Due
to the complexity, the model is developed for a steady-state problem assuming velocity
of the vessel is constant, seabed is flat and the effect of wind and wave is insignificant.
Non-dimensional system is used to scale the engineering parameters and grouped them
into only two main parameters which are the hydrodynamic drag of the fluid and the
bending stiffness of the cable. There are two solutions generated in this article; numerical
and asymptotic solutions. The result of these solutions suggests that the percentage of
slack can be reduced by the increase of the prescribed cable tension, and also the increase
in either the drag coefficient of the sea water or the bending stiffness of the cable, similarly
will result in lower slack percentage
In irradiation process, instead of traverse on the targeted cells, there is side
effect happens to non-targeted cells. The targeted cells that had been irradiated with
ionizing radiation emits damaging signal molecules to the surrounding and then, dam-
age the bystander cells. The type of damage considered in this work is the number of
double-strand breaks (DSBs) of deoxyribonucleic acid (DNA) in cell’s nucleus. By us-
ing mathematical approach, a mechanistic model that can describe this phenomenon is
developed based on a structured population approach. Then, the accuracy of the model
is validated by its ability to match the experimental data. The Particle Swarm (PS)
optimization is employed for the data fitting procedure. PS optimization searches the
parameter value that minimize the errors between the model simulation data and exper-
imental data. It is obtained that the mathematical modelling proposed in this paper is
strongly in line with the experimental data.