<?php
error_reporting(-1);
 function to32bit($t) {
 return(gmp_and($t,"0xffffffff"));
 }
 
 function gmp_shiftl($y,$n) {
 return(gmp_mul($y,gmp_pow(2,$n)));
 }
 
 function gmp_shiftr($y,$n) {
 return(gmp_div($y,gmp_pow(2,$n)));
 }
 
 
$MathExpectation = gmp_init(0);
$Dispersion = gmp_init(0);
$AverageSquareDifference = gmp_init(0);
$P = 10;
 
for ($i=1 ; $i <= $P; $i++) {
$x = gmp_init(mt_rand(120000000,123456789));
$y = gmp_init(mt_rand(360000000,362436000));
$z = gmp_init(mt_rand(510000000,521288629));
$c = gmp_init(mt_rand(7600000,7654321));
 
 
$x = to32bit(gmp_add(gmp_mul("69069",$x),"12345")); // Линейный конгруэнтный метод Xn+1 = (aXn +c)
$y = to32bit(gmp_xor($y,gmp_shiftl($y,13)));      // Xor-shift начало
$y = to32bit(gmp_xor($y,gmp_shiftr($y,17)));
$y = to32bit(gmp_xor($y,gmp_shiftl($y,5)));       // Xor-shift конец
 
$t = to32bit(gmp_add(gmp_mul("698769069",$z),$c)); 
$c = to32bit(gmp_shiftr($t,32));                  
$z = $t;                                         
 
$result = to32bit(gmp_add(gmp_add($x,$y),$z));
 
$result = gmp_init($M);
 
$MathExpectation = gmp_add($MathExpectation,$M);
 
echo $i . "|" . gmp_strval($result) . "<br></br>\n";
                   };             		   
 
$MathExpectation = gmp_div($MathExpectation,$P);
 
echo "<br> Математическое ожидание: {$MathExpectation} </br>";
echo "<br> Среднеквадратическое отклонение: {$AverageSquareDifference} </br>";
echo "<br> Дисперсия: {$Dispersion} </br>";
?>