public void C9_9_4_European()
{
var opt = _testOption;
var k = opt.T / numberOfSteps;
var discounting = Math.Exp(-opt.r * k);
var binParams = new CRRStrategy(opt.sig, opt.r, k);; // Factory
var bn = new BinomialMethod(discounting, binParams, numberOfSteps);
bn.ModifyLattice(opt.K);
// Phase III: Backward Induction and compute option price
Vector <double> RHS = new Vector <double>(bn.BasePyramidVector());
if (binParams.BinomialType == BinomialType.Additive)
{
RHS[RHS.MinIndex] = opt.K * Math.Exp(numberOfSteps * binParams.DownValue);
for (int j = RHS.MinIndex + 1; j <= RHS.MaxIndex; j++)
{
RHS[j] = RHS[j - 1] * Math.Exp(binParams.UpValue - binParams.DownValue);
}
}
var pay = opt.PayoffVector(RHS);
var pr = bn.GetPrice(pay);
Assert.AreEqual(pr, 2.8307, 0.01);
}
public void C9_9_ExcelOutput()
{
var opt = _testOption;
var steps = 50;
var S = opt.K;
double k = opt.T / steps;
double discounting = Math.Exp(-opt.r * k);
var binParams = new CRRStrategy(opt.sig, opt.r, k);; // Factory
var bn = new BinomialMethod(discounting, binParams, steps, opt.EarlyImpl);
bn.ModifyLattice(opt.K);
// Phase III: Backward Induction and compute option price
Vector <double> RHS = new Vector <double>(bn.BasePyramidVector());
if (binParams.BinomialType == BinomialType.Additive)
{
RHS[RHS.MinIndex] = S * Math.Exp(steps * binParams.DownValue);
for (int j = RHS.MinIndex + 1; j <= RHS.MaxIndex; j++)
{
RHS[j] = RHS[j - 1] * Math.Exp(binParams.UpValue - binParams.DownValue);
}
}
var pay = opt.PayoffVector(RHS);
var pr = bn.GetPrice(pay);
// Display lattice in Excel
var file = Path.GetTempFileName();
file = Path.ChangeExtension(file, "xlsx");
ExcelMechanisms exl = new ExcelMechanisms(file);
try
{
// Display in Excel; first create array of asset mesh points
int startIndex = 0;
Vector <double> xarr = new Vector <double>(steps + 1, startIndex);
xarr[xarr.MinIndex] = 0.0;
for (int j = xarr.MinIndex + 1; j <= xarr.MaxIndex; j++)
{
xarr[j] = xarr[j - 1] + k;
}
string sheetName = "Lattice";
exl.printLatticeInExcel(bn.GetOptionLattice, xarr, sheetName);
exl.Save();
}
catch (Exception e)
{
Console.WriteLine(e);
}
Assert.AreEqual(pr, 3.0732, 0.01);
}
