Page 1 of 21
Journal for Studies in Management and Planning
Available at http://internationaljournalofresearch.org/index.php/JSMaP
e-ISSN: 2395-0463
Volume 01 Issue 07
August 2015
Available online: http://internationaljournalofresearch.org/ P a g e | 126
N-dimensional stochastic Runge-Kutta scheme for Stochastic
Differential Equation Model
1E.O.Ogbaji, 2E.S.Onah, 3T.Aboiyar and 4A.R.Kimbir
1Department of Mathematics and Statistics, Federal University Wukari
(Ogbajieka@yahoo.com)
2,3,4Department of Mathematics/Statistics/Computer Science University Of Agriculture, Makurdi
Abstract
We formulated a four-compartment stochastic
differential equation which follow geometric
Brownian motion model. We showed the
existence and uniqueness of the solution
.Stochastic Runge-Kutta scheme was used to
solve the model.N- dimension stochastic Rung- Kutta scheme ,existence and uniqueness of
solution of stochastic differential equation was
established.
Key Words: geometric Brownian motion
model, stochastic differential equation,
Stochastic Runge-Kutta scheme
1.0 Introduction
The concept of Stochastic Differential
Equation (SDE) has been initiated by Einstein
in 1905, [56]. In his article he presented a
mathematical connection between microscopic
random motion of particles and the
macroscopic diffusion equation. Later it has
been seen that the stochastic differential
equation (SDE) model plays a prominent role
in a range of application areas such as physics,
chemistry, mechanics, biology,
microelectronics, economics and finance.
Earlier the SDEs were solved by using Ito
integral as an exact method which is discussed
in, [41]. But using exact method it is noticed
that there occur some difficulty to study
nontrivial problems and hence approximation
methods are used. In this context various
authors have given their contribution in these
field but we have mentioned which are directly
related to this problem. In [55] defined general
Runge-Kutta approximations for the solution
of stochastic differential equations and there
was given an explicit form of the correction
term. This work was carried out and then [37]
are discussed about the numerical solutions of
stochastic differential equation in detail. [49]
added discrete time strong and weak
approximation methods for the numerical
methods to find the solution of stochastic
differential equations. Next, [29] gave a major
contribution in this field to solve the
approximate solutions of stochastic differential
Page 2 of 21
Journal for Studies in Management and Planning
Available at http://internationaljournalofresearch.org/index.php/JSMaP
e-ISSN: 2395-0463
Volume 01 Issue 07
August 2015
Available online: http://internationaljournalofresearch.org/ P a g e | 127
equations and discussed few problems. Further
[30] investigated nonlinear stochastic
differential equations numerically. They
presented two implicit methods for Ito
stochastic differential equations (SDEs) with
Poisson-driven jumps. The first method is a
split-step extension of the backward Euler
method and the second method arises from the
introduction of a compensated, martingale,
form of the Poisson process. In this context
different authors have tried for various other
diffusion and application based problems. [28]
solved stochastic point kinetic reactor
problem. They modelled the point stochastic
reactor problem into ordinary time dependent
stochastic differential equation and studied the
stochastic behaviour of the neutron flux.
2.0 Formulation of Stochastic Differential
Equation Model
We formulate a four-compartment stochastic
differential equation for stock markets model.
Stochastic differential equation consists of
deterministic differential equations and white
noise term ‘dwt’. The white noise term is
stochastic in nature. Our formulation follow
geometric Brownian motion model which is
one form of Ito stochastic differential
equation. We use Figure1 Schematic diagram
of stock market interaction for the formulation.
Definition of Variables and Parameters
St Stock price
Rt
Return on investment
Vt
Stock price volatility
Ut Return on investment volatility
μ Rate of return on investment with reinvestment,
X Mean reversion of Vt,
Y Volatility of Vt,
Page 3 of 21
Journal for Studies in Management and Planning
Available at http://internationaljournalofresearch.org/index.php/JSMaP
e-ISSN: 2395-0463
Volume 01 Issue 07
August 2015
Available online: http://internationaljournalofresearch.org/ P a g e | 128
J Mean reversion of Ut,
N Rate of return on investment without reinvestment
Y
X
N
J
Figure1 Schematic diagram of stock market interaction
The stochastic differential equation model with respect to the stock price compartment is given as,
dS
t
dt
=
2
1
St +
2
1
S
t Vt
ξ(t)
where
0
Vt
St
Rt
Ut
