Real life phenomena found in various fields such as engineering, physics,
biology and communication theory can be modeled as nonlinear higher order ordinary
differential equations, particularly the Duffing oscillator. Analytical solutions for these
differential equations can be time consuming whereas, conventional numerical solutions
may lack accuracy. This research propose a block multistep method integrated with a
variable order step size (VOS) algorithm for solving these Duffing oscillators directly.
The proposed VOS Block method provides an alternative numerical solution by reducing
computational cost (time) but without loss of accuracy. Numerical simulations
are compared with known exact solutions for proof of accuracy and against current
numerical methods for proof of efficiency (steps taken).
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, we study the numerical method for solving second order Fuzzy
Differential Equations (FDEs) using Block Backward Differential Formulas (BBDF)
under generalized concept of higher-order fuzzy differentiability. Implementation of
the method using Newton iteration is discussed. Numerical results obtained by BBDF
are presented and compared with Backward Differential Formulas (BDF) and exact
solutions. Several numerical examples are provided to illustrate our methods.
A new method to construct the distinct Hamiltonian circuits in complete
graphs is called Half Butterfly Method. The Half Butterfly Method used the concept
of isomorphism in developing the distinct Hamiltonian circuits. Thus some theoretical
works are presented throughout developing this method.
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.
In this paper, the application of the method of lines (MOL) to the Forced
Korteweg-de Vries-Burgers equation with variable coefficient (FKdVB) is presented.
The MOL is a powerful technique for solving partial differential equations by typically
using finite-difference approximations for the spatial derivatives and ordinary differential
equations (ODEs) for the time derivative. The MOL approach of the FKdVB
equation leads to a system of ODEs. The solution of the system of ODEs is obtained
by applying the Fourth-Order Runge-Kutta (RK4) method. The numerical solution
obtained is then compared with its progressive wave solution in order to show the
accuracy of the MOL method.
In this paper, we propose a method how to manage the convergence of
Newton’s method if its iteration process encounters a local extremum. This idea establishes
the osculating circle at a local extremum. It then uses the radius of the
osculating circle also known as the radius of curvature as an additional number of
the local extremum. It then takes that additional number and combines it with the
local extremum. This is then used as an initial guess in finding a root near to that
local extremum. This paper will provide several examples which demonstrate that the
proposed idea is successful and they perform to fulfill the aim of this paper.
Homotopy continuation methods (HCMs) can be used to find the solutions
of polynomial equations. The advantages of HCMs over classical methods such as the
Newton and bisection methods are that HCMs are able to resolve divergence and starting
value problems. In this paper, we develop Super Ostrowski-HCM as a technique to
overcome the starting value problem. We compare the performance of this proposed
method with Ostrowski-HCM. The results provide evidence of the superiority of Super
Ostrowski-HCM.
In this paper we consider a harvesting model of predator-prey fishery in which
the prey is directly infected by some external toxic substances. The toxic infection is
indirectly transmitted to the predator during the feeding process. The model is a modified
version from the classic Lotka-Volterra predator-prey model. The stability and bifurcation
analyses are addressed. Numerical simulations of the model are performed and bifurcation
diagrams are studied to investigate the dynamical behaviours between the predator and
the prey. The effects of toxicity and harvesting on the stability of steady states found in
the model are discussed.
In this paper, Maxwell fluid over a flat plate for convective boundary layer
flow with pressure gradient parameter is considered. The aim of this study is to compare
and analyze the effects of the presence and absence of λ (relaxation time), and also the
effects of m (pressure gradient parameter) and Pr (Prandtl number)on the momentum
and thermal boundary layer thicknesses. An approximation technique namely Homotopy
Perturbation Method (HPM) has been used with an implementation of Adam and Gear
Method’s algorithms. The obtained results have been compared for zero relaxation time
and also pressure gradient parameter with the published work of Fathizadeh and Rashidi.
The current outcomes are found to be in good agreement with the published results.
Physical interpretations have been given for the effects of the m, Pr and β (Deborah
number) with λ. This study will play an important role in industrial and engineering
applications.
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.
Numerical simulation of the behaviour of blood flow through a stenosed bifurcated
artery with the presence of single mild stenosis at parent artery is investigated. The
flow analysis applies the incompressible, steady, three-dimensional Navier-Stokes equations
for non-Newtonian generalized power law fluids. Behaviour of blood flow is simulated
numerically using COMSOL Multiphysicsthat based on finite element method.The
results showthe effect of severity of stenosis on flow characteristics such as axial velocity
and its exhibit flow recirculation zone for analysis on streamlines pattern.
Monthly data about oil production at several drilling wells is an example of
spatio-temporal data. The aim of this research is to propose nonlinear spatio-temporal
model, i.e. Feedforward Neural Network - VectorAutoregressive (FFNN-VAR) and FFNN
- Generalized Space-Time Autoregressive (FFNN-GSTAR), and compare their forecast
accuracy to linearspatio-temporal model, i.e. VAR and GSTAR. These spatio-temporal
models are proposed and applied for forecasting monthly oil production data at three
drilling wells in East Java, Indonesia. There are 60 observations that be divided to two
parts, i.e. the first 50 observations for training data and the last 10 observations for
testing data. The results show that FFNN-GSTAR(11) and FFNN-VAR(1) as nonlinear
spatio-temporal models tend to give more accurate forecast than VAR(1) and GSTAR(11)
as linear spatio-temporal models. Moreover, further research about nonlinear spatiotemporal
models based on neural networks and GSTAR is needed for developing new
hybrid models that could improve the forecast accuracy.
Simulation is used to measure the robustness and the efficiency of the forecasting
techniques performance over complex systems. A method for simulating multivariate
time series was presented in this study using vector autoregressive base-process. By
applying the methodology to the multivariable meteorological time series, a simulation
study was carried out to check for the model performance. MAPE and MAE performance
measurements were used and the results show that the proposed method that consider
persistency in volatility gives better performance and the accuracy error is six time smaller
than the normal hybrid model.
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.
The selection criteria play an important role in the portfolio optimization
using any ratio model. In this paper, the authors have considered the mean return as
profit and variance of return as risk on the asset return as selection criteria, as the first
stage to optimize the selected portfolio. Furthermore, the sharp ratio (SR) has been
considered to be the optimization ratio model. In this regard, the historical data taken
from Shanghai Stock Exchange (SSE) has been considered. A metaheuristic technique
has been developed, with financial tool box available in MATLAB and the particle swarm
optimization (PSO) algorithm. Hence, called as the hybrid particle swarm optimization
(HPSO) or can also be called as financial tool box particle swarm optimization (FTBPSO).
In this model, the budgets as constraint, where as two different models i.e. with
and without short sale, have been considered. The obtained results have been compared
with the existing literature and the proposed technique is found to be optimum and better
in terms of profit.
Markov map is one example of interval maps where it is a piecewise expanding
map and obeys the Markov property. One well-known example of Markov map is the
doubling map, a map which has two subintervals with equal partitions. In this paper, we
are interested to investigate another type of Markov map, the so-called skewed doubling
map. This map is a more generalized map than the doubling map. Thus, the aims of this
paper are to find the fixed points as well as the periodic points for the skewed doubling
map and to investigate the sensitive dependence on initial conditions of this map. The
method considered here is the cobweb diagram. Numerical results suggest that there exist
dense of periodic orbits for this map. The sensitivity of this map to initial conditions is
also verified where small differences in initial conditions give different behaviour of the
orbits in the map.
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.
A mechanistic model has been used to explain the effect of radiation. The
model consists of parameters which represent the biological process following ionizing
radiation. The parameters in the model are estimated using local and global optimiza-
tion algorithms. The aim of this study is to compare the efficiency between local and
global optimization method, which is Pattern Search and Genetic Algorithm respectively.
Experimental data from the cell survival of irradiated HeLa cell line is used to find the
minimum value of the sum of squared error (SSE) between experimental data and sim-
ulation data from the model. The performance of both methods are compared based on
the computational time and the value of the objective function, SSE. The optimization
process is carried out by using the built-in function in MATLAB software. The parameter
estimation results show that genetic algorithm is more superior than pattern search for
this problem.
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.