@Override
    public double variance() {
        return k * theta * theta;
    }
@Test
    public void testVariance() {
        System.out.println("variance");
        GammaDistribution instance = new GammaDistribution(3, 2.1);
        instance.rand();
        assertEquals(13.23, instance.variance(), 1E-7);
    }