%% Use SABR-BraceGatarekMusiela Model (Rebonato Parametric Form with Null Forward-Volatility Correlation) and IRMonteCarlo Pricer to Price Fixed Bond Instrument
% This example shows the workflow to price a FixedBond instrument when you use a SABRBraceGatarekMusiela model in Rebonato parametric form with null Forward-Volatility correlation and an IRMonteCarlo pricing method.
%% Create FixedBond Instrument Object
% Use fininstrument to create a FixedBond instrument object.

FixB = fininstrument("FixedBond",'Maturity',datetime(2021,7,1),'CouponRate',0.021,'Period',2,'Basis',1,'Principal',100,'Name',"fixed_bond_instrument")

%
%% Create Rebonato Form with Null Forward-Volatility Correlation SABRBraceGatarekMusiela Model Object
% Use finmodel to create a SABRBraceGatarekMusiela model object that is a SABR-BGM model in Rebonato parametric form with null Forward-Volatility correlation.

Alpha = [0.4;0.34;0.31;0.28];
Beta = [0.5;0.5;0.5;0.5];
numRates = 5;

SABRBGMVolVolFunc = @(a,t) (a(1)*t + a(2)).*exp(-a(3)*t) + a(4);
SABRBGMVolVolParams = [.3 -.02 .7 .14];

SABRBGMVolFunc = @(a,t) (a(1)*t + a(2)).*exp(-a(3)*t) + a(4);
SABRBGMVolParams = [.2 -.01 .8 .16];
VolFunc(1:numRates-1,1) = {@(t) SABRBGMVolFunc(SABRBGMVolParams,t)};

VolVolFunc(1:numRates-1,1) = {@(t) SABRBGMVolVolFunc(SABRBGMVolVolParams,t)};

CorrFunc = @(i,j,B) exp(-B*abs(i-j));
FwdFwdCorrelation = CorrFunc(meshgrid(1:numRates-1)',meshgrid(1:numRates-1),.08);
VolVolCorrelation = CorrFunc(meshgrid(1:numRates-1)',meshgrid(1:numRates-1),.04);

SABRBGM_Rebonato = finmodel("SABRBraceGatarekMusiela",'Alpha',Alpha,'Beta',Beta,'VolatilityofVolatility',VolVolFunc,'Volatility', VolFunc,'FwdFwdCorrelation', FwdFwdCorrelation, 'VolVolCorrelation', VolVolCorrelation)

%% Create ratecurve Object
% Create a ratecurve object using ratecurve.

Settle = datetime(2019,1,1);
Type = 'zero';
ZeroTimes = [calmonths(6) calyears([1 2 3 4 5 7 10 20 30])]';
ZeroRates = [0.0052 0.0055 0.0061 0.0073 0.0094 0.0119 0.0168 0.0222 0.0293 0.0307]';
ZeroDates = Settle + ZeroTimes;
 
myRC = ratecurve('zero',Settle,ZeroDates,ZeroRates)

%% Create IRMonteCarlo Pricer Object
% Use finpricer to create an IRMonteCarlo pricer object and use the ratecurve object for the 'DiscountCurve' name-value pair argument.

simDates = datetime(2019,7,1)+calmonths(0:6:24);
outPricer = finpricer("IRMonteCarlo",'Model',SABRBGM_Rebonato,'DiscountCurve',myRC,'SimulationDates',simDates)

%% Price FixedBond Instrument
% Use price to compute the price and sensitivities for the FixedBond instrument.

[Price,outPR] = price(outPricer,FixB,["all"])
outPR.Results

% You can access the simulated interest-rate Paths in the PricerData output.

outPR.PricerData

% Copyright 2021 The MathWorks, Inc.