properGeneralizedInverseGaussianMean

auto expectation     = properGeneralizedInverseGaussianMean!double(2.0, 3.0, 4.0);
auto secondRowMoment = properGeneralizedInverseGaussianMean!double(2.0, 3.0, 4.0, 2);