Option pricing with the Black-Scholes model. Code not running

Question

(*initialisation of variables*)
n = 10;(*expected value greater than 3*)
r = 0.03;
\[Sigma] = 0.2;
sMax = 2;
sMin = -2;

(*define \[CapitalDelta]x.*)
\[CapitalDelta]x = (sMax - sMin)/(n - 1);

(*adding vector matrix... any assignment of v value should be 
included here*)
vecMatrix = Array[v, n];

(*reusable internal variables*)
p = 0.5 (\[Sigma]^2)/(\[CapitalDelta]x^2);
q = (2 r - (\[Sigma]^2))/\[CapitalDelta]x;

x = (2 p + r) - 0.75 q;
y = 5 p + q;
a = p - 0.25 q;
b = -2 p - r;
c = p + 0.25 q;

(*introducing matC for one time initialisation*)
matC = SparseArray[{Band[{1, 1}, {n, n}, n - 1] -> x, 
    Band[{1, 2}, {n, n}, {n - 1, n - 3}] -> y, 
    Band[{2, 1}, {n - 1, n}] -> a, Band[{2, 2}, {n - 1, n}] -> b, 
    Band[{2, 3}, {n - 1, n}] -> c}, {n, n}];
 MatrixForm [matC]*
 MatrixForm[vecMatrix]

Answer 1

I've done a lot of option pricing with MMA, but I can't relate your code to plain old Black-Scholes:

d1[spot_,strike_,ir_,div_,vol_, T_] = (Log[spot/strike] + (ir-div+vol^2/2) T) / (vol T^0.5]);  
d2[spot_,strike_,ir_,div_,vol_, T_] = (Log[spot/strike] + (ir-div-vol^2/2) T) / (vol T^0.5]);
N[z_] = (1 + Erf[z/Sqrt[2]])/2 ; (*cumulative normal density*)

bsCall[spot_,strike_,ir_,div_,vol_,T_] :=
   spot*E^(-div*T)*N[d1[spot,strike,ir,div,vol, T]] - 
   strike*E^(-r*T)*N[d2spot,strike,ir,div,vol,T]]

where

ir = interest rate  
div=dividends paid  
vol=volatility of the stock as measured by the standard deviation of its price
T = time to expiry  
Log[] = the natural logarithm

(I added this simple stuff to help others understand your question.) I don't see where the Black-Scholes part comes into your code? It would help educate others if you put your question in a context, for example, black-scholes option pricing is easy if you have one option, but if you have 100,000 at varying strikes, etc. But your question seems to really be about matrix algebra and sparce arrays.

Answer 2

You have a syntax error in the last line. I suspect you mean:

matC.vecMatrix // MatrixForm

which yields:

0.19875 v[1]+0.55125 v[2]
0.09 v[1]-0.2325 v[2]+0.1125 v[3]
0.09 v[2]-0.2325 v[3]+0.1125 v[4]
0.09 v[3]-0.2325 v[4]+0.1125 v[5]
0.09 v[4]-0.2325 v[5]+0.1125 v[6]
0.09 v[5]-0.2325 v[6]+0.1125 v[7]
0.09 v[6]-0.2325 v[7]+0.1125 v[8]
0.09 v[7]-0.2325 v[8]+0.1125 v[9]
0.09 v[8]-0.2325 v[9]+0.1125 v[10]
0.55125 v[9]+0.19875 v[10]

As an aside, comments that just say (* define \[Delta] x *) aren't that helpful to the next person looking at the code. I would recommend using commenting to explain what the variables are for.

As a further comment on your code, all your variables are global to your session, so the comment "reusable internal variables" is misleading. If you want to scope variables as internal to a particular calculation, you need to use scoping constructs like Module, Block and With. You might find this question helpful.