How to create a Math Function raster dataset

'Create the Raster Function object and Function Arguments object for the first operation.
Dim rasterFunction1 As IRasterFunction
rasterFunction1 = New MathFunction()
Dim mathFunctionArguments1 As IMathFunctionArguments
mathFunctionArguments1 = New MathFunctionArguments()

'Open raster datasets and specify as inputs for the "Plus" operation.
mathFunctionArguments1.Operation = esriGeoAnalysisFunctionEnum.esriGeoAnalysisFunctionPlus
mathFunctionArguments1.Raster = inRasterDataset01
mathFunctionArguments1.Raster2 = inRasterDataset02

'Create and initialize first function raster dataset with the Raster Function object and its arguments object.
Dim functionRasterDataset1 As IFunctionRasterDataset
functionRasterDataset1 = New FunctionRasterDataset()
functionRasterDataset1.Init(rasterFunction1, mathFunctionArguments1)

'Create the Raster Function and the Function Arguments object for the second operation.
Dim rasterFunction2 As IRasterFunction
rasterFunction2 = New MathFunction()
Dim mathFunctionArguments2 As IMathFunctionArguments
mathFunctionArguments2 = New MathFunctionArguments()

'Specify input rasters for the "Divide" operation.
'Use the output function raster dataset from the 1st operation as one of the input.
mathFunctionArguments2.Operation = esriGeoAnalysisFunctionEnum.esriGeoAnalysisFunctionDivide
mathFunctionArguments2.Raster = functionRasterDataset1
mathFunctionArguments2.Raster2 = inRasterDataset03

'Create and initialize the second function raster dataset.
Dim functionRasterDataset2 As IFunctionRasterDataset
functionRasterDataset2 = New FunctionRasterDataset()
functionRasterDataset2.Init(rasterFunction2, mathFunctionArguments2)

//Create the Raster Function object and Function Arguments object for the first operation. 
IRasterFunction rasterFunction1 = new MathFunction();
IMathFunctionArguments mathFunctionArguments1 = new MathFunctionArguments()as
    IMathFunctionArguments;

//Open raster datasets and specify as inputs for the "Plus" operation.
mathFunctionArguments1.Operation =
    esriGeoAnalysisFunctionEnum.esriGeoAnalysisFunctionPlus;
mathFunctionArguments1.Raster = inRasterDataset01;
mathFunctionArguments1.Raster2 = inRasterDataset02;

//Create and initialize first function raster dataset with the Raster Function object and its arguments object.
IFunctionRasterDataset functionRasterDataset1;
functionRasterDataset1 = new FunctionRasterDataset();
functionRasterDataset1.Init(rasterFunction1, mathFunctionArguments1);

//Create the Raster Function and the Function Arguments object for the second operation.
IRasterFunction rasterFunction2 = new MathFunction();
IMathFunctionArguments mathFunctionArguments2 = new MathFunctionArguments()as
    IMathFunctionArguments;

//Specify input rasters for the "Divide" operation.
//Use the output function raster dataset from the first operation as one of the input. 
mathFunctionArguments2.Operation =
    esriGeoAnalysisFunctionEnum.esriGeoAnalysisFunctionDivide;
mathFunctionArguments2.Raster = functionRasterDataset1;
mathFunctionArguments2.Raster2 = inRasterDataset03;

//Create and initialize the second function raster dataset.
IFunctionRasterDataset functionRasterDataset2;
functionRasterDataset2 = new FunctionRasterDataset();
functionRasterDataset2.Init(rasterFunction2, mathFunctionArguments2);