RHMC sampling error

[William]

When I try to use the RHMC algorithm for flux sampling, I receive an error. Below is my code and the output:

Code

[polytope, flux_samples] =  sampleCbModel(casM_kidney_12_week_model, [], 'RHMC', options);

Output

Unrecognized property 'extraHessian' for class 'TwoSidedBarrier'.

Error in Polytope (line 111)
o.barrier.extraHessian = opts.extraHessian;

Error in sample (line 69)
polytope = Polytope(problem, opts);

Error in sampleCbModel (line 309)
o = sample(P, nPointsReturned, opts);