public function test_gammaln() { $this->assertEquals(round(log(1), 3), round(Stats::gammaln(1), 3)); $this->assertEquals(round(log(1), 3), round(Stats::gammaln(2), 3)); $this->assertEquals(round(log(1.3293326), 5), round(Stats::gammaln(2.5), 5)); $this->assertEquals(round(log(2), 5), round(Stats::gammaln(3), 5)); $this->assertEquals(round(log(6), 5), round(Stats::gammaln(4), 5)); $this->assertEquals(round(log(24), 5), round(Stats::gammaln(5), 5)); $this->assertEquals(round(log(120), 4), round(Stats::gammaln(6), 4)); }