nssvie
A python package for computing a numerical solution of stochastic Volterra integral equations of the second kind
where
- is an unknown process,
- is a continuous function,
- are continuous and square integrable functions,
- is the Brownian motion (see Wiener process) and
- is the Itô-integral (see Itô calculus)
by a stochastic operational matrix based on block pulse functions as suggested in Maleknejad et. al (2012) [1].
nssvie
is distributed under the terms of the GNU GPLv3 license.
Install
Install using either of the following two methods.
1. Install from PyPi
The nssvie
package is available on PyPi and can be installed using pip
pip install nssvie
2. Install from Source
Install directly from the source code by
git clone https://github.com/dsagolla/nssvie.git
cd nssvie
pip install .
Dependencies
nssvie
uses
- NumPy for many calculations,
- SciPy for computing the block pulse coefficients and
- stochastic for sampling the Brownian Motion
Usage
Consider the following example of a stochastic Volterra integral equation
so
>>> from nssvie import StochasticVolterraIntegralEquations
>>> # Define the function and the kernels of the stochastic Volterra
>>> # integral equation
>>> def f(t):
>>> return 1.0
>>> def k1(s,t):
>>> return s**2
>>> def k2(s,t):
>>> return s
>>> # Generate the stochastic Volterra integral equation
>>> svie = StochasticVolterraIntegralEquations(
>>> f=f, kernel_2=k1, kernel_1=k2, T=0.5
>>> )
>>> # Calculate numerical solution with m=20 intervals
>>> svie_solution = svie.solve_method(m=20, solve_method="bpf")
The parameters are
-
f
: the function . -
kernel_1
,kernel_2
: the kernels . -
T
: the right hand side of [0,T). Default is1.0
. -
m
: the number of intervals to divide [0,T). Default is50
. -
solve_method
: the choosen method based on orthogonal functions. Default isbpf
.
for the stochastic Volterra integral equation above.
Citation
[1] | Maleknejad, K., Khodabin, M., & Rostami, M. (2012). Numerical solution of stochastic Volterra integral equations by a stochastic operational matrix based on block pulse functions. Mathematical and computer Modelling, 55(3-4), 791-800. |