docs: public static RandomVariable bachelierOptionValue(
  			RandomVariable forward,
  			RandomVariable volatility,
  			double optionMaturity,
  			double optionStrike,
  			RandomVariable payoffUnit)
  	{
  		if(optionMaturity < 0) {
  			return forward.mult(0.0);
  		}
  		else
  		{
  			RandomVariable integratedVolatility = volatility.mult(Math.sqrt(optionMaturity));
  			RandomVariable dPlus	= forward.sub(optionStrike).div(integratedVolatility);
  
  			RandomVariable valueAnalytic = dPlus.apply(NormalDistribution::cumulativeDistribution).mult(forward.sub(optionStrike))
  					.add(dPlus.apply(NormalDistribution::density).mult(integratedVolatility)).mult(payoffUnit);
  
  			return valueAnalytic;
  		}
  	}